Nonlocal magnetization dynamics in ferromagnetic heterostructures 

Yaroslav Tserkovnyak/ Arne Brataas,^ Gerrit E. W. Bauer, ^ and Bertrand I. Halperin^ 

^ Lyman Laboratory of Physics, Harvard University, Cambridge, Massacliusetts 02138, USA 

^ Department of Physics, Norwegian University of Science and Technology, N-7491 Trondheim, Norway 

^Kavli Institute of NanoScience, Delft University of Technology, 2628 CJ Delft, The Netherlands 

''Lyman Laboratory of Physics, Harvard University, Cambridge, Massachusetts 02138, USA 

(Dated: February 2, 2008) 

Two complementary effects modify the GHz magnetization dynamics of nanoscale heterostructures 
of ferromagnetic and normal materials relative to those of the isolated magnetic constituents: On 
the one hand, a time-dependent ferromagnetic magnetization pumps a spin angular-momentum 
flow into adjacent materials and, on the other hand, spin angular momentum is transferred between 
ferromagnets by an applied bias, causing mutual torques on the magnetizations. These phenomena 
are manifestly nonlocal: they are governed by the entire spin-coherent region that is limited in size 
by spin-flip relaxation processes. We review recent progress in understanding the magnetization 
dynamics in ferromagnetic heterostructures from first principles, focusing on the role of spin 
pumping in layered structures. The main body of the theory is semiclassical and based on a 
mean-field Stoner or spin-density— functional picture, but quantum-size effects and the role of 
electron-electron correlations are also discussed. A growing number of experiments support the 
theoretical predictions. The formalism should be useful to understand the physics and to engineer 
the characteristics of small devices such as magnetic random-access memory elements. 
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A ferromagnet is a symmetry-broken state in which 
a majority of electrons point their spin into a certain 
common direction below critical temperatures as high as 

I III. Spin emission bv coherentlv-precessing ferromagnet£ll4 IQOO K. The robustness of the magnetic Order and the 

permanence of a given magnetization direction against 
elevated temperatures and external perturbations have 
been employed in applications as diverse as compass nee- 
dles, refrigerator-door stickers, and memory devices. 

In spite of its stability, ferromagnetism is neither rigid 
nor static. Depending on sample size and anisotropics 
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ferromagnet is often unstable with respect to a domain 
structure that lowers the macroscopic magnetic energy. 
Thermal fluctuations reduce the macroscopic moment 
until it completely vanishes at the critical temperature 
Tc. At temperatures sufficiently below Tc, the inter- 
nal dynamics of the ferromagnet are dominated by low- 
energy transverse fluctuations of the magnetization, so- 
called spin waves or magnons, that are the magnetic 
equivalence of phonons in a lattice. Classical course-grain 
computer simulations of the detailed position- and time- 
dependent magnetizatio n ( "micromagneti s m" ) describe 
these phenomena well (jBrown. Jrl 119631: iMiltat et all 
120021) . 

When magnetic grains become sufficiently small, the 
exchange stiffness renders domain structures energeti- 
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cally unfavorable and a single-domain picture is ade- 
quate. When the ferromagnet is exposed to a uniform 
driving field, the macroscopic magnetization dynamics 
may then be dominated by a collective precession of the 
entire ferromagnetic order parameter. Restricting the 
ferromagnetic degrees of freedom to this mode is often re- 
ferred to as the "macrospin model." In infinite ferromag- 
netic media, low-energy spin waves resemble symmetry- 
restoring Goldstone modes, but in real life, the spin- 
rotational symmetry is broken by magnetic anisotropics 
caused by magnetic dipolar fields or by crystal-field spin- 
orbit interactions. In thermodynamic equilibrium, the 
macrospin then points in a certain fixed direction with 
small thermal fluctuations around it. The ferromagnet 
can still be coerced into motion by applying external 
magnetic fields at a finite angle to the magnetization 
direction. The system then moves in response, trying 
to minimize its Zeeman energy. The compass needle, a 
freely suspended single-domain ferromagnet with a suf- 
ficiently high anisotropy (coercivity), does this by align- 
ment of its lattice. In this review, we are interested 
in mechanically-fixed magnets whose magnetic moments 
move in the presence of external and anisotropy effec- 
tive magnetic fields, as well as applied electric currents. 
Viscous damping processes are required to achieve a re- 
orientation (switching) of the magnetization, if, for ex- 
ample, a magnetic-field direction is suddenly changed. 
Minimization of this finite switching time by engineering 
magnetic anisotropics and magnetization-damping rates 
is an important goal in the design of fast magnetic mem- 
ories. When the applied magnetic fields are large enough 
to surmount the anisotropics, the magnetization can be 
reversed, often by large amplitude and complex trajec- 
tories, even in the simple macrospin model. At finite 
temperatures, the magnetization reorientation becomes 
probabilistic and is des cribed by a Fokke r-Planck equa- 
tion on the unit sphere l)Brown. Jrl Il963|) . 

In the last two decades, a new subdiscipline in the field 
of magnetism has risen that is devoted to the studies of 
heterostructures of ferromagnets (F) with normal metals 
(N) and, to a lesser extent, semiconductors and supercon- 
ductors. Especially magnetoelectronics, the science and 
technology directed at understanding and utilizing the 
transport properties of layered structures of ferromag- 
netic and normal metals, has grown into a mainstream 
topic of condensed-matter physics. Its attraction derives 
from large effects at room temperature that can be un- 
derstood easily in terms of transparent physics and that 
have found already numerous applications. Two crucial 
discoveries in magnetic multilayers still reverberate in re- 
cent res earch, viz., the nonlocal oscillatory exchange cou- 
pling bv lGriinberg e^'aT. '1986' and the giant magnetore- 
sistance (GMR) by Baibich et ai, 1988; Binasch et iil.. 
Il989i The exchange coupling through a metallic spacer 
favors an antiparallel coupling between ferromagnetic 
layers for certain spacer thicknesses, depending on the 
occupation of spin-polarized quantum-well states. It is 
therefore a quantum-interference effect sensitive to de- 



fect scattering, which vanishes exponentially with in- 
creasing spacer-layer thickness. GMR stands for a sig- 
nificant reduction of the resistance of multilayers when 
the magnetic configuration is forced from antiparallel 
to parallel by an applied magnetic field. In disor- 
dered multilayers, it is a semiclassical transport effect 
that can be understood in terms of a diffusion equation 
(Camlev and Barnas. 1989 ^ Valet and Pert, 1993). In a 
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transport can be mapped on a two-channel (spin-up and 
spin-down) resistor model, in which interface and bulk 
resistances for a fixed spin are simply connected in se- 
ries. The spin-relaxation processes are usually modeled 
by a local finite-resistance link connecting the spin-up 
and spin-down circuits. 

Initially, the community focused its attention on sta- 
tionary magnetic states, like those responsible for the 
magnetoresistance in metallic and tunneling structures 
with applied dc bias. This has changed drastically in re- 
cent years. The main catalyst was the experimental v eri- 
fications of an earlier prediction bv lSlonczews and 
lBergeilll99fil that electric currents can cause a reorienta- 
tion of the ferrom agnetic order in multilayer structures. 
iTsoi et "aHi . Il998l experimentally demonstrated magneti- 
zation precession in (Go|Gu)Ar multi layers by currents 
injected by a point contact, whereas iMvers et at, 199^ 
observed switching in the orientation of magnetic mo- 
ments in Go|Cu|Go sandwich structures by perpendicu- 
lar electric currents ("GPP spin valves"). Much earlier, 
a coupling between a dynamic ferromagnetic magneti- 
zation and spin accumulatio n in adjacent norma l met- 
als has been postu lated by iJanossv and MonodL Il976t 
ISilsbee et all Il979l These authors demonstrated that 
microwave transmission through normal-metal foils facil- 
itated by conduction-electron spin transfer is significantly 
enhanced by ferromagnetic-layer coating. 

This review covers the developments in the under- 
standing of the magnetization dynamics in heterostruc- 
tures of ferromagnets and normal conductors in the last 
five years or so. We believe that the time is ripe, since 
from a microscopic point of view, much of the basic 
physics is well understood. A consistent and coherent 
picture has evolved that is based on the diffusion equa- 
tion for the bulk transport in metallic ferromagnets and 
normal conductors with quantum-mechanical boundary 
conditions at possibly sharp interfaces between them. 
NoncoUinearity of the magnetization directions in struc- 
tures with more than one magnet is an essential ingre- 
dient in order to understand the physics. We focus here 
on explicitly dynamic effects, referring to a separate arti- 
cle (Brataas et ai, 2 005a,) for the static transport prop- 
erties of magnetoelectronic circuits and devices. The 
present review follows a self-consistent approach extend- 
ing the static magnetoelectronic circuit theory to time- 
dependent phenomena. It provides a framework to in- 
clude on equal footing two physical effects that are two 
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sides of one coin, viz. , the spin-tr a nsfer torque induced 
by appHed currents ijSlonczewskl Il996() and the spin 
pum ping by moving ferromagne ts into adjacent conduc- 
tors (iTserkovnvak et oil 1200231) . The theory is derived 
from microscopic principles and the material-dependent 
input parameters are thus accessible to ah initio calcula- 
tions. We concentrate on quasi-one-dimensional models 
corresponding to, e.g., layered pillar structures. With 
few exceptions, we do not attempt accurate modeling of 
concrete device structures and deviations of the magneti- 
zation dynamics from the macrospin model, although the 
theory discussed in this review can be readily generalized 
to treat such situations. We focus on "adiabatic" effects 
to lowest order in the characteristic Larmor frequency. 
In spite of these limitations, the agreement with various 
experiments is found to be gratifying. 

Effects beyond such a model certainly may cause ob- 
servable phenomena. For example, the quantum inter- 
ference that leads to inversion of the magnetoresistance 
in high-quality tunnel junctions cannot b e treated by 
the diffusion equation ijYuasa et al\ . Nonlinear- 
ities require numerical simulations or a stability anal- 
ysis based on the t heory of dyn amic systems that are 
outside our scope ijValetl l2004|) . High temperatures 
and currents can best be treate d by stochastic meth- 
ods beyond the present re view ijApalkov and Visscheii 
120041: iLi and Zhand. l2004b[) . but the input parameters of 
such approaches are provid ed here. The current-induced 
dynamics of domain walls ( Barnes and Mackawql l2005t 
iLi and ZhailS l2004at pTat ara and Kohno, 2004 and ref- 
erences therein) are also beyond the macrospin consider- 
ations central to this review. Whereas the spin-transfer- 
torque-induced dynamics are a crucial ingredient, space 
constraints force us to abandon complete coverage of the 
numerous experiments published recently. 

Throughout the review, we focus on self-consistent ef- 
fects arising from the time-dependent ferromagnetic ex- 
change field felt by itinerant carriers in the mean-field 
picture. We take the spin-orbit interaction into ac- 
count only in terms of a phenomenological spin-flip re- 
laxation time. Sees. IIIIVl but consider it more seriously 
in Sees. IVI.BI and IVII.BI Most results are not mate- 
rial specific, but, unless specified otherwise, we have het- 
erostructures of transition-metal ferromagnets (and its 
alloys) with noble or other simple normal metals in mind. 

The main body of this review is organized as fol- 
lows: Sec. U introduces several basic concepts that we 
rely on in the remainder of the article, which is pri- 
marily aimed at nonspecialists. Sec. ^] is a brief but 
in-depth discussion of the magneto electronic circuit the- 
ory (see also Brataas et alV f2005aj) . which is then gen- 
eralized in Sec. IIIII to dynamic problems by means of 
the spin-pumping concept. Sees. HVl and Ivl respectivelv 
discuss Gilbert damping and dynamic ferromagnetic ex- 
change in heterostructures, which are mediated by spin 
pumping and spin-transfer torques. Sec. I VII is devoted to 
an alternative linear-response view of the nonlocal mag- 
netization dynamics, and Sec. IVIll treats several special 



topics before we conclude the article with summary and 
outlook in Sec. IVTTT1 

B. Nonlocal exchange coupling and giant 
magnetoresistance 

The discovery that the energy of magnetic multilayers 
made from alternating ferromagnetic and normal-metal 
films depends on t he relative direction o f the individ- 
ual magnetizations ijGriinberg et aZ.Lri986() is perhaps the 
most important in magnetoelectronics. The existence 
of the antiparallel ground-state configuration at certain 
spacer-layer thicknesses was essential for the subsequent 
discovery of the giant magnetoresistance ijBaibich et alV 
1988; Binasch et al., 1989J. Adjacent ferromagnetic lay- 
ers in such structures are coupled by nonlocal and, as 
a function of normal- metal layer thickness, oscillatory 
ijParkin et al\ . Il990j) exchange interaction that can be 
qualitatively under stood by perturbation theory analo- 
gous to the RKKY li|kasuvaKll956HRuderman and Kittel 
11954 F^sidalll95'^ exchange coupling between magnetic 
impurities in a normal-metal host. Different oscillation 
periods, which can be resolved in measured magnetiza- 
tion configuration as a function of spacer thickness, are 
well explained in terms of the normal- metal Fermi surface 
calipers in the growth direction. The magnetic ground- 
state configuration is, at least in principle, accessible 
to first-principles electronic-structure calculations in the 
spin-density-functional theory formalism, and that is ba- 
sically the end of the story. However, in order to make 
connection to the main topic of this review, we briefly 
discuss the formulation of the equil ibrium exchange cou- 
pling in terms of scattering theory ijErickson et al I ll993t 
ISlonczewskL .1989. 1993), that can also be formulated 
from fir st princ i ples and calculated by density-functional 
theorv ijBrunol il995HStileslll999() . Another advantage 
of a scattering-theory formulation is that effects of disor- 
der can be understood employing the machinery of meso- 
scopic physics, such as random-matrix l)Beenakkeilll997f) 
or diagrammatic perturbation theory. 

Let us consider a layered N|F|N|F|N spin valve with an- 
gle between the magnetizations and an N-spacer with 
thickness L, see Fig. ^ schematic. Suppose we can view 
the F|N|F trilayer as some spin-dependent scatterer em- 
bedded into a normal-metal medium. The trilayer gives 
a 6'-dependent contribution to the tot al ground-state en- 
ergy, given by a standard formula ijAkkermans et all , 
[1993): 

E(L,e) = —[ e^lndets(L,6l,e)de, (1) 
2711 de 

in terms of the energy-dependent scattering matrix 
s{L,6,e) of the trilayer. The scattering matrix is made 
up from the matrices r and t of the reflection and trans- 
mission coefficients for a basis of spin-resolved incoming 
states at energy e from the right normal metal whereas 
the primed ones (r' and t') are defined for states coming 
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from the left normal metal: 



r t' 
t r' 



(2) 



The scattering matrix of the total system can be com- 
posed out of the transmission and reflection coefficients 
of the F|N interfaces as well as of the bulk layers by well- 
known concatenation rules. Quantum- well states and 
resonances are formed by multiple reflections at inter- 
faces caused by potential steps and electronic-structure 
mismatches. The angle and thickness dependence of the 
total energy |^ can be understood in terms of the vari- 
ation of the interference pattern of the spin-dependent 
electron waves in and close to the normal- metal spacer. 
By varying L and 0, quantum-well states enter or leave 
the Fermi sea with abrupt changes in the total energy 
that can be large for small L. The minimum energy for a 
given L is usually found at = and tt, i.e., either par- 
allel or antiparallel configurations are favored. Though 
exponentially suppressed, the coupling between the mag- 
netic la yers persists when the insertion is an insulating 
barrier ijBrunol Il995t ISlonczewskl Il989|) . 
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FIG. 1 Spin- valve schematic: two monodomain ferromagnets 
(F) separated by a normal (N) spacer and attached to nor- 
mal leads. In this CPP configuration, the current flowing be- 
tween two normal reservoirs sequentially traverses two mag- 
netic layers. Also shown are two effective circuits discussed 
in Sec. 111. AI for the semiclassical regime where the interlayer 
exchange coupling vanishes. 

Of special interest is the asymptotic dependence of en- 
ergy E{L, 9) for large L. In this limit, the sharp jumps in 
energy by the population or depopulation of individual 
quantum- well states as a function of L and 9 become less 
prominent. The quasi-continuous angular dependence of 
the energy in this limit is well described by the lowest 
term in the expansion into Legendre polynomials: 



EiL, 



= cos I 



i^^sin {qa±L 



which is a sum of contributions from each critical caliper 
at the Fermi surface of the normal-metal spacer labeled 
by a. The par ameter s Jn a nd c/jq, are model and mate- 
rial dependent ijStilesl . Il999() . qaj_ is the critical Fermi- 
surface-spanning wave vector in the layering direction, 
which determines the caliper in reciprocal space. Note 
that whereas we in principle require the scattering ma- 
trix for all occupied states in Eq. Q , Eq. Q is governed 
by the scattering coefficients at the Fermi energy only, 
just as the transport properties at low temperature and 
bias. In practice, Eq. © can often be used for all but 
the most narrow spacer layers. 

At configurations that are not at equilibrium, the 
derivative 



(4) 



(3) 



does not vanish. A finite r is therefore interpreted as 
an exchange torque acting on the magnetizations, pulling 
them into an energetically-favorable configuration. Phys- 
ically, this torque is an angular-momentum transfer that 
is carried by the electron spin. A spin valve that is 
"strained" by a relative misalignment of the magneti- 
zation directions from their lowest-energy configuration 
therefore supports dissipationless spin currents. The sit- 
uation is quite analogous to the Josephson junction in 
which a difference of the superconducting phase over a 
weak link induces a supercurrent. We note that at finite 
temperatures, the ground-state energy E{L, 9) should be 
replaced with the free energy F{L, 9) in Eq. 

Essential for the existence and the magnitude of the 
nonlocal exchange coupling and the corresponding per- 
sistent spin currents is the phase coherence of the wave 
functions in the normal spacer. An incoming electron 
in the spacer with information of the left magnetization 
direction has to be reflected at the right interface and in- 
terfere with itself at the left interface in order to convey 
the coupling information. This implies strong sensitivity 
to the effects of impurities, since the diffuse scattering 
destroys the regular interference pattern required by a 
sizable coupling. This qualitative notion has been for- 
mulated by Wainta l et aL , _2000 in the scattering-theory 
formalism invoking the "isotropy" condition for valid- 
ity of the random-matrix theory. Isotropy requires dif- 
fuse transport, viz., that L is larger than the mean free 
path due to bulk and interface scattering. It can then 
be shown rigorously that the equilibrium spin currents 
vanish on average with fluctuations that scale like N~^, 
where N stands for the number of transverse transport 
channels in the normal- metal spacer. In layered metallic 
structures, N is large and the static exchange coupling 
and spin currents can safely be disregarded in the diffuse 
limit. On top of the suppression by disorder, the magni- 
tude of the coupling scales like even in ballistic sam- 
ples, see Eq. Q. Experimentally, even the presumably 
best Co|Cu|Co samples indeed do not show an apprecia- 
ble coupling beyond a spacer-layer thickness of about 20 
atomic monolayers. 
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In section^ we discuss the magnetization dynamics of 
multilayers and superlattices. We find that, on top of the 
equilibrium spin currents that communicate the nonlocal 
static exchange coupling, a dynamic exchange interaction 
with a much longer range becomes important. Any signif- 
icant coupling at equilibrium in the dynamic studies can 
be represented approximately in terms of parametrized 
conservative forces that react to deviations from equilib- 
rium without interference with nonequilibrium spin cur- 
rents. 

The GMR has been originally discovered in a con- 
figuration in which the current flows in the plane of 
the film (CIP). Multilayer pillar structures in which 
the current flows perpendicular to the planes (GPP) 
llGiis and Baueilll997trGiis et a/.lll993tlPratt. Jr. all 
Il99l|) are more relevant in the present context. As- 
suming diffusive transport, the GPP GMR is easily un- 
derstood in terms of a two-channel series-resistor model 
l|Valet and Fertlll993(l : In the parallel configuration, the 
charge current is short-circuited by the low-electrical- 
resistance spin channel. The charge and spin transport in 
the intermediate configurations in which the magnetiza- 
tions vary between parallel and antiparallel is important 
in the context of the present review due to the emer- 
gence of transverse spin currents that are absorbed by 
the magnetizations and contribute as a driving torque 
to the dynamics. This dissipative channel for transverse 
spins also modifies the angular magnetoresistance of spin 
valves. These and other nonco Uinear magnetoelectro nic 
dc phenomena are reviewed by iBrataas et all l2005ai 



C. Landau-Lifshitz-Gilbert phenomenology 

At temperatures well below the ferromagnetic crit- 
ical temperature Tc, the equilibrium magnetization of 
a bulk ferromagnet saturates at some material-specific 
value Mg. Symmetry- restoring Goldstonc modes of this 
broken-symmetry state are spin waves (magnons) men- 
tioned in Sec. II. Al which can condense into macroscopic 
transverse magnetization dynamics. In itinerant ferro- 
magnets, there are also longitudinal spin excitations, the 
Stoner modes. However, since macroscopic variation of 
the magnetization magnitude at T <C Tc is very costly 
in free energy, we focus on the purely transverse mo- 
tion of the position-dependent magnetization direction 
m = M/Ms, with a fixed magnitude |M| = M^. 

A traditional starting point in studying the trans- 
verse magnetization dynamics in a ferromagnetic medium 
is based o n the phenomenolog ical Landau-Lifshitz (LL) 
equation l|Landaii et all Il98fl() . The magnetization di- 
rection m(r, t) is treated in this approach as a classical 
position- and time-dependent variable obeying equations 
of motion which are determined by the free-energy func- 
tional F[M] for degrees of freedom coupled to the magne- 
tization distribution M(r) (such as the electromagnetic 
field or itinerant electrons experiencing a ferromagnetic 



exchange field): 

5fm(r, t) = -7m(r, t) x Hoff(r) , (5) 

where 7 is (minus) the gyroniagnetic ratio and 

Heff(r) = ~dMF[M] (6) 

is the so-called "effective magnetic field." Gorresponding 
to the respective contributions to the free energy, the 
effective field can usually be decomposed into the ap- 
plied, dipolar demagnetization, crystal-anisotropy, and 
exchange fields. In the case of free electrons, 7 — 
2fiB/fi > 0, in terms of Bohr magneton hb and Planck's 
constant h = 2iTh, and it is usually close to this value in 
transition-metal ferromagnets. 

It is easy to see that the LL equation with ef- 
fective field ®, describes transverse- magnetization dy- 
namics preserving the free energy ^"[M]. Definition ® 
with the effective field depending on the instantaneous 
magnetic configuration assumes that the magnetization 
dynamics are very slow on the scale of the relevant micro- 
scopic relaxation processes. However, some slow degrees 
of freedom may not respond sufficiently fast to the mag- 
netization motion, making the effective field dependent 
on the history of the magnetization dynamics M(r, i). 
This should be associated with dissipation of energy into 
the degrees of freedom that are coupled to the magneti- 
zation. 

As a specific example, consider the magnetization dy- 
namics (O described by the effective field 

Heff(r,t) = -aM(i/(M)), , (7) 

where iJ(M) is the many-body Hamiltonian for itinerant 
electrons, parametrized by a mean-field magnetic config- 
uration M(r, t) of, e.g., some localized magnetic orbitals 
(as in the s — d model), and Qt evaluates its expec- 
tation value for the many-body state (or ensemble) at 
time t. Setting the many-body ensemble at time t to its 
thermal-equilibrium configuration determined by M(t) 
reproduces the Landau-Lifshitz (LL) definition In 
the opposite extreme, when electrons do not respond at 
all to the fast magnetic dynamics, ()t « ()o, the effective 
field is determined by the history-independent functional 
(i?(M))o instead. In the intermediate regime, a finite 
time lag in the response of the itinerant electrons to the 
varying magnetization causes dissipation of the magnetic 
energy, as discussed in Sec. IVI.BI To lowest order in fre- 
quency (i.e., keeping only terms linear in dt), such damp- 
ing can be described by an additional torque in Eq. ^ 
(lGilbertL[l955l 12004) : 

dtui = —7m X Hoff -I- am x dtui , (8) 

where a is the dimensionless Gilbert constant and Hcff 
is an effective field depending only on the instanta- 
neous magnetic configuration. (Partial time derivatives 
imply here the possibility of spatial variations of the 
magnetization, as, e.g., in the case of spin waves; the 
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full time derivatives are reserved for the monodomain 
dynamics.) The Gilbert term in Eq. JSJ has been 
obtained for various microscopic formulations of the 
magnetiz ation dynamic s, see, e. g., iH einrich et al!, '1967"; 
Korenm an and Prange l Il972t iKunes and K amberskv. 
2nn2|:lT.ntovinov a,nd T^eizeTlll 97flHRa,fonov a,nd Bertra,niL 
20001: ISinova et a/.ll200l iTserko vnvak e t a^J . l2002ai 

Energy dissipation implied by Eq. (|HJ preserves the lo- 
cal magnitude of the magnetization. For example, for a 
constant H^s obeying Eq. © and a = 0, m precesses 
around the field vector with frequency uj — 7i?cff- When 
damping is switched on, a > (assuming positive 7, as in 
the case of free electrons) , the precession spirals down on 
a time scale of (auj)^^ to a time- independent magnetiza- 
tion along the field direction, i.e., the lowest-free-energy 
state. Close to an equilibrium axis with rotational sym- 
metry, the Landau-Lifshitz-Gilbert (LLG) equation ^ 
is obeyed by a small-angle damped circular precession, 
while in the presence of anisotropics, small-angle trajec- 
tories are elliptic and the damping is in general a tensor. 
For most of our purposes, simple circular precession with 
a scalar damping a suffices (but see Sec. IVII!b1i . It is 
sometimes convenient to work with a different Gilbert 
parameter 



G = ajMs 



(9) 



It can be made explicit that magnetization dynamics 
described by Eq. |(SJ) dissipate energy at a rate determined 
by a. To this end, suppose for simplicity that 7Hoff = 
LUoZ, ujo > 0, is uniform throughout a monodomain fer- 
romagnetic sample, so that Eq. (jS)) describes a damped 
macrospin circular precession around the z axis. Small- 
angle dynamics around the z axis can thus be resonantly 
excited by a (right-hand) circularly polarized rf field with 
a small amplitude h^i and frequency uj close to wq, that is 
h-{t) = hx{t) — ihy{t) = /irf exp(— iwt). The magnetic re- 
sponse to such a field is 5M-{lo) = x v['^)h-{i^), where 



ll974HHeinrich and CochranLll993(l . FMR spectra are ob- 
tained by placing the sample into a microwave cavity and 
sweeping the external dc field. 7Hcff then determines 
the position and a the width of the resonance absorp- 
tion peak. The FMR lincwidth can have an additional 
contribution due to inhomogeneities in Hcs, loosely cor- 
responding to a finite range of the resonance frequency 
loq in Eq. Hll|l . For example, small disorder by sur- 
face roughness or a nonuniform surface field in exchange- 
biased thin films, contributes to the resonance broad- 
ening by ( in quantum-mechanical t erms) two-magnon 
scattering ("M ills and Rezend3 . l200^ . The inhomoge- 
neous linewidth broadening is associated with dephas- 
ing of the global precession that in general does not 
conserve the magnitude of the magnetization. Whereas 
the Gilbert damping predicts a strictly linear depen- 
dence of FMR linewidths on frequency, the inhomoge- 
neous broadening is usually associated with weaker fre- 
quency dependence as well as a zero-frequency contri- 
bution. Another common technique in studying long- 
wavelength spin wa ves is Brillouin light scattering (BLS) 
(see, e.g.. FPcmokrit ov and Tsvmb al. 1994). Both FMR 
and BLS probe magnetic excitations close to the surface, 
i.e., within the corresponding skin depth of the order of 
100 n m for FMR and 10 nm for BLS ( Mill s and Rezendd. 
'200^. In contrast to FMR, BLS excites spin waves 
with finite wavelengths in the surface plane (in the range 
of that of the visible lig ht), bearing consequences for 
the signal linewidths, see iMills and Rezend I I2OO3I and 
Sec. IVII.CI In closing this subsection, we remark that 
ferromagnetic magnetization dynamics and, in partic- 
ular, magnetization relaxation processes are collective 
many-body phenomena that continue to fascinate in 
spite of decades of theoretica l and experimental efforts 
to understand them, see , e.g.. lDobin and VictoraL 120031 
iQian and VignaleLl200l 



[uoq — uj) — iaui 



(10) 



is the transverse magnetic susceptibility. The linear- 
response expression for the energy-dissipation power per 
unit of volume 



P = t^Imx-+(c^)/irf 



(wo — + {auo)"" 



(11) 



does not depend on the microscopic origin of a, as long 
as Eq. © holds. For a steady precession, one can also 
show that P = — h(i) • rh{t)Ms is the work done by the 
rf field h(i). The stability of the system, P > 0, requires 
that 07 > 0. 

Eq. © has been found to very successfully de- 
scribe the dynamics of ultrathin ferromagnetic films 
as well as bulk materials in terms of a few material- 
specific parameters that are accessi ble to ferromagnetic- 
resonance (FMR) experiments l|Bhagat and LubitzL 



D. Current-induced magnetization dynamics 

It has been realized only relatively recently that in 
magnetic multilayers the magnets can be excited by other 
me ans than exter nal magnetic fields. ISlonczewskil Il996l 
and lBergeRll996l predicted that in CPP spin- valve struc- 
tures a dc current in the right direction can excite and 
even reverse the magnetization of a magnetic layer. This 
can be observed by monitoring dl /dV which depends on 
the magnetic configuration, as in the GMR. The pre- 
dictions have by no w been ampl y confi rmed by many 



dictions nave by no w been ampl y conn rmed by many 
recent experiments (iJi et all 12003,: Katine et all I200C 



Kiseley eii aZl.l2003tlKrivorotoy et a/.il2005HMvers et al. 

9nn9l I1QQC* n]\ lOnn^l IPnfoll t,^ nl\ 1900^ 



2002L 11999 11 'Ozvilmaz et alK '20031 'Puf all aZ.1,' 
Sun et qZj. 12005 : Tsoi et at. 2000: Urazh din 6^/]. 



2004; 



2003; 



Wegrowe et The current-induced magnetic 

dynamics have been found t o also affect current- 
noise spectra dCovingtoiL^ "aD. I2004ll. The predic- 
tion l|Polianski and Brouweit 120041: IStiles et all l2004|) 

that even a single dc-current-driven ferromagnetic 
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layer may undergo a resonant finite-wave-vector spin- 
wave excitation has b een expe rimentally confi rmed by 
lOzvilmaz etdl 120041 see also I.Ti et all hClOA Conse- 
quently, higher-wave-vector spin-wave e xcitations may 
in som e situations succe s sfully compete (iBrataas et all 
l2005hl : lOzvilmaz et all l2005^ with current-induced 
macrospin motion considered below. The current- 
induced magnetization dynamics pose a challenging 
physics problem that requires understanding of the cou- 
pling of nonequilibrium quasiparticles with the collective 
magnetization dynamics. It carries technological poten- 
tial as well. In small structures, the current-induced mag- 
netization reversal may be more power efficient to write 
information into magnetic random-access memories com- 
pared to switching by Ampere magnetic field. The gener- 
ation of microwaves by exciting stable magnetization or- 
bits with dc bias-de pendent frequencies may also satisfy 
technological needs l|Kiselev et aR l2003t iRippard et all 
12004'! . 

Current-induced magnetization dynamics are a con- 
sequence of the spin-dependent transport in F|N het- 
erostru ctures. For example , Slonczewski's magnetization 
torque l)SlonczewskiL Il996j) is equivalent to absorption of 
an incident spin current with a pol arization component 
perpendicular to the magn e tization JBr ataas et oiJ . l2000t 
IStiles and Zangwill 120021: IWaintal et all l200(]() . The 
component of the electron spin perpendicular to the mag- 
netization is not a constant of the motion in a ferromag- 
net. On the other hand, neglecting the effects of spin- 
orbit coupling (other than the macroscopic anisotropy 
already included in Hcs) and other spin-flip processes, 
the total spin angular momentum is conserved. The spin 
angular-momentum difference between an electron enter- 
ing and leaving a ferromagnet is therefore transferred to 
the magnetization. Under a sufficiently large angular- 
momentum transfer, the magnetization starts to move. 
The component of the net spin angular-momentum flow 
out of the ferromagnet parallel to m vanishes, since 
the outflow cancels the inflow for the parallel component 
(assuming spin along the uniform magnetization direc- 
tion is conserved). The spin-transfer torque t — —Is 
should be accounted for as an additional source term in 
the equation motion of the magnetization. In the pres- 
ence of spin-flip scattering, a component parallel to m 
must be projected out to represent the torque that drives 
the transverse magnetization dynamics: 



T = m X I, X m . 



(12) 



An electron injected into a ferromagnet at the Fermi 
energy and transverse polarization is not in a momentum 
eigenstate, but should be described by a linear combina- 
tion of majority and minority spin eigenstates associated 
with different Fermi wave vectors, kp and kp. The linear 
coefficients of up and down spins carry out oscillations as 
a function of time and position, equivalent to a spin pre- 
cession around the exchange magnetic field. Fermi-level 
states entering the ferromagnet at different angles pre- 
cess on different length scales perpendicular to the inter- 



face, depending on the perpendicular component of the 
spin-up and spin-down wave-vector difference. In ferro- 
magnets with a large cross-section area, a large number 
of transverse modes with different spin-precession lengths 
contribute to the total spin current. The destructive in- 
terference of numerous states with different phases cor- 
responds to the absorption of the transverse spin current 
inside the ferromagnet on the scale of the so-called trans- 
verse spin-coherence length 



Asc — 



k^p 



(13) 



Asc ~ Af (the Fermi wavelength), an atomistic length 
scale for, e.g., transition-metal ferromagnets like Co, Ni, 
and Fe, or their alloys. The smallness of penetration 
depths Asc a posteriori justifies the implicitly assumed 
clean limit. Age ^ A (the mean free path). It should 
be noted that Asn s ets a length scale of a power-law 
ijStiles and Zang will'. '2002^. not an exponential suppres- 
sion of the transverse spin current. 

After transmission through a ferromagnetic film much 
thicker than Age, electrons are completely polarized along 
the magnetization direction. When refiection at the F|N 
boundary may be disregarded, Is on the right-hand side 
of l|12|) is simply the negative of the transverse spin cur- 
rent incident on the ferromagnet. When reflection cannot 
be neglected, the transverse polarization of the reflected 
electrons should be taken into account: although the re- 
flected electrons hardly penetrate the ferromagnet (over 
the Fermi wavelength) , the strong exchange fleld can still 
induce a significant precession of the reflected component 
(^iles and Zanewill, 2002J. This can lead to a reaction 
torque on the ferromagnet as an effective magnetic field 
oriented parallel to the spin accumulation in the normal 
metal. However, at interfaces to transition-metal ferro- 
magnets, positive and negative contributions to the ef- 
fectiv e field typically average out to be small (|Xia et all 
I2002D . 

The dynamics of a monodomain ferromagnet of volume 
V and magnetization that is subject to the torque 
(|12|l arc modified by an additional source term on th e 
right-hand side of the LLC equation ijSlonczews kiL ll99(l) : 



9m 



torque 



-m X I, X m . 



(14) 



For a fixed current density, Eq. (|14(l is proportional to the 
interface area and inversely proportional to the volume 
of the ferromagnet. Current-induced magnetization dy- 
namics are usually realized in perpendicular spin valves 
with one hard (highly coercive) ferromagnet that acts 
as a static polarizer and a second soft ferromagnet that 
responds sensitively to the spin-transfer torque. 
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E. Spin emission by excited ferromagnets 

When seeking a consistent theory of the magnetiza- 
tion dynamics in heterostructures, the current-induced 
magnetization torque discussed above is only one side 
of the coin. A moving magnetization in a ferromagnet 
that is in electric contact with normal conductors emits 
( "pu mps" ) a spin current into it s nonmagnetic environ- 
ment l|Tse rkovnva k et a^J ■ l2002a^ ■ giving a contribution 
to Is in Eq. H14|) . The spin pumping thus leads to an ad- 
ditional source term in the LLG equation, also when the 
magnetization dynamics are induced by external mag- 
netic fields and not by applied current bias. In typical 
biased systems with current-induced dynamics, the spin 
pumping is of the same order as the current-driven torque 
and should be treated on equal footing, as explained in 
SeclEO 

The spin pumping by a precessing ferromagnet is, in 
some sense, the reverse process to the current-induced 
magnetization dynamics. When the pumped spin angu- 
lar momentum is not quickly dissipated to the normal- 
metal atomic lattice, a spin accumulation builds up and 
creates reaction torques due to transverse-spin backflow 
into ferromagnets. The interplay between the magneti- 
zation dynamics and the nonequilibrium spin-polarized 
transport in heterostructures is the central topic of this 
review. The conversion of magnetization movement into 
spin currents and vice versa at a possibly different loca- 
tion is what we mean by the "nonlocality of the magne- 
tization dynamics" in our title. In the remainder of this 
subsection, we put this topic into a historic perspective. 

Nonlocality of the magnetization dynamics can be in- 
terpreted as a nonlocal exchange coupling with explicit 
time depende nce. A first s tep into this direction was 
carried out by iBarnesl Il974f who generalized the RKKY 
theory for the static linear response of the electron gas 
to magnetic impurities to dynamic phenomena in order 
to understand the electron-spin resonance of localized 
magnetic moments embedded in a conducting medium. 
He showed that the dynamic part of the RKKY inter- 
action in diffuse media is limited by the spin-diffusion 
length. A related experimental observation of "giant 
electron-spin-resonance transmission" through a Cu foil 
implanted wit h mag netic Mn ions on one or both sides 
l)Monod et all Il972() showed that precessing impurity 
magnetic moments cause nonequilibrium spin diffusion. 

Subsequently, ISilsbee et all Il979l observed a strong 
enhancement of the microwave transmission through a 
Cu foil with a thin ferromagnetic layer evaporated on 
one side, when the ferromagnetic and Cu conduction- 
electron-spin resonances are tuned into a collective 
mode. This is related to the enhancement of "Lar- 
mor waves" in nonresonant electron-spin transmission 
throu gh normal-metal foils coate d with a ferromagnetic 
layer ijjanossv and Monodl . ll97'6|) . The experiments were 
interpreted by postulating a phenomenological spin in- 
terdiffusion through the F|N interface by nonequilibrium 
components of the magnetization/spin accumulation on 



both sides. These authors concluded that the precess- 
ing magnetic moments can be a source of nonequilib- 
rium spin accumulation diffusing through the nonmag- 
netic conducting medium. Vice versa, the nonequi- 
librium spin accumulation can be transferred into the 
magnetizat ion motion. This picture was investig ated 
further by iJanossvl 119801: iParks an d Silsbeel Il987l and 
was invoke d later to qualit atively i nterpret the exper- 
imen t s bv iHurdeaui nt. 1991; Hurd eauint and Dunifeil 
ll988HHurdeauint and Maloucha 119911 

The discussion of the dynamic coupling between a 
precessing magnetization and itinerant electrons in lay- 
ered F I N stru ctures has been (independently) revived by 
iBergeil Il996i He predicted an enhanced Gilbert damp- 
ing in thin ferromagnetic films in trilayer F|N|F config- 
urations, relying on an elementary quantum process of 
magnon annihilation associated with electron spin flip. 
A very different ap proach to the problem was taken 
by iTserkovnvak et al.. .2QQ2a Th e y used the formalisrn 
of pa rametric pumping ijBrouweil Il998t iBiittiker et all 
\1994} developed in the context of mesoscopic scatter- 
ing problems in order to show that the time-dependent 
magnetization induces spin emission into the itinerant 
degrees of freedom, see Sec. IIIII [A host of other 
mesoscopic spin pumps have been proposed in recent 
years, see, e.g., ISharma and Broiiweii 1200,^ and refer- 
ences therein, at least one of which has been realized ex- 
perimentally l|Watson et all l2003(l .] The spin-pumping 
picture enables us to discuss several topics of this re- 
view in a unified manner and is easily rendered quanti- 
tative. More r ecently, a linea r-response formalism sim- 
ilar to that of 'Barnes, '1974" has been put forward by 
^manck and Hcinrich, 2003. This alternative point of 
view has the advantage to be more familiar to many in 
the magnetism community, but it is much less suited for 
quantitative comparison with experiments, as discussed 
in Sec. ED 



II. SCATTERING-THEORY APPROACH TO 
MAGNETOELECTRONICS 

A. Magnetoelectronic dc circuit theory 

Electron spin and charge transport in F|N heterostruc- 
tures with static magnetic configurations has attracted a 
considerable attention following the discovery of GMR. 
Most of the activities in recent years, including the 
work reviewed here, concentrated on the CPP geom- 
etry in which the electrons pass sequentially through 
magnetic and nonmag netic elements of the circuit, see 
iGiis and Baueil Il997l for a review. A systematic and 
quite general, yet easy-to-handle, semiclassical approach 
to this problem — the magnetoelectroni c circuit the- 
ory — is reviewed by iBrataas et all l2005ai In the follow- 
ing we want to give a brief account of that theory before 
extending it to the dynamic magnetic configurations in 
Sec.lml 
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A basic element of the magnetoelectronic circuit the- 
ory is a magnetic "interconnector" between two nor- 
mal nodes or reservoirs attached to the ferromagnet via 
Ohmic leads, as shown in Fig. |21 Physically, this could 
be realized, e.g., as a layered pillar N|F|N structure. The 
normal nodes are chaotic to the extent that the nonequi- 
librium transport through the leads can be expressed 
in terms of the energy-dependent distribution functions 
/(e) in each node (averaged over orbital states at energy 
e), which are 2x2 energy-dependent matrices in spin 
space of spin-1/2 electrons. (We make a convention of 
denoting such 2x2 matrices in spin space by hats.) In 
thermal equilibrium /(e) = /fd(£)o'o, where /FD(e) is 
the Fermi-Dirac distribution for electrons and ctq is the 
2x2 unit matrix. In linear response, it is convenient to 
define local electrochemical potentials Hr. an d "spin accu- 
mulations" /x^ in all nodes: l|Brataas et a/.L 12000 . 2001) 



N 



deTr[/(e)] , 



rfeTr[^/(e)] 



(15) 



(16) 



choosing a reference energy eo that lies below the Fermi 
energy by much more than the thermal energy and 
voltage biases, but which is arbitrary otherwise, a = 
(ctk, CTy, CTz) is a vector of the Pauli matrices. Notice that 
in our convention, the "spin accumulation" is a vector 
with the direction determined by the total nonequilib- 
rium spin-imbalance density s and the magnitude given 
by the corresponding energy splitting of spins along this 
direction. In linear response and at low temperatures, 
these quantities are related through the Fermi-level den- 
sity of states (per spin and unit volume) in the node, 

M{eF): s={h/2)M{eF)l^s- 

Here we wish to calculate the dc charge and spin 
angular-momentum currents, Ic and Is, entering the 
nodes through the leads, which are induced by the 
nonequilibrium spin accumulations in the nodes and/or 
electrochemical imbalance between the nodes. It is con- 
venient to define the 2x2 tensor current 



(17) 



whose isotropic and traceless components determine re- 
spectively the charge and spin currents. Since, as dis- 
cussed above, spin currents are not conserved at F|N in- 
terfaces, we use the convention that it is calculated on 
the normal side (unless specified otherwise). The 2x2 
current operator /; for the /th lead {I — L, R) can be 
expressed in terms of operators ;(e) [5^ ;(e)] that an- 
nihilate a spin-tT electron with energy e leaving (entering) 
the Zth node through the nth quantum channel of the lead 
ijTserkovnv ak and Brataas, 2001): 



TCTCT 



--y 
h ^ 



dede' 



(18) 
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FIG. 2 Magnetic scatterer (interconnector) connecting two 
chaotic normal-metal (N) nodes via ballistic leads that sup- 
port a quantized number of transverse channels at the Fermi 
energy. The scatterer includes a ferromagnetic (F) region 
characterized by a uniform magnetic direction m (but not 
necessarily uniform magnitude of the exchange spin splitting 
along m), which is depicted as the gray box in the center. 
The dark-gray areas of the ferromagnetic region near both 
F|N interfaces mark the extent of the transverse-spin coher- 
ence characterized by Asc, Eq. Each of the two normal- 
metal regions is divided into a reservoir characterized by the 
electrochemical potential /ic and (vector) spin accumulation 
fj,^, a ballistic lead with a fixed number of transport channels, 
and possibly a disordered region incorporated in the intercon- 
nector (depicted by wavy lines), which accounts for relevant 
normal-metal scattering processes. The interconnector is de- 
scribed by spin-dependent reflection and transmission coeffi- 
cients forming the scattering matrix, Eq. 1211 . for quantum 
channels in the normal leads. The purpose of the formal- 
ism is to calculate the nonequilibrium charge and spin flows 
in the leads, Ic and Is, induced by spin accumulations and 
electrochemical-potential imbalance in the nodes. 



Suppose, as a starting point, that the momentum-space 
distribution in each node is isotropic, i.e., /(e) does not 
depend on orbital quantum numbers. This is true if 
the nonequilibrium currents do not cause significant drift 
contributions to the distribution function. (We will re- 
visit and drop this assumption in Sec. lILBl l For the ^th 
lead then 

[aUe)^a:;,,,{s')) = ff" {e)Su'S^^,5{e - e') , (19) 

and it is now straightforward to evaluate the expectation 
value {I^"' ) (also denoted simply by Ij^'^ ) of the current 
operator, after relating the scattered states to the in- 
coming states via the scattering matrix of the magnetic 
interconnector: 



bU^) 



E 

a' n'l' 



^nn' ,11 



(20) 



The scattering coefficients s'^, n, characterize reflection 
ii I = I' and transmission otherwise. Eq. (|20|l assumes 
that the entire interconnector is elastic, so that the elec- 
tron energy is conserved upon scattering between the nor- 
mal nodes. For a ferromagnet with an exchange spin 
splitting along unit vector m and vanishing spin-orbit 
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interaction in the system ijBrataas et a/.l . l2000(l . 



(21) 



in terms of the scattering coefficients for spins up (down) 
along m, s„„/ ;;/, and the projection matrices 



-,T(i) _ 



= - (o-o ± <T • m) 



(22) 



A consequence of the elastic-scattering approximation is 
a rigid exchange field: there are no magnons excited by 
the electron transport. If we furthermore assume a suffi- 
ciently low temperature, voltage imbalance, and spin ac- 
cumulations, so that the scattering-matrix variation on 
these energy scales is negligible, we can replace Snn',u'{£) 
by its value at the Fermi level ep- 

It is then convenient to group the conductance parame- 
ters into two pairs of 2 x 2 matrices. For electrons incident 
from the right lead, we define 



nn \ nn J 



(23) 
(24) 



where the index n' is summed over the channels in the 
right lead, and n runs over the channels in the right lead 
in Eq. and left lead in Eq. The coeflicients 

r^„/ and ij^^, are reflection and transmission amplitudes, 
i.e., elements sj^^, of the scattering matrix |[?T)) with 
I 7^ I' and I = r, respectively. The range of summa- 
tion for n' in Eqs. H23|l and (|24|l . i.e., the total number 
of transverse quantum channels in the right lead (for a 
given spin species at the Fermi level) is called the Sharvin 
conductance, g^^, a quantity which will be useful later. 
For electrons incident from the left lead, we denote the 
reflection and transmission amplitudes by primed (quan- 
tities, and we similarly define primed matrices g'"'"' and 
t'"'"' in terms of the primed scattering amplitudes. We 
denote the Sharvin conductance of the left lead by g'^^. 
It should be understood that all scattering coefficients 
and corresponding conductance parameters are evaluated 
at the Fermi level. Putting ab ove equations together 
llBrataas et 1200(1 12OO11 120031) . 



{2(gTT + giX) (^,3-^,,l) 



r(0) , 

_L|2(gTT „gU)(^^j^_^^ J, 



(25) 



t(0) 

+ 2gl^m X /X, R X m + 2.gJV^,R x m 
-2<Tim X /X, L X m - 2<f x m} , (26) 

and the currents through the left lead are obtained by 
interchanging L ^ R, g°'"' ^ g""' and t'^'^ ^ t'^K 



[The superscript (0) introduced here denotes currents 
with static magnetizations.] By unitarity of the scatter- 



mg matrix, g 



9" 



t''^'^ , if the spin compo- 



nent along the magnetization direction is conserved. The 
dc transport in the two-terminal geometry is then deter- 
mined by two real-valued spin-dependent conductances 
5'^'^ and four complex-valued (spin-)mixing parameters 
gll — gU _|_ igji-^ ftl — ^Ti _|_ j^Ji subscrlpts r and i 
respectively denoting the real and imaginary parts), g'^^ 
and t''^K For a mirror-symmetric structure, g^^ ~ g'^^ 
and i^J- = t'^^. 

We have now also access to microscopic expressions for 
the spin-transfer torques H12|l . The torque on the right 

surface of the ferromagnet tr = — m x 1^°^ x m: 



t;Tim X . X m - V, t x m ) , (27) 



"^R =^ I al^^ X /X, R X m-f 5I Vs.R X m 



is proportional to the spin-mixing (i.e., off-diagonal) com- 
ponents of the conductance matrices H23() and (|24|1 . The 
first two terms in Tr involve reffection at the right F|N 
junction and the last two terms — transmission through 
the entire N|F|N structure. The latter terms can thus 
be disregarded when the ferromagnet is much thicker 
than the transverse-spin coherence length Asc, Eq. H13|l . 
since transmitted electrons accumulate phases differing 
by more than tt for opposite spins (along m). In that 
limit, the first term, proportion al to m x p x m, is sim- 
ilar to the torque introduced bv lSlonczewsk that is 
responsible for instability leading to magnetization pre- 
cession or reversal. The second term, proportional to 
/Xg p; X m, acts as an effective magnetic field collinear 
with the spin accumulation in the right normal node. 
In transition-metal ferromagnets, gl'^ < Q.lgl^, see, e.g.. 
Table n and Sec. IVII.A.Tl so that the effective magnetic 
fleld can be disregarded in many practical situations. 

For ferromagnetic films much thicker than Asc, the re- 
maining mixing conductances g^^ and g'^^ are insensitive 
to scattering processes deep inside the ferromagnet (i.e., 
in the light-gray area of the ferromagnet in Fig.|2J) and are 
determined by the scattering potential of a thin slice of 
the ferromagnet near the interfaces (the dark-gray areas) 
and eventually the normal metal (the wavy-line area). 
The mixing conductance is then a property of the isolated 
F|N junction rather than the entire N|F|N scatterer. In 
this limit, we may introduce a ferromagnetic node at a 
sufficient distance from the interfaces and consider with- 
out lack of generality the two interface separately. This 
is allowed, since in the bulk of the ferromagnet the spin 
accumulation becomes a well-defined semiclassical dis- 
tribution function collinear with the magnetization. An 
analysis entirely analogous to the one leading to Eqs. H25|) 
and (|26|l results in the same equations for the currents 
Ic and Is (on the normal-metal side of the right inter- 
face). Only the terms proportional to t'^^ and t'^^ in 
Eq. (|26|l drop out now because of the coUinearity of /Xs,l 
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and m, left node no\y being : assum ed ferromagnetic, see 
Fig. 151 ijBrataas et all I2OOOI l200l|) . Naturally, when the 
F-layer thickness is much larger than Age , the mixing con- 
ductance for the normal lead g^'^ in the F|N structure of 
Fig. O is the same as g^'^ for the right lead in Fig. [21 
(assuming same F|N contacts including dark-gray and 
wavy-line areas, at the right lead). On the other hand, 
the spin-up and spin-down conductances of the N|F|N 
structure are not identical to the conductance parame- 
ters of a single F|N interface, being dependent on two 
F|N junctions and the bulk F layer. 




FIG. 3 Contact between a ferromagnetic and a normal node. 
The notation is analogous to Fig.|2| Here, the charge and spin 
currents in the normal lead depend on the conductance matrix 
g'^'^ defined in terms of the spin-dependent reflection coeffi- 
cients from the normal- metal side, as before, see Eq. H2.'i|l . 
The nonequilibrium spin accumulation in the ferromagnetic 
node is coUinear with the magnetization m. 

So far, we focused on the spin and charge flow through 
a single resistive element, i.e., a single N|F|N or F|N junc- 
tion. Equations H25|l and (|26|l . which can be viewed 
as a gene raliza t ion o f the Landauer-Biittiker formulas 
(see, e.g.. Ilmrvl Il997j) to describe charge and spin cur- 
rents in a two-terminal geometry with possibly nonequi- 
librium spin accumulations in the nodes. These are basic 
building blocks of the magnetoelectronic circuit theory. 
We also need the Kirchhoff laws generalized to consider 
the spin currents and spin accumulations on equal foot- 
ing with the usual charge currents and voltage biases 
|Brataas et al . 2005a). The properties of a given device 
or circuit can be calculated by first prudently separating 
it into reservoirs, nodes, and resistors (interconnectors), 
where the latter are the current-limiting elements. The 
nodes are supposed to have a negligibly small resistance 
and their choice may depend on the problem at hand. In 
a disordered multilayer, for example, it is convenient to 
imagine (fictitious) nodes at both sides of an interface, 
treating the latter as a separate resistive element (whose 
conductance parameters may need to be redefined, how- 
ever, as explained in Sec. III.B]l . Reservoirs represent the 
"battery poles," i.e., large thermodynamic baths at ther- 
mal equilibrium with a constant electrochemical poten- 
tial. The electrochemical potentials and spin accumula- 
tions in the nodes are not known a priori, but determined 
from the generalized Kirchhoff laws based on the conser- 
vation of spin and charge. For example, disregarding 



spin-flip scattering, in the stationary dc state, all spin 
and charge flows into a normal-metal node must vanish: 

^/,,,=0, ^ls,;=0, (28) 
I I 

summing over all leads attached to the node. The net 
spin flow into a ferromagnetic node projected onto its 
magnetization also vanishes, whereas the transverse cur- 
rents are absorbed at the interface, as discussed above. 
We have also seen that the spin and charge currents in 
each lead can be calculated as a function of the distri- 
butions on adjacent nodes/reservoirs in terms of well- 
deflned conductance parameters. The spin- and charge- 
current conservation laws H28|l then allow computation of 
the circuit properties as a function of, e.g., the voltages 
applied to the reservoirs. The protocol for calculating 
the current-voltage curves may be summarized as: 

1. Divide the circuit into resistors, nodes, and reser- 
voirs. 

2. Specify the control variables, e.g., the voltages 
(electrochemical potentials) applied to the reser- 
voirs. Parametrize the electrochemical potential 
and spin accumulation of each node. 

3. Compute the currents through the resistors as a 
function of the distributions in the adjacent nodes, 
which requires the spin-dependent and spin-mixing 
conductances defined earlier. 

4. Use the spin- and charge-current conservation laws 
H28|l at each node. If there is spin decoherence, use 
a modified continuity equation with spin relaxation. 

5. Solve the resulting system of linear equations to ob- 
tain all currents as a function of the electrochemical 
potentials of the reservoirs. 

The charge and spin currents through the resistors, the 
net spin torques on the ferromagnets, and the spin accu- 
mulations anywhere in the circuit can be computed this 
way. 

As a specific application, consider the layered spin- 
valve structure and its effective circuit models sketched 
in Fig.Q] The reservoirs are described by electrochemical 
potentials /ic,L and fJ-c.K- We can define four nodes on the 
normal side of each F|N interface, with electrochemical 
potentials /ic,i and spin accumulations /x^j. There are 
correspondingly three normal resistors in the problem, 
each described by a single real-valued conductance pa- 
rameter g, and two magnetic resistors, corresponding to 
two ferromagnets, each described by two real- valued spin- 
dependent conductances and g^^ and four complex- 
valued spin-mixing conductances g'^^, g'^^, and t^^, t''^K 
In the absence of spin-orbit coupling, the charge and 
spin currents in each node depend on the electrochem- 
ical potentials and spin accumulations on each side of 
the respective ferromagnet according to Eqs. H25|l and 
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(|26|l . Transport across normal resistors is described by 
simpler equations that could be easily obtained by set- 



ting 



9" 



t" 



g in Eqs. ^ and (EEJ. 



By following the steps outlined above, these equations 
can be used to self-consistently determine fic,i and /x^ ^ 
and then charge and spin currents into each node, as a 
function of the applied bias /ic,L — /^c,r- Obviously, the 
same procedure can be carried out for current-biased in- 
stead of voltage-biased systems. When the conductance 
parameters are to be evaluated from first principles, the 
definitions 1)23(1 and ((24(1 have to be corrected for kinetic 
effects when the interfaces are highly transmitting, as 
discussed in the next subsection. Placing of the nodes is 
to some extent arbitrary but different choices should not 
lead to contradictory results. For the system in Fig.Q] for 
example, it might be more convenient to replace nodes 2 
and 3 with a single one somewhere in the middle of the 
normal spacer and eliminate nodes 1 and 4 altogether, 
while redefining the conductance matrices of the mag- 
netic regions to include scattering on the normal-metal 
sides. In such a case, the entire structure would consist 
of two magnetic scatterers connecting each reservoir with 
the middle of the normal spacer. On the other hand, for 
sufficiently thick ferromagnets, we may introduce four ad- 
ditional nodes on the ferromagnetic side of the interfaces, 
and define pure F|N- interface resistances and ferromag- 
netic bulk resistances. 



Circuit theory assumes momentum scattering in each 
node but not necessarily inelastic scattering. It is only 
required that the chemical-potential gradients/drops are 
small enough such that linear- response theory holds, and 
the energy dependence of the scattering matrix may be 
disregarded. Inelastic scattering in the nodes per def- 
inition does not affect the transport properties. It is, 
however, often stated that the scattering approach to 
transport is valid only when the orbital-dephasing length 
is sufficiently longer than the dimensions of the scat- 
tering region. This is a relevant statement only when 
phase coherence is essential for the physics under con- 
sideration. In magnetoelectronic systems, we have to 
worry about orbital interference only at sharp interfaces. 
The spin-mixing conductances g^^, g'^^ and t^^, t'^^ in 
Eq. 1(26(1 govern transverse spin currents under the con- 
dition that » Asc. Also, the description of con- 
ventional spin-dependent conductances by first-principles 
band-structure calculations assumes wave-function co- 
herence on atomistic length scales as well. These con- 
ditions are obviously hardly restrictive and assumed to 
be valid up to high (room) temperature. When spin-flip 
scattering is strong in diffusive regions, the spin-diffusion 
equation should be solved their, for which Eqs. 1(25(1 and 
provide the boundary conditions. Sec. lIV.Bl demon- 
strates how this is carried out in practice. Although we 
illustrated here circuit theory using simple linear struc- 
tures, it can be easily applied to more general multiter- 
minal devices, such as spin transistors. 



B. Interfacial and thin-film conductance matrices 

In the previous subsection, the electron states in the 
nodes were assumed to be occupied according to energy 
and spin, but without any regard to their momentum. 
Physically, this "isotropy" in momentum space implies 
that net currents in the nodes may be disregarded. This 
is allowed only when the incoming and outgoing cur- 
rents do not significantly disturb the isotropic momen- 
tum space distribution. When the contacts to the nodes 
are relatively small (point contacts) or highly resistive 
(tunnel junctions), this approximation holds. Highly- 
conductive metallic multilayers that are the main sub- 
ject of this review, do not satisfy such a condition. The 
node interconnectors are then intermetallic interfaces or 
thin films. The isotropy of the distribution functions in 
the nodes can be significantly distorted by the current 
induced by a given voltage bias. In that limit, simple 
Kirchhoff laws with Landauer-Biittiker conductances pa- 
rameters, Eqs. (|^ and (|^ . do not apply. This sub- 
section summarizes how to rescue ci rcuit-theory by only 
modifying conduct ance parameters ijBauer et all l2003bt 

IS£hie_iLfiD,llM3)- 

We first illustrate the issue for a nonmagnetic metallic 
pillar with a uniform cross section connecting two reser- 
voirs. Transport through a ballistic pillar is governed by 
the Sharvin conductance g^^, i.e., the number of prop- 
agating transport channels. Let us introduce M inter- 
faces (e.g., grain boundaries) in series that scramble the 
transverse-momentum distribution of incident electrons 
without any significant backscattering. The total con- 
ductance must then still am ount to g^^, since w excluded 
any reflection l)lmrvl . [r997f) . If we naively carried out the 
circuit-theory protocol, we could assign a conductance 
9i = 9^^ to each interface and insert M ~ 1 nodes be- 
tween them. We would trivially obtain that the total 
resistance l/g is the sum of the individual ones: 



M 



9 



M 

"Sh 



(29) 



which is obviously wrong. This breakdown of the circuit 
theory can be "flxed" by renormalizing the individual 
interface resistance from scattering theory 1 / gt by sub- 
tracting the Sharvin resistance 1/g^^: 



1 
91 



1 

9i 



9' 



1 

sir 



(30) 



At the same time, the contact resistances between pillar 
and reservoirs are assigned half of their Sharvin resis- 
tance. We then arrive at the total resistance 



1 



M 



1 



1 



Sh 



25 



Sh 



9 



1 

sir 



(31) 



In this way, we enforced a vanishing local voltage drop 
over each interface, corresponding to jji = oo, and a fi- 
nite Sharvin resistance for the entire structure that is 
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governed by the geometrical sample cross section. The 
total conductance is g^^, as it should. 

ISchep et all Il997l justified such "renormalized" resis- 
tances, deriving a more general result for interfaces be- 
tween different materials: 



renormalized as 
1 



1 

9i 



Sh 



Sh 



(32) 



where gf^ and gf^ are Sharvin conductances on two 
sides of the interface, thus allowing to compare first- 
principles calculations of interface resistances with exper- 
imental results on diffuse multilayer systems that access 
gi. Eq. H32|l requires global diffusivity that separately 
randomizes the momentum distributions for right- and 
left-moving electrons and destroys quantum interference 
between consecutive contacts. 

The renormalization of resistances by subtracting 
Sharvin contributions is closely related to early contro- 
versies about resistance in mesoscopic systems. The 
"bare" Landauer-Biittiker conductance g = X^nn' l^nn'P 
is suitable for the description of two-point measurement 
in which reservoirs are current as well as voltage sources 
(see, e.g.. Imrv. 1997). In an idealized four-point mea- 
surement where the voltage drop is measured directly 
across the scatterer (e.g., an interface), the conductance 
should be renormalized since there is no geometric con- 
tribution. In a single- mode quantum wire with g^^ — 1, 
the renormalized conductance reduces to the "o ld" Lan- 
dauer formula g = g/{l — g) ijLandaueil Il97(]() . When 
the scatterer is embedded in a diffuse environment and 
circuit theory applies, the situation is analogous to the 
four-point measurement and Eq. (|32|l should be used. In 
coUinear magnetic structures, Eq. (I32II should be applied 
for each spin channel separately (Schep et al, 1993)- 

The "kinetic" corrections for magnetic structures with 
noncoUinear magnetizations require additional thought. 
As a starting point we consider a single F|N interface as 
sketched in Fig. |2| with conductance matrices 3°'°' and 
disregard spin-flip processes. In the presence of nonequi- 
librium currents, the effective electrochemical potentials 
and spin accumulations, Eqs. (fTK|l and lfTB|) . on either side 
of the interface are different for left- and right-moving 
modes. The discussion of Sec. III.AI must thus be modi- 
fied to account for this asymmetry in the left- and right- 
moving distributions ("drift"). Like in lSchep et all\l99^ 
the central underlying assump tion that allows an easy im- 
plementation of this program (Baue r et all |2003aD) is a 
globally diffuse or chaotic system. In that case, the distri- 
bution of electrons incident on each interface from either 
side is isotropic in momentum space. Such a random- 
ization is likely to be provided in realistic structures by 
interfacial disorder and bulk scatter ing. By extens i on of 
either the random matrix th eory of IWain tal et all ^2 0001 
or circuit theory, the analysis ijBauer et a t. 2003a b) then 
boils down to a set of simple rules. Eqs. H25|) and H2t)|) 
remain unchanged but the interfacial conductances are 



-o -sl7 + -^h (33) 



9f 



Sh 



whereas the average electrochemical potentials and spin 
accumulations in the nodes have to be found self- 
consistently as before. (We recall that the conductance 
parameter t''^^ drops out of the discussion for a single F|N 
interface.) Additionally, the leads to real reservoirs must 
be assigned half of the respective Sharvin resistances, for 
each spin a. 

The renormalization H33I) is highly significant for un- 
derstanding of the measured resista nces of metallic mul- 
tilaye rs in the CPP configuration ijBass and Pratt. Jrl 
I1999D . as well as other m easurement s, by ab ini- 
tio calcul a tions ite^^ere^^ . 2003b: Schep et all Il997t 
IXia et a.l\ . bool 120011 IZwierzvcki et all l200,'^^. In Ta- 
bleUwe quote the theoretical results (jZwierzvcki et'cdl . 
2005) for two representative N|F material combinations: 
Au|Fe(001) and Cu|Co(lll), the former routinely used 
by the Simon- Eraser grou p (e.g., iHeinrich et al\ . l2002l 
l2003bt lUrban et all I2OOII) and the latter by the C ornell 
group fe.g.. lKatine et a/.l. l2000HMvers et a/.l. ll999'). The 
large difference between the g's and g's is evident, when 
using Eq. (|33|l . as the bare interfacial conductances are 
comparable with the Sharvin conductances. 
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Ti 


Ti 

gi 


5n 




gpi 


Aufcc Febcc 


clean 


10.3 


2.1 


12.0 


0.1 


11.9 


21.4 


11.9 


(001) 


alloy 


10.1 


4.6 


11.9 


0.1 








CUfcc|COfcc 


clean 


10.8 


9.8 


14.1 


0.4 


15.0 


11.9 


27.9 


(111) 


alloy 


10.8 


8.5 


14.6 


-1.1 









TABLE I Calculated interface conductances (in units of 
quantum channels per nm^). The results are shown for clean 
and disordered interfaces. The latter are modeled by two 
atomic monolayers of 50% alloy. FromfZwicrzvcki et al. , 2005|. 

Eq. (|33|l can be used only when the distribution func- 
tion in the ferromagnetic layers is well defined. This is 
not the case anymore when the magnetic-film thickness 
is of the order or smaller than the transverse-spin co- 
herence length II13II. I n the latter situation, calculations 
llBauer et al I l2003al b') that lead to Eq. ^ have to be 
repeated for two normal-metal layers separated by a thin 
magnetic film as in Fig. [3 We find that the basic circuit- 
theory equations (|25|l and (|26() hold for all film thick- 
nesses after replacing the conductance parameters for the 
N|F|N sandwiches with the renormalized ones, g°'°' and 
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(34) 
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and the same after interchanging g ^ g' and t ^ t'. 

In the presence of weak (compared to momentum 
scattering) spin-flip scattering during bulk diffusion, the 
same renormahzations, Eqs. (|34|l . still hold for in- 

terfacial transport. In the opposite limit of high spin-flip 
rates, the layers can act as ideal sinks for spin currents. 
In particular, in the regime of pure spin transport, such 
layers are fully equivalent to reservoirs, bearing all con- 
sequences and rules outlined in this subsection for the 
reservoirs. The regime of intermediate-strength spin de- 
phasing cannot be characterized by such simple state- 
ments and is omitted from our discussion. 



C. Time-dependent theory 

In this review we are mostly interested in the dynamic 
phenomena stemming from a slow variation in the mag- 
netization direction m of ferromagnets that are part of 
an Ohmic circuitry. We focus on the adiabatic response 
of the itinerant carriers to the time-dependent m that 
is driven by external magnetic flelds or applied (spin) 
currents. It turns out that the magnetoelectronic (dc) 
circuit theory discussed above contains already all the 
necessary parameters in the adiabatic regime. An adi- 
abatic approximation is applicable when the frequency 
of the magnetization modulation is much smaller than 
the characteristic ferromagnetic exchange spin splitting, 
which is safely fulfilled for transition-metal-based struc- 
tures. The total current is a sum of the currents induced 
by the bias applied via the reservoirs (including the spin- 
transfer torques), viz., Eqs. iPsjl . and the "pump- 
ing" component proporti onal to the rate of va riation of 
the scattering potentials teiittiker et al\ . \l994} : 



(35) 
(36) 



We will make a convention to include in /i"' and li"' 
the currents driven by any charge and spin imbalance 
brought about by the pumping, as well as by the applied 
bias. 

In Sec. mil we derive the currents pumped by a time- 
dependent magnetization direction m into adjacent nor- 
mal nodes: 



(37) 



introducing a new complex- valued parameter 

A^^ = Al^ + iAj^ = g^^ - t'Ti (39) 

that determines the magnitude of the spin pumping 
as a fu nction of device parameters l)Tserkovnvak et al\ . 
l2002a|) . Magnetism (or another source of spin-dependent 
scattering) is essential for a nonvanishing pumping pa- 
rameter A^^ ■ The spin current into the left lead is similar 



to Eq. (|2H1) and governed by A''^^ = g''^^ - f^. When the 
ferromagnet is much thicker than the coherence length 
(|13|l . the quantities t'^^ and t'^-^ vanish and the mixing 
conductances g^^ and g'^^ do not depend on the thick- 
ness of the ferromagnet; the spin pumping currents orig- 
inate from the interfaces. The spin pumping H38|) thus 
does not depend on spin-flip processes in the ferromagnet 
when far from the F|N interface on the scale of the co- 
herence length H13(l or, in other words, when the spin-flip 
scattering rate is small compared to the exchange split- 
ting. This is typically the case in real materials. The 
spin (de)coherence in the attached normal conductors is 
crucial, however, since it affects the self-consistent reac- 
tion torque exerted on a slowly-precessing monodomain 
ferromagnet, as explained in Sec. II VI 

It can be shown that the rcnormalizations (|34|l in- 
troduced in the dc circuit theory for layered structures 
must be also applied to the mixing conductances in A^^, 
Eq. H39|l. for multilayers with diffuse spin backscatter- 
ing, but not for contacts to spin-sink reservoirs. Using 
Eqs. (|S3|), we find for A^^ = - i'^^ 



iTi 5Ti(2g^Sh„^^,U)„iU(25^Sh_,,Ti) 



23: 



Sh 



{2gl,^ ~ gn)i2g'^i^ - gni) ^ ttlt'U ' 



(40) 



and the same after interchanging primed and unprimed 
quantities. For mirror-symmetric structures, Eq. (|40|l re- 
duces to simply 



1 



1 

A^ 



(41) 



The spin-pumping expression l|38|) sets the stage for 
several interesting developments reviewed in the follow- 
ing. We therefore devote the entire section IIIII to de- 
rive Eq. and discuss its phys ical content, see also 
iTserkovnv^ ^ni\.\2m2fiH.\2()()U 



III. SPIN EMISSION BY COHERENTLY-PRECESSING 
FERROMAGNETS 

A. Parametric spin pumping 

The adiabatic spin-pumping expression H38|) can 
be derived by t ime-de pendent scattering theory 
l|Tserkovnvak et all l2002aj) . Consider an N|F|N junc- 
tion as in Fig. 12 where the two normal nodes are 
now assumed to be large reservoirs in a common 
thermal equilibrium. Without a voltage bias, no spin 
or charge currents flow when the magnetization of the 
ferromagnet is static. When it is moving, however, 
the time dependence of the scattering matrix in spin 
space can induce nonequilibrium spin currents in the 
nonmagnetic leads. The current I(t) pumped by the 
precession of the magnetization into the right and 
left paramagnetic reservoirs can be calculated in the 
adiabatic approximation, since typical precession fre- 
quencies are several orders of magnitude smaller than 
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the ferromagnetic exchange field that sets the relevant 
energy scale for spin-dependent transport. The adiabatic 
charge-current response in nonmag netic systems has 
been derived by IBiittiker et al\ . \l994 . The generalization 
to the 2x2 spin- and ch arge-current matrix Ijl?!) by 
iTserkovnvak et aD . I20fl2albl is explained in the following. 

The 2x2 current operator /; for the Ith lead {I = L, R) 
is in general given by Eq. (|18(l . When the scattering 
matrix sJ^J^, (t) of the ferromagnetic layer varies slowly 
compared to the relevant microscopic time scales of the 
system, an adiabatic approximation may be used, mean- 
ing that the energy of the scattered states is assumed 
to be weakly modulated with respect to the energy of 
the incoming states by the oscillating part of the scat- 
tering matrix: The state annihilated by a^, (e, t) — 
'^n' (e)e~'^*/'' is partitioned into states in the mth chan- 
nel of lead I' with energies determin ed by the time depen- 
dence of s^ri'/,i;'(£' (^'O; see IBiittiker et all Il99l 
The scattering amplitude at a given energy shift is de- 
termined by the Fourier transform of sJ^J^/ (e, t) in time 

space. The expectation value of the current operator /; is 
evaluated similarly to the dc limit discussed in Sec. Ill.AI 
When the scattering matrix depends on the real-valued 
parameter X{t), the Fourier transform of the current ex- 
pectation value Ii{u!) — J dte^'^* Ii{t) can be written as 



Ii{uj) = gx,i{uj)X{uj) 



(42) 



in terms of the frequency lo and X-dependent parameter 
9x,i- 



5x,i(w) 



de 



dhviie) ds 
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nn'l' 
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dX 



H.c. 



nn',ii'(^) 

(43) 



Eq. (|42|l is the first-order (in frequency) correction to 
the dc theory of Sec. Ill.AI At sufficiently low tempera- 
tures, one can approximate — 9e/FD(£) by a J function 
centered at the Fermi energy. The expectation value 
of the 2x2 particle- number operator Qi{uj) [defined by 
Ii{t) — dQi{t)/dt in time or by Ii{lo) = ^iujQi{Lu) in 
frequency domain] for the Ith reservoir is then given by 



\ nn'l' 



H.c. X{lj) , (44) 



where the scattering matrix is evaluated at the Fermi 
energy. Because the prefactor on the right-hand side of 
Eq. 144() does not depend on frequency the equation is 
also valid in time domain. The change in particle number 
SQi{t) is thus proportional to the modulation SX{t) of 
parameter X and the 2x2 matrix current reads 



dhi dX{t) 
'dX dt 



(45) 



introducing the matrix "emissivity" into lead I: 



dm 



1 
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nn' ,11' 



H.c. 



(46) 



The right-hand side of Eq. H45|) should be in general 
summed over all parameters, if any, that modulate the 
scattering matrix. When the spin-flip scattering is disre- 
garded, the scattering matrix s can be written in terms 
of the spin-up and spin-down scattering coefficients sT(i) 
using the projection matrices, see Eqs. 1)21(1 and H22|l . 
The spin current pumped by a magnetization precession 
is obtained by identifying X{t) = f{t), where ip is the az- 
imuthal angle of the magnetization direction in the plane 
perpendicular to the instantaneous rotation axis. For ex- 
ample, suppose the magnetization rotates around the z 
axis: m = (cos(p, sin(^, 0). Using Eq. H21|l. it is then easy 
to calculate the emissivity (146(1 for this process: 
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AI^{<Jx sinip — ay cos ip) — A\^(Tz 



(47) 



Inserting this into Eq. ((45() and comparing the result with 
the definition ((171) . one finally finds that the charge cur- 
0, and the spin current 



rent vanishes, Z^™'' 



Trpump 



[-AY sin ip, AY cos LP, Al^^ 



h dip 
47:~dt 



(48) 



can be rewritten as Eq. ((38(1 . Since the spin current trans- 
forms as a vector in spin space, Eq. ((38(1 is also valid in 
the case of a general (slow) motion of the magnetization 
direction. 

Even though the mathematics of the scattering ap- 
proach to adiabatic spin pumping is entir ely analogous to 
the charge-pumping theory developed bv lBrouweilll998t 
IBiittiker et a^J . Il994l there is an important difference in 
the physics. In the case of a nonmagnetic scattering re- 
gion, the average pumped charge current has the same 
direction in the two leads, because charge cannot be ac- 
cumulated or depleted for long; the net charge entering 
the scattering region through one lead must leave it via 
the other within a period of the external-gate modula- 
tion. On the other hand, the total conduction-electron 
spin angular momentum does not have to be conserved 
since the magnetization can act as a sink or a source. The 
preceding analysis shows that a precessing ferromagnet 
polarizes adjacent nonmagnetic conductors. The phe- 
nomenon can be called as well a "spin well" or "spin 
fountain" rather than spin pump: An excited ferromag- 
net emits spins into all adjacent conductors, spending its 
"own" angular momentum rather than pumping spins 
from one lead into the other. If the lost angular momen- 
tum is not replenished by an external magnetic field, the 
precession invariably will die out. 

Per revolution, the precession pumps an angular mo- 
mentum into an adjacent normal-metal layer which is 
proportional to Al^, in the direction of the precession 
axis. At first sight, it may be surprising that a net dc 
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pumping can be achieved by a si ngle parameter varying 
in time, whereas it can be shown (|Brouweilll998|) that at 
least two periodic parameters are required for that. How- 
ever, there are actually two periodic parameters (out of 
phase by 7r/2) hidden behind ip{t), viz., the projections 
of the unit vector defined by ip in the plane perpendicular 
to the axis of precession. 



B. Rotating-frame analysis 

It turns out that Eqs. H38|l . H39|l for the noncquilibrium 
spin currents pumped out of a precessing ferromagnet can 
be understood as an equilibrium property by a transfor- 
mation to a rotating frame. We prove this for a ferromag- 
netic film sandwiched between two normal-metal nodes 
at a common thermal equilibrium, /ic,L = A*c,R) as shown 
in Fig. 13 Let 



H{t) = Ho + H'{t) 



(49) 



be the Hamiltonian experienced by conduction electrons, 
where Hq is the sum of the kinetic and potential energies 
and H'(t) is a time-dependent exchange Hamiltonian in 
the ferromagnet: 



m(r,i) ^-^'(0 



(50) 

where Vx is the local exchange interaction along the mag- 
netization direction m(r, t) = m(i) that is assumed to 
be uniform in the ferromagnet. ^'o- is the spin-tr electron 
field operator. 

At time t < 0, the entire system is in a common- 
thermal-equilibrium state corresponding to a time- 
independent magnetization direction m(0). At t — 0, 
m{t) starts to rotate with frequency uj about an axis de- 
noted z that is perpendicular to m, m = x m. The 
time-dependent magnetization drives the electron system 
out of equilibrium. 

Let us introduce a (many-body) state ^'{t) of the elec- 
trons in the rotating frame of reference for the spin vari- 
ables, which is related to the solution of the Schrodinger 
equation ih<^{t) = H{t)^{t) by the unitary transforma- 
tion $(t) = U{t)<i>'{t). Here, U = cxp{-iS^ujt/h) is a 
spin-rotation operator in terms of the total-spin z-axis 
projection S^. Since H'{t) = U(t)H' (0)UHt ), <^'{t) sat- 
isfies a Schrodinger equation l)Slichteii Il99(]() : 



ih 



d$'(t) 
dt 



[F(0) - SM $'(*) 



(51) 



provided [Hq,Sz] — 0, i.e., spin-orbit coupling or other 
spin- flip processes are disregarded. This is the equation 
of motion for electrons subjected to a static magnetic 
field Lo/^ in the z direction, in addition to a static ex- 
change interaction (|50|l along m(0) in the ferromagnet. 
A steady-state solution for the system in the rotating 



frame is thus characterized by spin-polarized normal- 
metal nodes with (sz) = and (sx) = (sy) = 0, 
where x is their static (isotropic) spin susceptibility. 

In the laboratory frame, the spin imbalance must be 
identical to the polarization in the rotating frame along 
the z axis. Since there is no magnetic field in the labo- 
ratory frame, this spin imbalance can only be accounted 
for by a noncquilibrium spin accumulation, i.e., differ- 
ence between the chemical potentials for spins parallel 
and antiparallel to z, of magnitude 



(52) 



According to the magnctoelectronic dc circuit theory, the 
spin accumulation in turn drives a spin current 



rback 



1 

47r 



+A^V. X m 



(53) 



from a normal- metal node back into the ferromagnet, 
where the coefficients A^^ have been defined in Eq. I|39|l . 
Here, we used Eq. ()26|l with fXc,L — Mc,r = and /x^ l — 
A*s.R — Ms- According to Eq. /x^ = huz = hm x rh. 
Since in the steady state the backflow spin current H53|) 
must cancel the spin current pumped by the ferromag- 
net into the normal metal, above arguments thus provide 
an alternative derivation of the spin pumping l|38|) . For 
finite-size nodes and slow rotation, equilibrium may be 
established in the rotating frame in a time short com- 
pared to the period of rotation. In the adiabatic limit, 
the equality of the spin currents into and out of the fer- 
romagnet thus holds instantaneously, and not just over 
the average over a period. The assumption of a fixed ro- 
tation axis is also not essential, as long as this axis moves 
slowly in time. In the limit of vanishing exchange split- 
ting, A^'^ ~ and the spin current (|38|) vanishes, as it 
should. 

It is instructive to scrutinize the energy and spin- 
angular-momentum conservation laws of the combined 
F|N system in the lab oratory frame of reference 
ijTserkovnvak et aTLl2002bj ). Suppose the rotation axis z 
is oriented along the static effective magnetic field Hes, 
and let us focus on the noncquilibrium dynamics in the 
normal metal node, see Fig. ^ After starting the mag- 
netization rotation a,t t — 0, net spin currents flow into 
the normal mode in a finite time interval before reach- 
ing the steady state. We wish to account for the energy 
and angular-momentum transfer between the ferromag- 
net and normal metal in this transient regime. If the nor- 
mal node is sufficiently small, the steady state is reached 
after a vanishingly small transfer of spins, with only lit- 
tle effect on the ferromagnet. Ns spins oriented along 
the z axis transferred into the normal metal correspond 
to an excess energy AE^ — and angular mo- 

mentum ALat = Nsh/2. The conservation laws dictate 
AEp = — Ai?jv and ALp = — Aijv for the ferromagnet. 
Using Eq. lO, the magnetic energy AEp — jALpHcS 
and we find NsHs/^. = ^Ns{h/2)Hcs. It follows that 
/is — hjHes = So;, where co = 7-ffcff is the Larmor pre- 
cession frequency in the effective field. Using energy and 
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momentum conservation, we thus derived a nonequilib- 
rium spin accumulation that agrees with the result of the 
rotating- frame analysis. 
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FIG. 4 The ferromagnetic magnetization m rotates around 
tiie effective field Heft. In the steady state, the spin pumping 
jpump (,a,ncels the spin backflow I^^'^'' that accompanies the 
nonequilibrium spin accumulation fi^ = ?i7Hcff. 



C. FMR-operated spin battery 

In the preceding subsection, we showed that precessing 
ferromagnets inject a spin current into adjacent conduc- 
tors via Ohmic contacts. In this subsection, we discuss 
how this opens the way to create a pure spin source ( "spin 
battery") by ferromagnetic resonance (FMR) and esti- 
mate the spin current and spin bias for different material 
combinations (Brataas et aL, . .2002.} . The spin source can 
be realized as a ferromagnetic layer at resonance with an 
rf field. Pure spin-current injection into low-density con- 
ductors should allow experimental studies of spintronic 
phenomena in mesoscopic, ballistic, and nanoscale sys- 
tems. The combination of a ferromagnet at the FMR in 
Ohmic contact with a conductor can be interpreted as a 
spin battery, with analogies and differences with charge 
batteries. For example, charge-current conservation dic- 
tates that a charge battery has two poles, plus and minus. 
A spin battery requires only one pole, since the spin cur- 
rent does not need to be conserved. Furthermore, the 
polarity is not a binary, but a three-dimensional vector. 

Devices made from metallic layered systems dis- 
playing the giant llBaibich et aL. .1983) and t unnel 
i|Mivazaki and TezukaLll995HMoodera et a/.Ul995(l mag- 
netoresistance have been proven useful for read-head 
sensors and magnetic random-access memories. Inte- 
gration of such devices with semiconductor electron- 
ics is desirable but difficult because a large conduc- 
tivity mismatch between magnetic and normal mate- 
rials is detrimental to spin injection ijSchmidt et all 
12000,') . Spin injection into bulk semiconductors has 
so far been reported only in optical pump- and- probe 
experiments (iKikkawa an d Awschal om. '1999*1 a nd with 
high- r esistance ferromagn etic injectors ( Fi edt;rling et all 
Il999t lOhno et all Il999l) or Schottky /tunnel barriers 
l)Monsma et all Il998t IZhu et all l200ll[) . In these cases, 
the injected spin-polarized carriers are hot and currents 
are small, however. Desirable are semiconductor devices 
with an efficient all-electrical cold-electron spin injection 



and detection via Ohmic contacts at the Ferm i energy, 
just as has been realized bv lJedema et a^J . f200ll for metal- 
lic devices. The spin battery discussed in this subsection 
is an alternative conceptual approach to accomplish such 
electrical spin injection. 

In the absence of spin-orbit scattering, the spin-battery 
idea can be understood in terms of the rotating-frame 
analysis of the previous subsection. A magnetization un- 
dergoing a circular precession around the z axis with fre- 
quency Lo induces a spin- imbalance density s = xa;z/7 
in the normal metal or semiconductor adjacent to the 
ferromagnet. The spin susceptibility x includes possi- 
ble effects of electron-electron interactions in the normal 
metal/semiconductor. The chemical-potential difference 
(|52|l between spin-up and spin-down electrons along the 
z axis, on the other hand, is universal for the FMR- 
operated spin battery with vanishing spin-orbit coupling. 
Spin-orbit scattering limits the efficiency of the spin bat- 
tery in real structures, as discussed in the following. 

The important parameters of a charge battery are the 
maximum voltage in the absence of a load, as well as 
the maximum charge current that can be drawn from 
it. In the following we discuss estimates for the analo- 
gous characteristics of the spin battery, closely following 
iBrataas et all I2OO2I To this end, one should solve the 
dynamic problem of spin pumping (|38|l at the F|N con- 
tact as a boundary condition for the spin diffusion equa- 
tions in the normal conductor. When the ferromagnet is 
thicker than the ferromagnetic coherence length H13|l (a 
few Angstr0ms in transition metals such as Co, Ni or Fe), 
the spin current (|38|l emitted into the normal conductor 
is determined by the mixing conductance g'^, since 
vanishes. With some exceptions such as f erromagnetic in- 
sulators l)Huertas-Hernando et all l2002[l. the imaginar y 
part of the mixing conductance is small |xia et a^J . |2002|) . 
see also Tabled and will be disregarded in the following. 
The spin current emitted into the normal metal is then 
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When the spin current (|54|) flows freely into the normal 
metal, the corresponding loss of angular momentum in- 
creases the (Gilbert) damping of the magnetization dy- 
namics, as discussed in detail in the next section. Eq. H54() 
is the maximum spin current that can be drawn from 
the spin battery. The rotating-frame analysis indicates 
that a nonequilibrium spin accumulation builds up in the 
normal conductor when the spin-flip relaxation rate is 
smaller than the spin-injection rate. This in turn induces 
a "backflow" spin current if' that opposes the total spin 
current Is ~ IP""^P+li"''. The component of the backflow 

spin current li*^^ parallel to the instantaneous magneti- 
zation direction m is canceled by an opposite flow from 
the ferromagnet when the FMR frequency and spin-flip 
scattering rate in the ferromagnet are much smaller than 
the characteristic spin-injection rate into the ferromag- 
net. The surviving component of if , Eq. I|2t)|) . is per- 
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pendicular to m: 

li°) = -|^mxM,xm. (55) 

We note that the mixing conductance in Eqs. H54|l and 
(|55|l ought to be renormahzed in layered structures, as 
discussed in Sec. III.BI 

The relation between the spin accumulation /x^ and the 
total spin current Is in a normal diffuse conductor is gov- 
erned by the spin-diffusion equation, see Sec. IfV.BI Let 
us consider an F|N bilayer structure with cross section 
S and thickness L, diffusion coefficient D and character- 
istic spin-flip time Tgt for the N-layer. T he so lution for 
the spin accumulation is simple when ^yD/uJ <^ L <C 
\/Dts{ = Asd , requiring uj :s> t^^ l)Brataas et alL^OO^ . 
The spin accumulation in the normal layer is then nearly 
uniform and time- independent for a steady ferromagnetic 
precession cone with angle with angle 6: 

t^s ^ z . (56) 

sm tl + 1] 

•q = Ti/rgf is a reduction factor expressed in terms of the 
"injection rate" Tj"^. When the conductance of a normal- 
layer slab of thickness ^ D juj is much larger than the 
contact conductance gj^, Tj"^ = g]}- j {JiM SV), where M 
is the normal-metal one-spin density of states per unit 
volume. In the opposite limit, t^^ — \/Duj/L. The 
former regime is relevant for an opaque interface, such as 
a Schottky or tunnel barrier. In the latter regime, ^ 1 
since we have assumed ojTst ^ 1 and L ^ ^/Dr^f. Large 
systems have a smaller injection rate since more states 
have to be filled. When Tgf ^ cxd, 77 = and Eq. itS^ 
reduces to the result of the rotating- frame analysis (|52|l . 
ISchmidt et oil I2OOOI realized that efficient spin injection 
into semiconductors by Ohmic contacts is difficult with 
transition-metal ferromagnets since virtually all of the 
applied potential drops over the nonmagnetic part and is 
unavailable for spin injection. The present mechanism of 
spin injection is current- rather then bias-driven and thus 
docs not suffer from this conductivity-mismatch problem. 

When the spin-relaxation rate is smaller than the spin- 
injection rate and the precession-cone angle is sufficiently 
large, sin^ d > -q, the spin bias (|56|l saturates at its max- 
imum value 1)52(1 . that does not depend on the material 
parameters. At resonance, sm.6 = h-ri/ {aHcs) (assum- 
ing for simplicity a constant effective field Hoff). For a 
dc field of H^s = 1 T, rf field /iif = 1 mT and damping 
a = 0.01, sin^ 9 = 0.01. In order to realize a battery with 
maximum spin bias (|52|l . a suppression factor rj < 0.01 
is thus required. Epitaxially-grown clean samples with 
low spin-flip rates are needed for spin batteries operating 
at small precession angles. The precession-cone angle 6 
in FMR experiments is typically small, but large angles 
could in principle be achieved by sufficiently intense rf 
fields and a soft ferromagnet such as permalloy or, for 
transient dynamics, by transverse magnetic-field pulses. 



The dc component of the maximum spin-current bias 
(|54|l is suppressed by sin^ 6 upon averaging over one pre- 
cession period. For operation as a spin current source, 
a large-angle precession is required as well. The current 
can be increased by a larger F|N interface area. 

A potential practical problem of large-angle dynam- 
ics is the FMR energy dissipation. For a total spin 
Sp — MsV/^ of a ferromagnet of volume V , dissipation 
power P — aSp^'^ sin^ 9 has to be efficiently drained in 
order to avoid excessive heating. On the other hand, 
possible undesirable spin-precession and energy genera- 
tion in the nonmagnetic parts of the system is of no con- 
cern for material combinations with different g factors, 
as, e.g., Fe {g = 2.1) and GaAs [g — —0.4), or when the 
magnetic anisotropy shifts the resonance frequency with 
respect to electrons in the normal metal. 

Purified elemental metals such as Al and Cu are 
suitable materials, since their spin-diffusion lengths are 
very long: In Cu, Asd ~ 1 /^m at low temperatures 
that is reduced by a factor of 3 at room temperature 
( Jedem^^ al. , 2001) , and similar numbers hold for Al 
( Jedema et al.. . .2002,) . Indirect proof of a spin accumula- 
tion pumped into Cu is provided by th e FMR measure- 
ment s of permalloy! Cu|Pt structures ijMizukami et al\ . 
i2002affl) discussed in Sec. llV.ljl Semiconductors have 
the advantage of a larger ratio of spin bias to Fermi en- 
ergy. In lightly n-doped GaAs, the spin-flip relaxation 
time can be ver y long: Ts = 10~^ s at n = 10^ ^ cm^"^ 
carrier density ijKikkawa and Awschalornl Il998|) . This 
favorable number is offset by the difficulty to form Ohmic 
contacts to GaAs, however. Large Schottky barriers ex- 
ponentially suppress the interface mixing conductance. 
InAs has the advantage of a natural charge accumu- 
lation layer at the surface that avoids Schottky bar- 
riers when covered by high-density metals. However, 
the spin-orbit interaction in a narrow-gap semiconduc- 
tor such as InAs is substantial, which reduces Tgf. In 
asymmetric-confinement structures, the spin-flip relax- 
ation rate is governed by the Rashba-type spin-orbit 
interaction that vanishes in symmetr ic quantum wells 
l|Engels et all Il997t iNitta et all 119971) The r emaining 
D 'yak onov-Perel scattering rate ijP'vakonov and PcreJ 
119711) is reduced in narrow quasi-one-dimcnsional chan- 
nels of width w due to waveguide diffusion modes by a 
factor of (Ls/w)^, where i s the spin-precession length 
l)Mal'shukov and Chad |2000^ , which makes InAs and its 
alloys potentially interesting materials for a spin battery 
as well. In pure Si, the spin-flip relaxation time should 
be very long, since spin-orbit interaction is weak. Fur- 
thermore, the possibility of heavy doping allows control 
of Schottky barriers. Si therefore also appears to be a 
good candidate for spin injection into semiconductors. 

The FMR-generated spin bias can be detected nonin- 
vasively via a ferromagnet that is connected to the nor- 
mal metal by a tunnel junction. A voltage difference 
pfis should be measurable between parallel and antipar- 
allel configurations of the analyzing magnetization with 
respect to the spin accumulation in the normal metal. 
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Here p = {G^ — Gi)/{G-[ + Gi) is the relative polariza- 
tion of the spin-dependent tunnel conductance G"' of the 
contact. For a fixed analyzing magnetization, the same 
voltage difference can be generated by reversing the static 
magnetic field. The spin current, on the other hand, can 
be measured via the drop of spin bias over a known re- 
sistive element with weak spin-fiip scattering. 

Spin pumping into the normal metal can also have 
consequences for the nuclei via the hy perfine interaction 
betw een electrons and nuclear spins ijKawakami et aO . 
I2OOII) . An initially unpolarized collection of nuclear 
spins can be oriented by a spin-polarized electron cur- 
rent, which transfers part of its angular momentum by 
spin-fiop scattering due to the hyperfine interaction. A 
ferromagnetically-ordered nuclear-spin system, in turn, 
is felt by the electrons as the Overhauser magnetic field 
which indu ces an equilibrium spin density Sq in the nor- 
mal metal ()Overhauseilll95 3'). The spin-flop scattering 
can be included in electron-spin dynamics by adding an 
additional term to the total spin current, = JP"™? + 

to account for the exchange of angular momen- 
tum between electrons and nuclei (Br ataas et ai, 2002}. 
The nuclear spin polarization increases with the spin ac- 
cumulation and decreasing temperature. In bulk GaAs, 
the nuclear magnetic field^is iJ„ = 5.3 T when the nuclei 
are fully spin polarized ifPaeet et alLll99T^ . which should 
occur at thermal energies sufficiently smaller than the 
FMR frequency, as can be understood from the rotating- 
frame analysis. 

In conclusion, the spin battery is a source of spin, just 
as a conventional battery is a source of charge. Estimates 
of its performance for bilayers of metallic ferromagnets 
with either normal metals or doped semiconductors sug- 
gest that it is a feasible concept. 



pumping-induced contributions. Suppose initially that 
the normal metals act as ideal reservoirs in a common 
thermal equilibrium that are perfect spin sinks for the 
pumped spin currents, so that I^";' = 0. This model is 
valid when the spin current I^""^^ is completely drained 
by massive and highly conductive reservoirs or by a ma- 
terial with effective spin-flip processes that prevent any 
spin-accumulation build-up and backscattering into the 
ferromagnet. Using Eq. 
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Since this spin current is polarized perpendicularly to m, 
the torque H12|l becomes simply r = — mxIsXm = — I^. 
Substituting this into equation H14() . we recover the LLG 
equation JHJl with a new effective damping constant, aoff, 
and gyromagnetic ratio, 7off, deflned by 
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The spin pumping thus affects FMR experiments as a 
shift of the resonance magnetic fleld via .4J''' -I- A'^^^, 
whereas y^J^ + A'^-^ increases the relative resonance 
linewidth. Al^ defined by Eq. (|39|l is positive, since, by 
unitarity of the scattering matrix, it can be rewritten as 
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IV. GILBERT-DAMPING ENHANCEMENT 



A. Ideal spin sinks 



In section IIIII we showed that a moving magnetiza- 
tion emits spins into adjacent nonmagnetic conductors. 
This effect is necessarily associated with an energy loss 
for the ferromagnet, as explained in Sec. IIII.51 We now 
discuss spin pumping as a source of viscous damping of 
the magnetization dynamics in thin films or small par- 
ticles. Under quite general conditions, this damping is 
consistent with the Gilbert phenomenology (|SJ). 

Consider a ferromagnet sandwiched between two nor- 
mal metals labeled by / = L, R. By conservation of 
angular momentum, spins ejected out of the ferromagnet 
exert a transverse reaction torque H12|l . The total spin 
current 
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that determines the spin torque, i.e., the additional 
term (|14ll in the LLG equation, contains bias- and spin- 



The spin pumping by a moving magnetization therefore 
cannot reverse the sign of the effective damping parame- 
ter (without reversing the sign of the gyromagnetic ratio) , 
as is required by the LLG phenomenology, see Sec. II. CI 

Enhanced Gilbert damping leads to a broader reso- 
nance absorption peak (|ll|l of an rf magnetic field. By 
the fiuctuation-dissipation theorem, this should be man- 
ifest in increased fluctuations of the magnetization in 
thermodynamic equilibrium. This additional magnetiza- 
tion noise can arise from the torques exerted by Johnson- 
Nyquist spin-current fluctuations that are exchanged be- 
tween t he ferromagnet and reservoirs or other spin sinks. 
Indeed, iForos et~al\ . l2005al proved explicitly that mag- 
netic noise caused by the spin-current fluctuations is con- 
sistent with the dissipation power predicted by the spin- 
pumping theory for a monodomain ferromagnet in con- 
tact with ideal spin sinks. 

[Mizukami et ai, 2001a h, measured room-temperature 
FMR linewidths of sputtered N| permalloy (Ni8oFe2o, 
Py)|N sandwiches, and discovered systematic trends in 
the damping parameter as a function of Py-film thick- 
ness d for different normal metals N. Their data for the 
Gilbert parameter G, Eq. lO, are shown by symbols in 
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Fig. [S| In the following, we wish to compare these mea- 
surements with theoretical predictions based on Eq. (|60|l . 
For the Py films, d > 2 nm > A^c, so that t^^, t'^^ w 
and A^^, Eq. H39I) . is simply the mixing conductance g^^ 
of the Py|N interface, gj^ ^ gl^ for the interfaces listed 
in Tabled so that 

ffTi « gn . (62) 

Eq. H59|) then reduces to jcS ~ 7 and Eq. H60|l becomes 

Geff«G+^^. (63) 

Moreover, for the parameters in Tabled 

gl^^gf. (64) 

As explained by IZwierzvcki et alL l2005l Eqs. H62() and 
(|64|l are good approximations for intermetallic inter- 
faces because of the large phase differences between spin- 
up and spin-down reflection coefflcients (see, however, 
Sec. IVii.A.2l for ultrathin magnetic films, d < Age)- We 
will assume the validity of these approximations in much 
of the review. Sharvin conductances for different normal 
metals (Zwierzvcki et ai, 2005) are listed in Table ITU 
The data in Fig. |S1 are used to extract gl^ according 
to Eq. (|^ . A possible d-dependence of the bulk G 
is disregarded in our analysis. The G's in Fig. [S] are 
close to the previously meas ured bulk Gilbert damp- 
ing constant s of p ermalloy l)Bastian and Billeil Il976t 
iPatton et all Il975l) . We see that the extracted mixing 
conductance for the Pd|Py interface is similar to that ex- 
pected from the Sharvin conduc tance of Pd, while that of 
Pt|Py is about 1.5 times larger. iMizukami et a^J . l2001albl 
also reported FMR measurements on Ta|Py sandwiches 
that had only a small damping enhancement with respect 
to the bulk value. 
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TABLE II Sharvin conductances (in units of quantum chan- 
nels per nm'^) for bulk fee Cu, Pd, Pt and bee Ta, den- 
sity of states A/" at the Fermi level (in units of states per 
Ry, atom, spin) , and Stoner enhancement factor. From 
IZwierzvck i et all l2005l . Typical values of the exchange- 
correlation integral, 7xc, were taken from iGunnarsson. .1976. : 
|janaklll977l . 

According to the earlier arguments, the absence of 
a significant thickness dependence of the damping in 
the Cu|Py system could correspond to an unrealistically 
small conductance glK The explanation should however 
be sought in the long spin -flip relaxation times in clean 
Cu ijMeservev and Tedrowlll978^) : Since the accumulated 




i> ♦ ♦ 4^ <>■ 

q\ , \ , \ 

2 4 6 8 10 

d [nm] 

FIG. 5 Symbols represent the data points derived from FMR 
experiments on N|Py|N sandwiches with N=Pt, Pd, or Cu 
l|Mizukami et all l2001allbll . The lines show fits by G -I- a/d. 
G = 0.97, 1.06, and 1.31 x 10** s'^ , where following Eq. (15^ . 
a corresponds to gl^/S = 0.2, 15.3, and 25.8 nm~^ for Cu, 
Pd, and Pt, respectively , taking a (^-factor of 2.1 for Py films 
iMizukami et ai.ll2001allbt) . The horizontal axis uses a recip- 
rocal scale. 



spins drive a diffusion spin current opposite to the pump- 
ing current, the additional damping vanishes in clean 
Cu. On the other hand, Pd and especially Pt (which 
is below Pd in the periodic table) have much larger spin- 
relaxation rates, that can be rationalized by the higher 
atomic number and complex Fermi surfaces (Pd appears 
to also have an additional spin-decoherence mechanism 
du e to spin-flip scatter ing at magnetization fluctuations, 
see lForos et aZ.L l2005bj) . We note that the spin-fiip effi- 
ciency of a dirty normal metal is dete r mined by defects 
and impurities as well. iLubitz et all l2003l reported a 
significant damping contribution that scaled as 1/d as 
a function of the Fe-layer thickness d in polycrystalline 
Cu|Fe|Cu sandwiches in contrast to the experiments by 
IMizukami et all l2001bl on Cu|Py|Cu. This was inter- 
preted in terms of a larger Cu spin-relaxation rate in the 
Fe|Cu as compared to the Cu|Py system. The effect of 
spin flip in the normal metal on the excess damping in 
multilayers is treated quantitatively in the next subsec- 
tion. 

Ilngvarsson et alL l2002l also carried out FMR studies 
on Py films sandwiched by the normal metals Pt, Nb 
(which is above Ta in the periodic table) , and Cu, as well 
as by insulators. Pt|Py combinations displayed a signif- 
icantly stronger damping than other structures. They 
also identified an additional thickness-dependent damp- 
ing correlated with the disorder in the Py films that was 
not reported bv lMizukami et an . l2001albl . Damping pro- 
cesses that are intrinsic to the ferromagnetic layer cannot 
be addressed by the present theory, however. 

iLagae et aLLl2005l studied pulsed dynamics of a thin Py 
layer in contact with Cu or Ta layers. In addition, they 
investigated the role of a second, pinned ferromagnetic 
layer, attached to the free layer via a Cu spacer or a 
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tunnel barrier. The measured Gilbert damping was found 
to be consistent with the spin-pumping picture. We will 
discuss the effect of the second ferromagnetic layer in 
Sec. UvTH 

Finally, we comment on the role of electron correla- 
tions in the spin-pumping formalism, see also Sec. I VII .Dl 
In order to calculate the mixing conductance, the scat- 
tering matrix h a s bee n obtained by iXia e t al^, l2002 t 
IZwierzvcki et a/.l . l2005l for the effective potential from 
Kohn-Sham density-functional theory. This potential is 
calculated self-consistently and includes electron-electron 
interaction effects via the Hartree and an exchange- 
correlation potential in the local-spin-density approxi- 
mation. The interface parameters when computed self- 
consistently and nonperturbatively include the magnetic 
moments induced in the normal metal by the proximit y 
to the ferromagnet (also discussed bv ISimanekl l2003(l . 
For Cu and Au normal-metal films, this effect is very 
small due to the small densities of states at the Fermi 
energy. The question arises whether larger effects may 
be expected for materials such as Pd and Pt with a large 
Fermi-level density of states (see Table lll)l . These metals 
are known to be "almost ferromagnetic" with a Stoner- 
enhanced spin susceptibility x/xo = (1 — J^Ixc)^^ , also 
listed in Table UTI So, can the increased damping of the 
magnetization dynamics in con tact with thin layers of Ta , 
Pd and Pt as compared to Cu ijMizukami et aLLl2001albl) 
connected to their larger Fermi-level densities of states? 
In the spin-pumping theory, the quantity that governs 
the damping enhancement is the mixing conductance. 
For metallic structures, the band-structure calculations 
give results that are close to the Sharvin conductance 
of the normal metal, Eq. (|64|l . thus not sensitive to the 
magnetic moments induced in the proximity to the fer- 
romagnet. The Sharvin conductance for Cu, Ta, Pd and 
Pt is seen to change less than J\f. More significantly, 
the trend does not correspond to that observed exper- 
imentally for the damping enhancement: The Sharvin 
conductance is maximal for Ta, not Pd with the highest 
spin susceptibility, supporting our earlier conclusion that 
the spin dissipation in the normal conductors is crucial 
for the enhanced magnetization damping, as is explored 
in more detail below. 



B. Diffuse systems 

In the previous subsection, we concentrated on the ex- 
treme situation in which the normal-metal layer is an 
ideal spin sink that immediately absorbs the injected spin 
current IP"™p, either by relaxation through spin- flip pro- 
cesses or the absence of back scattering. Disregarding 
the total spin current through the contact, Eq. ( 
1^ . reduces to Ig = 1^"™^. In general, we have to take 
into account the spin accumulation in a diffuse normal 
metal that drives a spin current back into the ferro- 
magnet. As a technical note, we recall that in this regime 
the drift correction H34|l should be applied to the conduc- 



tance parameters entering Eqs. (|25|l . (|26|l . and (|38|l . 

As the simplest example, we first discuss an F|N bi- 
layer, in which a magnetization precession generates a 
pure spin but no charge current. The spin accumulation 
or nonequilibrium chemical potential imbalance fJ-six) 
[similar to Eq. (|l(j|l . but with spatial dependence] in the 
normal metal is a vector with magnitude that depends on 
the distance from the interface x, < x < L, where L is 
the thickness of the normal-metal film, see Fig. |B| When 
the ferromagnetic magnetization steadily precesses (in an 
easy plane) around the z axis, m x rh and the normal- 
metal spin accumulation fJ-six) are oriented along z, as 
depicted in Fig. We will see in the following that the 
instantaneous /x^ is always perpendicular to m even in 
the case of a precessing ferromagnet with time-dependent 
rotation axis, as long as the precession frequency w is 
smaller than the spin-flip rate t^^ in the normal metal. 




FIG. 6 Schematic view of an F|N bilayer in which the mag- 
netization direction m(t) of the ferromagnet F rotates in an 
easy plane, pumping a spin current IP"™p into the adjacent 
normal-metal layer N. The N layer here is not an ideal reser- 
voir but rather a film of the same cross section as the magnetic 
layer F. The spin pumping builds up a position (a:) dependent 
spin accumulation in N that either relaxes by spin-flip scat- 
tering or drives a spin current back into the ferromagnet as 



I 



(0) 



In frequency domain, a spin accumulation diffuses in a 
normal metal according to Ijohnspn and Silsbec , 1988) 



Tsf 



(65) 



where D = VpT/3 is the diffusion coefficient (in three 
dimensions) in terms of the Fermi velocity vp, and the 
transport mean free time r. It is assumed here that the 
frequency uj is smaller than the scattering rate t~^. The 
ratio of the momentum to spin-flip scattering time is an 
important parameter: 
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Tsf 
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Eq. H65|) holds when e <C 1, i.e., the spin- flip relaxation 
may be treated perturbatively. The boundary conditions 
are given by the spin current at the F|N interface x = 
and vanishing of the spin current at the outer boundary 
X — L: 



X =L : dxfJ-s = , 



(67) 
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where S is the interface area and Af the normal-metal 
one-spin density of states per unit volume. The solution 
of Eqs. (EHJ, EI) is 



cosh [k{x — L)] 2 
sinh kL hAfSDn 



(68) 



where k — yT+Tcj?^, recalling that Asd = VDr^f 
is the spin-diffusion length in the normal metal. In 
Sec. nil. CI we used similar arguments to calculate the 
spin accumulation generated by the precessing magne- 
tization. Whereas in Sec. IIII.CI the magnitude of the 
spin accumulation and its relevance for spintronics are 
discussed, we focus here on the effect of the spin accu- 
mulation on the ferromagnetic magnetization dynamics. 

We assume in the following that the precession fre- 
quency CO is smaller than the spin-flip relaxation rate, 
w <C ^ , so that K « A^^ . For a typical effective field of 
1 T, w ~ 10^^ s~^. The scattering rate corresponding to 
a mean free path of A ~ 10 nm is ^ 10^"' s'^. Con- 
sequently, for such metals, uj <C t^^ requires e > 10~^ 
This condition is easily satisfied ijMeservev and Tedrowl 
Il978f) for imp urities with higher atomic num bers Z [since 
e scales as Z4 lib rikosov and Gor 'kOTlll962D ]. The high- 
frequency limit a; > Tg^^, on the other hand, is relevant 
for systems with weak spin-flip scattering in the normal 
metal, and has been discussed in the context of the spin- 
battery concept in Sec. IIII.CI We will see that a sizable 
Gilbert-damping enhancement requires a large spin-flip 
probability e > 10^^, thus uj t^^ , unless the frequency 
is comparable with the momentum scattering rate in the 
normal metal. [The latter regime lies below our treat- 
ment based on Eqs. (ESJ, (|57|) .] 

It is convenient to define an effective energy-level spac- 
ing of the states participating in the spin-fiip scattering 
events in a thick (i.e., L ^ Asd) normal- metal film as 
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The spin-diffusion length (which sets the scale for spin- 
accumulation penetration into the normal layer) written 
in terms of the relevant scattering times is 



\ /TTsf 
Asd — Vf\I ■ 



(70) 



The spin-accumulation-driven spin current li*^^ through 
the interface is obtained by substituting ^ig{x — 0) from 
Eq. ^ into Eq. ^ to give 



(71) 

where the spin current returning into the ferromagnet is 
governed by the "backflow factor" /3: 



TsfSsd/h 

tanh(_L/Asd) 



(72) 



The last equality in Eq. H71|l is obtained from Eq. H26() 
by assuming that the nonequilibrium transport is lim- 
ited to spin currents that are polarized normal to the 
magnetization. This is allowed in the adiabatic regime 
when UJ <^ Tsf where the average angular momentum of 
the pumped and backscattered spins is at a given instant 
perpendicular to the magnetization. Consequently, as 
long as the ferromagnetic spin-relaxation length is larger 
than the transverse spin-coherence length Age, we may 
disregard spin relaxation inside the ferromagnet. (3 is 
given by the ratio between the energy-level spacing of 
the normal-metal film with a thickness Lgf = min(L, Asd) 
and the spin-flip rate. When the normal metal is shorter 
than its spin-flip diffusion length, L <C Asd, P TsfS/h, 
where S = {J^SL)~^ is the energy-level splitting. In 
the opposite regime of thick normal metals, L ^ Asd, 

P Ts[Ssd/h. 

By solving Eq. H71(l . the total spin current Is can be 
expressed in terms of the pumped spin current IP"™p, 
Eq. , to finally obtain 

^ , T I dm , 1 1 dm \ , , 

^- Re^I^.mx_+W:^- , (73) 
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which has the same form as the original Eq. (|38|) . but 
with effective "spin-pumping efficiencies" A_ 



Ti. 

off- 
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Here, i?sd = TsfSsd/h is the dimensionless resistance [in 
units of (2e^//i)^^] of the normal-metal layer of thick- 
ness Asd, which follows from Eq. H69(l and the Einstein 
relation a — e^DAf between conductivity a and diffusion 
coefficient D. When L ^ Asd, the effective spin pumping 
out of the ferromagnet is governed by the mixing con- 
ductance of the F|N interface in series with a diffusive 
normal- metal film of thickness Asd- 

By accepting the approximation (|62|l . we disregard gj'^ 
in the remainder of this section. Im^^Jg then also van- 
ishes and upon substitution into Eq. (|14ll the damping 
torque due to the spin current Is obeys the Gilbert phe- 
nomenology. The effective Gilbert-damping parameter in 
the diffuse model reads [cf. Ea. l|^ ]: 
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sd 
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The prefactor on the right-hand side of Eq. H75|l reffects 
the reduction effect on the Gilbert damping that is caused 
by the spin-diffusion back into the ferromagnet. This 
correction has been disregarded in Sec. IIV.XI where the 
normal metals were considered to be ideal spin sinks. 
Because spins accumulate in the normal metal polar- 
ized transversely to the ferromagnetic magnetization, the 
spin-accumulation-driven transport across the F|N con- 
tact is governed by the spin-mixing conductance (just like 
the spin-pumping current). The absence of spin accumu- 
lations or currents polarized coUinear to the magnetiza- 
tion (in the limit u) <C t^^) explains why the diagonal 
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components of the conductance matrix do not enter 
Eq. CZSJ. 

The numerical values of the parameters in Eq. H75|) 
can be estimated by the free-electron model for the nor- 
mal metal, for which Af — kp/{TThvF)- Using Eqs. H70() 
and (Enil, we find R^^ = h/{6sdTsf) = 4:^/7/3g^, where 
gf^ = Skp/{ATr) is the dimensionless Sharvin conduc- 
tance (number of transport channels) of the normal 
metal. Using approximation H64I) for the mixing con- 
ductance gj^, and Eq. H33|l to correct for the drift con- 
tribution, we obtain gl^ ~ 2(ji^^. We thus estimate 

1-1- [Vitanh(L/A,d)] 

A significantly increased damping therefore requires a 
high spin-flip to spin-conserving scattering probability, 
e > 10^^, and sufficiently thick normal layers, L > Asd 
{L > Asd alone is not a sufficient condition, however, 
contrary to what has sometimes been assumed in the lit- 
erature) . Poor spin sinks that do not modify the magne- 
tization dynamics are characterized by a large denomina- 
tor in Eq. JTCJl. Clean light metals with atomic number 
Z < 50, such as Al, Cr, and Cu, as well as heavier met- 
als with only s electrons in the conduction band, such 
as Ag, are less effective spin sinks because of a relatively 
small intrinsic spi n-orbit c oupling, typic ally in the range 
of £ < 10~^ (Bergmannl ^1982; Mcscr vev and Tedrowl 
ll978tlYang et adll994D .' Heavier elements with Z > 50 
and p OT d character of the conduction electrons, such 
as Pd, Pt, and Pb, are effective spin sinks with much 
larger e > 10~^ llMeservev and Tedrowl fl978j) . Au is in- 
termediate in e ijChiang et alU2Q04\ . presumably due to 
the s character of its outer orbital. These conclusions 
explain the hierarchy of t he observed Gilbert dam ping 
enhancement observed by iMizukami et all l2001albL see 
Sec. IIV.AI By doping a small-e matrix with high-Z or 
magnetic impurities (e.g., Cu lattice with Pt atoms), a 
good spin sink can be created. We conclude that the 
damping parameter (and thus the switching speed) of 
thin ferromagnetic films can be engineered flexibly by 
coating them with suitable normal metals. 

In the limit of a large ratio of the spin-flip to momen- 
tum scattering, e 1, the spin-diffusion equation and, 
consequently, Eq. (|7S|l do not hold. In this regime, the 
normal metal very effic iently relaxes the sp ins injected 
from the ferromagnet. iForos "oil 120050 investigated 
this scenario for Pd (in contact with Fe) that appears to 
have a spin-relaxation rate, arguably by paramagnons, 
that is faster than the bulk momentum-scattering rate. 
When thicker than the spin-decoherence length vpTsi, 
otherwise ballistic Pd films should thus be good spin 
sinks. The ferromagnetic relaxation by ideal spin sinks 
is determined by the bare spin-mixing conductance g^^, 
as discussed in Sec. IIV.AI The latter generally pro- 
vides an upper bound for the magnitude of the addi- 
tional Gilbert damping. However, in the present theory, 
spin-orbit coupling is treated phenomenologically in the 



diffuse-transport regime only. Strong spin-orbit coupling 
immediately at interfaces, for example, requires general- 
ization of spin-pumping and circuit theories beyond the 
scope of this review. 

Infinite vs vanishing spin-flip rates in the normal metal 
are two extremes for the F|N bilayer dynamics. In the 
former case, the damping parameter Gcfr is enhanced, 
and in the latter case unmodified. Both limits are ex- 
perimentally accessible by using Pt as a good or Cu a s a 
poor spin sink, as shown bv IMizukami et all . 1200 lalbl for 
N|Py|N sandwiches. 

Experiments on Py where both Cu a nd Pt were com- 
bined in a Py|Cu|Pt heterostructure ijMizukami et all 
_2002£^ provide a novel method to study spin-diffusion 
in the central Cu layer. The measured room- 
temperature magnetization damping in Cu|Py|Cu(L) 
and Cu|Py|Cu(L)|Pt structures as a function of Cu-film 
thickness L is shown by circles in Fig.EI The experiments 
can be understood by a slight extension of the diffusion- 
theory-based discussion for F|N bilayers. The damping 
is governed by the angular-momentum loss of the ferro- 
magnet through the normal-metal compound Ni|N2, as 
schematically shown in Fig. |S| Once injected into Ni, 
spins are either scattered back into the ferromagnet or 
relax in Ni or N2. When there is only a very weak spin- 
flip scattering in Ni compared to N2, the spins have to 
diffuse through Ni before they can relax in the second 
normal- metal layer N2 . For simplicity, we model N2 as an 
ideal sink that instantaneously relaxes incoming spins, as 
is appropriate for Pt. The analysis below illustrates that 
the ferromagnetic magnetization damping as a function 
of L provides important information about spin transport 
through Ni and the N1IN2 interface. The Cu substrate 
on the o ther side of the permalloy fllm l)Mizukami et all 
l2002albf l can be disregarded since it is a poor spin sink. 
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FIG. 7 Gilbert-damping constants (circles) in 
Cu|Py(3 nm)|Cu(L)|Pt (5 nm) and Cu|Pv(3 nm)|Cu(L) 
structures measured by IMizukami et al 1 l2002 a'>b! as a func- 
tion of Cu-layer thickness L. Solid lines are the theoretical 
results according to Eq. H80^ with parameters discussed in 
the text. The horizontal axis uses a logarithmic scale. 
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FIG. 8 Similar to Fig. El but the normal-metal system is 
an N1IN2 bilayer. Ferromagnetic magnetization precession 
pumps spins into the first normal-metal layer Ni. The spin 
accumulation in Ni either flows back into the ferromagnet F 
as spin current I^'J', relaxes in Ni, or leaves via the second 
normal-metal layer N2 as spin current I^j' • The spin accumu- 
lation in the ideal spin sink N2 vanishes. 



Similar to Eqs. H67|l . the boundary conditions for the 
spin accumulation in the normal metal Ni are now: 



x = : 5:eMs = 



hAfSD 
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(77) 



Isl and Is2 are the total spin currents through the left 
(x — 0) and right (x — L) interfaces, respectively. 
1,1 = IP™P + cf. Eq. GH) in the previous sub- 
section, is the sum of the pumped (|38|l and the spin- 
accumulation-driven (|26|) spin currents. 1,2, on the other 
hand, is entirely governed by the Ni— >N2 diffusion 



(78) 



Setting g = decouples the two normal-metal systems 
and reduces Eq. (|Sn|l to Eq. l(75)l. the da mping coefficient 
of the F |Ni bilayer. In the experiment l)Mizukami et all 
|2g02a'b), the permalloy thickness d = 3 nm is fixed and 
the Cu film thickness L is varied between 3 and 1500 nm 
as shown by the circles in Fig. [7| The theoretical result 
(|Hn|l is plotted in Fig. H for comparison, using the fol- 
lowing parameters: the bulk damping G = 0.7 x 10^ s~^ 
l)Bastian and Billeil ITqtI iPatton et al\ . \l97^ . the spin- 
flip probability e = 1/700 and the spin-diffusion length 
Asd = 250 nm for Cu (which correspond to the mean free 
path A — VSeAsd = 16 nm), in sa tisfactory agreement 
with values reported i n literatu re llJedema et all l200lt 
iMeservev and TedrowL 11978: Ya ng et aZlll994|) . dl^/S = 
16 nm~^ extracted from the experimental a ngular mag- 
netoresistancc of Py|Cu l|Bauer et QiJ . l200'3bi) . and g/S — 
35 nm^^ for the Cu|Pt contact. This g corresponds to 
the bare one-spin conductance g'^'^ / S = 16 nm~^, which 
is close to the Sharvin conductance of Cu, see Table 
Figure H shows a satisfactory agreement (within the ex- 
perimental error) between experiments and theory. This 
proves the diffusive nature of spin transfer in the Cu 
spacer. Whereas the detailed mechanism for spin injec- 
tion (relaxation) at the Py|Cu (Cu|Pt) interface cannot 
be deduced directly, the agreement on absolute scale ob- 
tained with parameters taken from other sources strongly 
supports the spin-pumping picture. 

It is illuminating to discuss Eq. (|80|l in the limit of 
vanishing spin-flip in the spacer layer Ni. Recalling the 
definitions for Xgd, Eq. (|70|l . and Ssd, Eq. (|69|l . and tak- 
ing the limit Tgt — > 00, we find that Eq. (|80|l reduces to 
Eq. (|63|l . only with 2g^^ (where the factor of two cor- 
responds to two F|N interfaces in the N|F|N trilayers) 
replaced by gig. 



where g is the effective spin conductance of the N1IN2 
interface: 
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Here, we corrected the bare single-spin resistance 
V5nT|N2 of the all-normal interface for the drift effect by 
subtracting the Sharvin contribution on the Ni side only, 
because the layer N2 is assumed to be an ideal spin bath. 
Solving the diffusion equation H65() with the boundary 
conditions l|77(l . we can find the spin current I^i, which 
is transferred to the magnetization as a torque, as dis- 
cussed in the preceding discussion for the bilayer. The 
Gilbert damping enhancement due to the spin relaxation 
in the trilayer is then found to be 
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where Asd and Rsd parametrize the spin-diffusion in Ni, 
gl^ is the renormalized mixing conductance (|33|l of the 
F|Ni interface, and g is the spin-transfer conductance 
(EHl of the N1IN2 interface. 
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(81) 



where R^i — {2e'^ /h){L/ Sa) is the dimensionless resis- 
tance of the Ni layer with conductivity a. The right-hand 
side of Eq. (|81|l is simply the inverse bare mixing conduc- 
tance of the diffuse Ni spacer in series wit h its two in- 
terface s, one with F and the other with N2 ijBauer et all 
l20031j) . In particular, when Ni is thick enough, the total 
mixing conducta nce qli is limited by the s pacer sepa- 
rating F and N2 l|Brataas et all l2000l I200H) . The spin 
pumping into layer Ni with subsequent spin-conserving 
diffusion and then spin absorption by the ideal spin sink 
N2 is thus equivalent with spin pumping (|38|l across an 
effective scatterer separating the ferromagnet F from the 
ideal spin sink N2. 

The general trends in Fig. [7| can be understood as fol- 
lows. Since Cu is a poor spin sink, a Py|Cu contact with 
a single Cu film only weakly increases the damping for 
all thicknesses. This enhancement saturates at L ^ Asd 
and vanishes in the limit L ^ Asd- When a Pt film, 
a very good spin sink, is connected to the bilayer and 
the Cu spacer is thinner than the transport mean free 
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FIG. 9 Gilbert damping measured in CujPy(3 nm)|Cu(L)|Pt(2 nm) and Cu|Py(3 nm)|Cu(L) structures for several Cu-layer 
thicknesses L at various temperatures. Solid lines are fits based on the spin-pumping theory. Adapted from lYakata et an.l20nfiL 



path, L ^ A, the spin accumulation is uniform through- 
out the Cu. The spin pumping is then partitioned: One 
fraction of the pumped spins is reflected back into the 
ferromagnet, while the remainder is transmitted to and 
subsequently relaxes in the Pt layer. Their ratio equals 
the ratio between the conductance gl^ of the Py|Cu and 
the conductance g of the Cu|Pt interfaces, and is of the 
order of unity. Since a significant fraction of the spin- 
pumping current is dissipated in Pt, a large magneti- 
zation damping is achieved. When L is increased, less 
spins manage to diffuse through the Cu spacer, and, in 
the limit L ^ Asd, the majority of the spins scatter back 
into the ferromagnet before sensing the presence of the 
Pt layer. In the intermediate regime, the spin pumping 
into the Pt layer decays algebraically on the scale of the 
transport mean free path and exponentially on the scale 
of the spin-diffusion length. 

The temperature dependence of the Gilbert damp- 
ing in such t rilayers has been m easured recently by the 
same group l)Yakata et adl2005|) . see Fig. El The spm- 
diffusion length extracted from the data by means of an 
analysis similar to the one presented here increases as 
the temperature is lowered, indicating a reduced role of 



phonon scattering at low temperatures. The good agree- 
me nt with result s frorn electrical-transport experiments 
by iJedema et al\ . l200ll establishes FMR experiments as 
an important tool to measure spin transport in magnetic 
heterostructures. 

The dependence of the damping on the Cu-layer thick- 
ness L in the Cu|Py|Cu(L)|Pt multilayers reflects the 
amount of accumulation in the normal metals. This spin 
accumulation, in turn, indicates that an excited ferro- 
magnet (as in the FMR experiment discussed here) trans- 
fers spins into adjacent nonmagnetic layers according to 
the spin pumping lp!H|l . The concept of the spin battery 
discussed in sect ion lin . C I relies on this effect. 

C. Enhanced Gilbert damping in spin valves: 
First-principles calculations vs experiment 

In F|N|F structures, the presence of two ferromagnetic 
layers can make damping possible for each individual 
layer even in the absence of spin-flip relaxation in the 
system (^erger, 199^. The point is that if one ferromag- 
net is excited while the other is static, the latter acts as a 
sink for transverse spin currents pumped by the former. 
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In this subsection, we are assuming a sufficiently tfiick or 
disordered normal spacer, so that it does not support any 
persistent spin currents or, in other words, the static ex- 
change interaction between the magnetic films vanishes, 
see Sec. II.BI Static exchange will be taken into account 
in Sec. in discussing coupled dynamics of two or more 
magnetic films. If the F|N|F magnetic structure is weakly 
excited from a collinear equilibrium state, then, defining 
the effective (complex-valued) "spin-pumping efficiency" 
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(82) 



we can summarize an analysis similar to that of the 
preceding section as follows. The total spin current Is 
through the normal spacer consisting of the pumped and 
backflow components, Eq. (|36|l . is given by the right-hand 
side of Eq. (|38|1 with A^-^ replaced by -4p|i^|p where m 
is the magnetization of the excited film (assuming one of 
the ferromagnets is static). Effective damping and gy- 
romagnetic ratio are then given by Eqs. I|59() and H60() 
with A}p\^\-p substituted for A''^ +A'''^. In Eq. a is 

the conductivity of the N spacer, L is its thickness, and 
S is the area of the trilayer. In the spirit of the theory 
discussed in Sec. III.BI Eq. (|S^ requires that the spacer 
(in series with the two interfaces) is diffuse. The trans- 
mission mixing conductance t^^ is disregarded, assuming 
sufficiently thick magnetic films, d ^ Age, or insulating 
substrate and cap for the F|N|F trilayer. Adding inverse 
spin-mixing conductances in series with the diffuse-spacer 
resistance in Eq. reflects partitioning of the pumped 
spin currents between the two magnetic layers, having 
di sregarded here spin relaxation in the spacer. 

lUrban et alV l200ll reported room-temperature (RT) 
observations of increased Gilbert damping for a system 
consisting of two epitaxially-grown Fe layers separated 
by a Au spacer layer. The complete structures were 
GaAs| (8,11, 16,21,31)Fe|40Au|40Fe|20Au(001), where the 
integers represent the number of monolayers (ML's), and 
the samples differ in the thickness of the thinner Fe film. 
The i nterface magnetic anisotropics allowed Urban et aL . 
1200 ll to separate the FMR fields of the two Fe layers 
with resonance-field differences that can exceed 5 times 
the FMR lincwidths. Hence, the FMR measurements for 
thinner F layer can be carried out with a nearly static 
thick layer: The FMR linewidth of the thin F layer in- 
creases in the presence of the second layer. The differ- 
ence in the FMR linewidths between the magnetic bi- 
layer and single-layer structures is nearly inversely pro- 
portional to the thin-film thickness d, suggesting that 
the additional damping occurs due to its F|N interface. 
Secondly, the additional linewidth is linearly dependent 
on microwave frequency for both the in-plane (the sat- 
uration magnetization parallel to the film surface) and 
perpendicular (the saturation magnetization perpendicu- 
lar to the film surface) configurations, strongly implying 
that the additional contribution to the FMR linewidth 



can be described str ictly as an interface Gilbert damping 
llUrban 61^ oIl . l200lD . 

The magnetization of the thin ferromagnetic layer pre- 
cesses in the external magnetic field, while the other 
static magnetic layer acts as a spin sink. No modification 
of the damping coefficient was measured for GaAs|Fe|Au 
structures without a second Fe layer. The latter find- 
ing is consistent with the prediction given by Eq. H75|) 
in the L <C Ash lim it, well fulfilled for thin Au films of 
lUrban et a/] . l200lL In the presence of the second Fe layer, 
Eq. should be used: Neglecting Im^p|j^|p leads to 
7cff — 7 and the damping enhancement 
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where a ~ 0.0046 is the dimensionless da mping measured 
for a single Fe layer. Using 7 = 2.1^b/^ l)Heinrich et al\ . 

and the values of the interface and Sharvin con- 
ductances from Table m Eq. (|83|l is compared with the 
experimental data in Fig. 1101 for various assumptions 
about (7 in l|82|) . In the low-temperature limit and ne- 
glecting the residual resistivity of the Au layer, ct ^ 00, 
Eq. H83|) yields the solid line which is seen to overesti- 
mate the damping enhancement compared to the mea- 
sured results. Using finite values of a will lead to lower 
and indeed, it was found experimen- 
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tally ijHeinrich et aL . 2003b) that lowering the tempera- 
ture (increasing the conductivity) increases the damping 
by as much as about 20% (open circle in Fig. [Tn|l . If 
one uses the RT conductivity due to phonon scattering 
in crystalhne bulk Au, cTph = 0.45 x 10^ Q~^m~^, the 
dashed line is obtained which, as expected, is closer to the 
RT m easurements. Measure ments of the sheet conduc- 
tivity l|lleinrich et all |2£)03b') indicate that the Au layers 
used in the experiments have non-negligible residual re- 
sistances. [We note, however, that the conductivity en- 
tering Eq. (|82|l does not include the interfacial-scattering 
contribution; the measurement of the sheet conductiv- 
ity therefore does not give us a direct information about 
CTros-] Assuming for example a^es ~ Cph would yield the 
K (with a = a^es) and RT (with 0-1 = a'l + a~^) 
lines in the close vicinity of the measured points. 

The theoretical results represented by the straight lines 
in Fig. llOl are based upon the asymptotic, single-interface 
value of gl^ for Au|Fe from Table ^ This approxi- 
mation needs to be relaxed in order to study possible 
size-dependent corrections in thin films. To estimate 
the variation which ca n result from finite-size effects, 
IZwierzvcki et all l2005l carried out a series of calcula- 
tions for a Au|Fe|vacuum syste m, using vacuum ins tead 
of GaAs as in the experiment l)Urban et all l200lt) for 
simplicity. The mixing conductance of the other, Fe|Au, 
interface in Eq. (|S3 was kept at its asymptotic value 
(Table nj. The calculated thickness ((i)-dependent mix- 
ing conductance gl^ was then converted into the Gilbert 
damping via Eq. (|83|l . The results for perfect (specular) 
structures, marked in Fig. ^| with black crosses (x) , ex- 
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FIG. 10 Enhancement of the Gilbert damping coefficient for 
an Fe|Au|Fe trilayer as a function of 1/d where d is the 
thickness of the excited Fe layer in monolayers (ML's). The 
filled circles ( •) are the room-temperature values measured by 
lUrban et all 1200 ll a nd the open oii e (o) is a low-temperature 
(77 K) value from iHeinrich eTail l2003bl The theoretical 
prediction based on Eq. H83|l for K (with a oo) is 
shown as solid and the RT-corrected (with phonon scatter- 
ing) one as dashed lines. The results of K calculations for 
a Au|Fe|vacuum system are given by crosses (x) and stars (*) 
for specular and disordered interfaces, respectively. The value 
of the Gilbert d amping a for a single F e film is marked with 
an arrow. From IZwierzvcki et aUr2005l 



hibit oscillations of non-negligible amplitude about the 
asymptotic values given by the solid line (arbitrarily tak- 
ing the low-temperature regime, i.e., tr ^ oo for refer- 
ence). The introduction of interface disorder (two ML's 
of 50%-50% alloy) yields values for the damping [stars (*) 
in Fig.^] essentially averaged back to the limit given by 
the single-interface calculations of Table ^ 

iLubitz et all l2003l also reported an increased mag- 
netization damping in polycrystalline Cu|Fe|Cu|Fe|Cu 
multilayers as compared to Cu|Fe|Cu structures. The 
additional damping scaled as 1/d with magnetic-film 
thickness d, but was considerably larg er than that 
for ep itaxially-grown systems reported by lUrban et al\ . 
1200 iL Besides, this damping was rapidly reduced by in- 
creasing the thickness of the Cu spa cer to only several 
nanometers, which was interpreted bv^L ubitz et a/J . l2003l 
to be due to possibly a very short spin-flip length in their 
Cu. We note that, based on Eq. H82() . this could be alter- 
natively explained by a short elastic scattering length in 
the polycrystalline Cu. .Lubitz et ai . 2003 found a mod- 
erate increase by about 10% in the additional damping 
on lowering the temperature to 77 K, which is roughly 
consistent with the decrease in the phonon contribution 
to the resistivity. 

This subsection demonstrated that direct first- 
principles calculations can produce values of the damping 
coefficient in the same range as those measured experi- 
mentally in good-quality structures. Moreover, by taking 



into account various other sources of scattering in the Au 
spacer and/or quantum-size effects, the calculations can 
be brought into close agreement with observations. A 
more conclusive comparison with experiments would re- 
quire a detailed knowledge of the microscopic structure 
of the experimental system, which is currently not avail- 
able. 



V. DYNAMIC EXCHANGE INTERACTION 

A. Magnetic bilayers 

The ground-state energy (free energy at finite temper- 
atures) of more than one magnetic layer embedded into 
a nonmagnetic medium depends on the relative orienta- 
tion of the magnetic moments. This is the essence of 
the static exchange coupling discussed in Sec. II. Bl Dis- 
order scrambles ballistic electron paths connecting mag- 
netic layers and exponentially suppresses the static in- 
teraction as a function of the nonmagnetic spacer thick- 
ness between magnetic layers and the inverse mean free 
path. We can picture the dynamic coupling in terms 
of a localized magnetic moment that suddenly changes 
its direction, thus creating a nonequilibrium local spin 
accumulation. The latter partly precesses in the local- 
moment effective field and partly diffuses away. Other 
magnetic moments at distances L experience this spin 
accumulation after a time td — L^/D, where D is the dif- 
fusion coefficient, as long as td < Tgt, the spin-flip time, 
or, equivalently, L < ^/DTsf, the spin-diffusion length. 
These other moments start to move by the torque they 
experience by absorbing part of the diffusing spin accu- 
mulation, and, in turn, emit spin currents by themselves. 
When the magnetic dynamics are slow on the scale of the 
diffusion time td, the retardation of spin diffusion may be 
disregarded and the dynamic exchange coupling is prac- 
tically instantaneous. This is the regime we treat in the 
present section. (A similar discussion applies to the bal- 
listic transport regime.) The dynamic coupling between 
moving magnetizations by the exchange of nonequilib- 
rium spin currents is affected by spin-conserving random 
scattering much less than the rapid suppression of the 
static coupling. 

When the normal spacer is much thinner than the spin- 
diffusion length, the dynamic spin exchange is governed 
by ^p|n|f' ^1- H'^2|l . which does not depend on the N- 
spacer width L for ballistic spacers. The coupling de- 
cays as 1/L for spacer widths larger than the scattering 
mean free path and is exponentially suppressed when L 
becomes larger than the spin-diffusion length. The dy- 
namic exchange interaction is less sensitive to disorder 
than its static counterpart because the latter relies on 
the orbital quantum interference between electron tra- 
jectories connecting magnetic moments/layers, whereas 
the former requires spin coherence on traversing the nor- 
mal spacer. The dynamic exchange coupling has been 
first addressed perturbatively in the context of electron 
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spin resonance bv IB ar nesl 1 1 9 74l who already pointed out 
its long-range nature as compared to the static coupling. 
However, the dynamic correlations between an individual 
magnetic moment and proximate conduction electrons 
(which, as described in Sec. lIVI lead to magnetization re- 
laxation) were overlooked at that time. In the context of 
ferromagnetic-resonance experiments, dynamic exchange 
coup ling has been studied by Hurdcquint and Malouc he, 
19911 and, more recently, by iHeinrich et alL l2003al "and 



Lenz et a/.l.l2004 



The spin-pumping picture of dynamic exchange cou- 
pling is impressively confirmed by rece nt experiments 
with s ufficiently large normal spacers l)Heinrich et al\ . 

as will be discussed in detail below. In the 
limit when the static coupling becomes appreciable, the 
spin-pumping-based circuit analysis leading to Eq. H82|) 
may not hold exactly, however, since the circuit model 
of Sec, in assumes that the interfaces and normal spacers 
scramble the electron distribution in momentum space. 
A sizable static exchange is thus in principle inconsistent 
with the key assumption underlying the circuit model. A 
rigorous treatment of the dynamic coupling in the regime 
of significant static exchange coupling is difficult, as ex- 
plained in Sec. IVII.A.2I Here we rather heuristically as- 
sume that the dynamic coupling as described by spin- 
pumping theory still holds in the presence of a resid- 
ual static exchange interaction, which turns out to be 
sufficie nt to make a conn ection with some recent experi- 
ments l|Lenz et al 1 12004 - This picture can be justified by 
the general observation that in intermetallic heterostruc- 
tures quantum size effects often modulate the semiclas- 
sical transport contribution only weakly. 

In the following, we consider the collective dynamics 
of an F|N|F structure, i.e., two magnetic films separated 
by a normal metal. A possible static exchange through 
very thin spacers is taken into account phenomenologi- 
cally by postulating a Heisenberg-type contribution (per 
cross-sectional unit area) 



Ex = — Jmi • m2 



(84) 



to the magnetic free-energy functional that determines 
effective fields ^ experienced by the ferromagnets. J is 
the Heisenberg coupling constant which is assumed to be 
small compared to the magnetic bulk exchange stiffness 
A divided by the magnetic-film thickness d, \J\ <C A/d. 
This assumption is necessary in order to treat each in- 
dividual magnetic layer as a macrospin pointing along a 
unit vector mi. J depends in an oscillatory fashion on the 
spacer-layer thickness and favors either parallel (J > 0) 
or antiparallel (J < 0) orientation of the magnetic layers, 
as discussed in Sec. II.BI In nanostructured pillars, the 
magnetostatic interaction can also favor an antiparallel 
coupling of the form (|84|l , even at spacer layer thicknesses 
at which the exchange coupling vanishes. 

We consider magnetic films that are thinner than v4/| J| 
but thicker than Asc [Eq- H13|) ]. so that they completely 
absorb transverse spin currents. (Note that in typical 
metallic structures, A/|J| ^ Xp and Asc ~ Af-) A 



precessing magnetization vector mi of ferromagnet 
pumps spin angular momentum at the rate H38|) deter- 
mined by the spin-pumping efficiency »4p|j^|p [Eq. ijH^Jl] 
into the normal-metal spacer. We concentrate here on 
small- angle excitations of a collinear magnetic equilib- 
rium configuration. Eq. I|82|l applies to the typical situ- 
ation that the magnetization dynamics are slow on the 
characteristic time scales for electron transfer across the 
spacer. When one ferromagnet is stationary, see the left 
drawing in Fig. 1111 the dynamics of the other film are 
governed by the LLG equation ^ but with the effec- 
tive damping parameter enhanced with respect to the 
intrinsic value, as given by Eq. (|83|l . When both mag- 
netizations are allowed to move, see the right sketch in 
Fig. ^] the coupled LLG equations expanded to take 
into account the spin torques 114|) read 



dt 



= - ■jiin, X ( Hcff,j 
dmi 

a) I m. X 



Jm, 



dm, 



aiUii X 



dt 



dm.j 
— m , X — 

dt dt 



(85) 



where Heff ^ are effective fields not including the exchange 
contribution a\ = ?i7iRey4,p|i^|p/(47rMs^idiS') is 

given by Eq. H83|l . having disregarded Imylp|j^|p, and 
j = 1(2) for i = 2(1). As a first simple example, con- 
sider the parallel equilibrium configuration, m^"'' — mj^"*, 
with zero static coupling, J = 0, and matched reso- 
nance conditions: a\ = a2 — a. and 7iHoff,i = 72Hcff^2, 
with magnetization-independent effective fields Hcft,i. 
After linearizing Eq. I|85|l in terms of small deviations 

Ui(t) = mi(t) — m^^-* of the magnetization direction m^ 

from its equilibrium value raf^ , we immediately see that 
u = (uisi + U2S2)/(si -I- Si), where si = jiMg^idi, viz., 
the symmetric mode, is damped with the intrinsic Gilbert 
parameter a, whereas the difference Au = Ui — U2, viz., 
the antisymmetric mode, relaxes with enhanced damp- 
ing constant a — a + a'l + a'l- This demonstrates that 
the dynamic interaction can lead to nontrivial collective 
magnetization dynamics even when the static interaction 
vanishes. 

Let us now analyze a more general case of a coupled 
magnetic bilayer undergoing a collective circular preces- 
sion near a parallel equilibrium configuration m^*^-* = 

10.^2 '' — For simplicity, we still focus on a nearly sym- 
metric structure by setting 7^ = 7, jJ/{Ms.idi) = ojx, 
Ui = a, and a'^ = a', but allow the effective fields 
7Hoff,i = ojS to differ, uji ^102- Eq. (|85|) then reduces 
to 



dm, 



dm,: 



dt dt 



a m.i X 



drxii 
'dT 



dm, 
m,- X — 

^ dt 



(86) 



The linearized equations of motion in the absence of 
driving force are solved by the form oc exp{iLut) with 
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FIG. 11 A cartoon of the dynamic-coupling phenomenon. 
In the left drawing, layer Fi is at resonance and its precess- 
ing magnetic moment pumps spin current into the spacer, 
while F2 is detuned from its FMR. In the right drawing, both 
films resonate at the same external field, inducing spin cur- 
rents in opposite directions. The short arrows in N indicate 
the instantaneous direction of the spin angular momentum 
oc mi X liii carried away by the spin currents. Darker ar- 
eas in Fi around the interfaces represent the narrow regions 
in which the t ransverse spin momentum is absorbed. From 
iHeinrich et aLLl2003a. . 



two complex-valued natural frequencies uj and definite 
circular polarization. Solid lines in the main panel of 
Fig. 1121 show the real part of these w's for various ratios 
ujijiox after setting lo^ = 0, and dashed lines after setting 
LOx/^i = 0.01, i.e., introducing a finite static exchange 
coupling. In both cases, a = a' = 0.02. The lower in- 
set shows the corresponding normalized imaginary part 
of (jj. For LOx — Q and resonance frequencies uji and W2 
well separated on the scale of the enhanced damping a', 
viz., |w2/wi — 1| ^ 2a' (assuming a' ^ 1), the dynam- 
ics of two ferromagnets decouple and the spin pumping 
can be accounted for by simply adding a' to the effective 
Gilbert parameter of each F layer. When, on the other 
hand, |w2/wi — 1| ^ 2a' ^ the spin pumping locks the 
collective dynamics to independent symmetric (acoustic) 
and antisymmetric (optic) normal modes with frequen- 
cies that are nearly degenerate and close to [oji +uj2)/2. 
It then follows from Eq. (|86|1 that the symmetric mode is 
weakly damped, with Im(a;)/Re(aj) close to a, whereas 
the antisymmetric mode experiences an enhanced Gilbert 
damping a -\- 2a' . 

The presence of a static exchange interaction lifts the 
frequency of the antisymmetric mode by 2ujx^ not affect- 
ing the symmetric mode. When the frequencies ui are 
well separated, the static interaction lifts both of them 
by ujx ■ One of the modes then acquires some symmetric 
and the other some antisymmetric character even when 
\ijj2/oJi — l]^^ 2a', as illustrated by the damping of both 
modes (dashed lines in the lower inset of Fig. I12|l . In 
the upper inset of Fig. ^1 the calculations are repeated 
for finite uox, but setting the intrinsic damping of one of 
the films to zero, ai = 0, while a2 = a'^ = 0.02. The dy- 
namic locking into symmetric/antisymmetric modes close 
to the resonance crossing has disappeared. The modes 
are harder to synchronize when the participating modes 
have different amplitudes. A significant locking may be 
expected only in a symmetric bilayer with the individual 
films having similar resonant modes near parallel equi- 
librium axes. In particular, a bilayer in an antiparal- 
lel equilibrium configuration is not disposed to dynamic 
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FIG. 12 The main panel displays the real part of the two 
frequencies uo that solve Eq. 18611 for the linearized dynamics 
of a symmetric magnetic bilayer close to the parallel mag- 
netization configuration. Solid lines are obtained by setting 
the static coupling to zero, lOx = 0, and dashed lines for 
uOx = COlciJi. The lower inset shows the corresponding imag- 
inary part of UJ. When ajx = and near the frequency cross- 
ing, uj\ ~ U2, the symmetric mode has a smaller imaginary 
part (less damping) and only slightly larger real part. For 
uJx < 0, the symmetric mode has a larger (real part of the) 
frequency than the antisymmetric mode, whereas the anti- 
symmetric mode acquires a higher frequency for sufficiently 
large and positive lox- In all cases, the antisymmetric mode 
remains to be stronger damped. The upper inset shows the 
results for lux ~ O.Oloji but setting ai = while keeping a 
finite Q2. This inset illustrates the effect of a bilayer asym- 
metry. 



locking, since individual layers have excitations with op- 
posite circular (or elliptic, in the presence of anisotropics) 
polarizations. 

We now turn to a discussion of the consequences for 
the observable in FMR experiments, viz., the energy- 
dissipation power P of the dynamically-locked collective 
dynamics as plotted in Fig. E| as a function of the fre- 
quency to of the uniform and circularly-polarized trans- 
verse rf driving field. The system parameters are the 
same as those used to generate the dashed lines in Fig. 1121 
For the smallest ratio ^2/^1 in Fig. ^1 the two bi- 
layer excitations are recognized as separated Lorentzian 
peaks with halfwidths close to (a -I- a')uj, see Eq. (|II|I . 
When the elementary bilayer excitations become locked 
at 1 072/(^1 — 1| ^ 2a', the rf radiation excites only the 
symmetric mode, and the two Lorentzians merge into 
a single sharper Lorentzian with halfwidth of only acu. 
This narrowing is explained by the cancellation of spin 
currents and can be quite dramatic when a <^ a' . 

In FMR experiments, the applied magnetic field is 
swept, whereas the rf frequency uj is fixed by a resonant 
cavity. In Sec. | IV.G| we discussed FMR experiments by 
lUrban et al\ . 12001 on magnetic bilayers in which one of 
the layers remains close to equilibrium while the other 
is resonantly excited. The collective magnetization dy- 
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FIG. 13 Dissipation power P of a symmetric bilayer in the 
presence of a uniform circularly-polarized rf driving field at 
frequency uo. The parameters are the same as those used to 
generate the dashed lines in Fig. 1121 



namics were measured by iHeinrich et al\ . l2003al in the 
same system by making use of an accidental crossing of 
the resonance frequencies when the static field is reori- 
ented relative to the crystal anisotropy axes. The result- 
ing spectra could be quantitatively explained in terms 
of the dynamic coupling in the limit of a vanishing static 
exchange interaction. We summarize their findings in the 
following. 

The MBE-grown structure of IHeinrich et a4l2003al in- 
corporates two ferromagnetic films, a thinner 16 ML (Fi) 
and a thicker 40 ML (F2) Fe film, separated by 40 ML's 
of Au, grown on GaAs and capped with Au, i.e., the 
stacking order is GaAs|16Fe|40Au|40Fe|20Au(001). The 
uniaxial magnetic anisotropy in Fi at the GaAsjFe inter- 
face can be used to intentionally tune the resonance fields 
for Fi and F2 into a crossover by rotating the static mag- 
netic field by an angle with respect to the (001) crystal 
axis. In a finite interval of near the crossover (shaded 
area m Fig. Hljl, the two FMR fields clearly "stick" to 
each other, a phenomenon explained above. When the 
resonance fields are identical. Hi = Hi^ the rf magne- 
tization components of Fi and F2 are moving in phase 
as depicted in the right drawing in Fig. 1111 For similar 
trajectories of Fi and F2, the total spin current through 
the spacer vanishes resulting in zero excess damping for 
both films, as follows from Eq. H85I) . The locked collec- 
tive motion is then hindered only by the intrinsic local 
damping. This is experimentally verifi ed, as shown in 
Fig- El For a theoretical analysis. Heinrich et a^J . [200331 
solved Eq. , taking into account the ellipticity of the 
magnetic motion caused by the anisotropies. Using pa- 
rameters derived from measurements on the uncoupled 
layers, they calculated the total FMR signal as a function 
of the difference between the resonance fields, Hi — Hi, 
without additional fitting parameters. The predictions 
are compared with the measurements in Fig. 1151 The re- 
markably good agreement between experiment and the- 



ory provides strong evidence that the dynamic exchange 
coupling not only contributes to the damping but leads 
to a new collective beh avior of magnetic heterostructures. 
IHeinrich et all l2003al additionally carried out measure- 
ments on samples with Au spacer thicknesses between 
14 and 100 monolayers: The reported weak dependence 
of the FMR response on the spacer thickness proves the 
long range of the dynamic interaction. 
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FIG. 14 Dependence of the static FMR fields Hi (circles) 
and H2 (triangles) of the thin Fe film (Fi) and the thick Fe 
film (F2), respectively, on the angle ip of the static magnetic 
field with respect to the Fe [100] crystallographic axis. The 
sketch of the in-plane measurement in the left inset shows 
how the rf magnetic field (double-pointed arrow) drives the 
magnetization (on a scale grossly exaggerated for easy view- 
ing). The right inset shows the measured absorption peaks 
for layers Fi and F2 at = 60 Deg. The absorption power 
is given by the imaginary part of the susceptibility of the rf 
magnetization component along t he rf driving field, whic h is 
denoted in the figure by x" ■ From IHeinrich et ttZ.ll2003al . 

iLenz et'cUl . l2004l investigated collective FMR dynam- 
ics in Ni|Cu|Ni and Ni|Cu|Co structures with Cu thick- 
nesses down to 2 ML's. Such thin Cu spacers sup- 
port a sizable static coup l ing between mag netic films. 
Like IHeinrich et all l2003al iLenz et ai\ . l2004l observed a 
sharp drop in the linewidth near an FMR-field crossing. 
Far from it, the optic-mode resonance is systematically 
broader than that of the acoustic mode, consistent with 
the spin-pumping mechanism. The difference between 
the optic and acoustic linewidths exhibits an oscillatory 
dependence on the Cu-spacer thickness, which roughly 
follows the predicted exchange-coupling constant J. As 
shown in the lower inset of Fig. 1121 a stronger coupling 
J would indeed result in a larger asymmetry in the two 
linewidths, away from the crossing, which could lead to 
oscillations in the linewidth difference. However, an in- 
homogeneous spread of the coupling strengths may also 
contribute to the broadening which depends nonmono- 
tonically on the spacer thickness, providing an alterna- 
tive explanation for the linewidth oscillations. 
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FIG. 15 Comparison of theory (solid lines) with room-temperature measurements (symbols) close to and at the crossover of 
the FMR fields, marked by the shaded area in Fig. 1141 The left and right frames show FMR signals for the field difference, 
H2 — Hi, of -78 Oe and -1-161 Oe, res pectively. The theor etical results are parametrized by the full set of magnetic parameters 
which were measured independently llUrban et all l200lfl . The magnitude of the spin- pump current was determined by the 
linewidth at large separation of the FMR peaks and agrees well with that predicted theoretically, see Sec. liV.CI The middle 
frame displays the effective FMR linewidth of magnetic layers for the signals fitted by two Lorentzians as a function of the 
external field. At Hi = H2, the FMR linewidths reached their minimum values at the level of intrinsic Gilbert damping of 
isolated films. The calculations in the middle frame did not take small variations of the intrinsic damping with angle (p into 
account, which resulted in deviations between theory and experiment for larger \Hi — H2\- Notice that AHi first increases 
before attaining its minimum, which is due to a contribution of the antisymmetric collective mode. As a side comment, it 
should be noted that although fitting the absorption signal by two Lorentzians is a legitimate approach to comparing theoretical 
calculation with experimental curves, the analysis does not imply that the signal is alw ays well approximated by the sum of 
two Lorentzians, which may not be the case very close to the FMR-field crossover. From lHeinrich et a/.. .2003a . 



In the regime of very strong static exchange coupling, 
I J| > A/d, opposite to what was assumed so far, magneti- 
zation gradients across the bulk magnetic layers alleviate 
the energy cost of the discontinuity in the magnetic ori- 
entation between two magnetic layers. This reduces the 
spin pumping through the spacer, and thus the additional 
broadening of the antisymmetric mode, because the adja- 
cent magnetizations are better locked by the strong static 
exchange interaction. 



B. Magnetic superlattices 

Magnetic multilayers display a rich pattern of physi- 
cal properties that ha,ve be en well investigated (see, e.g., 
ICamlev and Starnol . Il993l for a review). However, the 
relevance of the dynamic exchange coupling discussed in 
the previous subsection on the spin-wave dispersions and 
lifetimes appears to have not been recognized yet. In 
the following, we present a simple model description that 
should suffice to estimate the order of magnitude of the 
predicted effects that hopefully will stimulate new exper- 
iments. We also remark that superlattices can serve a 
theoretical purpose as toy models for describing certain 
features of bulk magnetism. 

Consider a periodic stack (in the x direction) of alter- 
nating F and N layers forming a two-component superlat- 
tice. We treat the model depicted in Fig. 1161 in which an 
F|N bilayer forms the unit cell with thickness b = L + d, 
where the normal-metal spacer of width L separates the 



magnetic films of thickness d ^ Age. Translational in- 
variance is assumed in the lateral directions. We consider 
here collective spin-wave excitations, taking both static 
and dynamic exchange couplings into account, in exactly 
the same fashion as in the previous subsection on mag- 
netic bilayers, modeling each magnetic layer as a single 
macrospin. 

b=L+d 



t 



t 



FIG. 16 A schematic view of the superlattice and its geomet- 
ric parameters as considered in the text: An F|N bilayer is 
repeated along the x axis, with either ferromagnetic or anti- 
ferromagnetic alignment of the consecutive magnetic layers. 
The system is translationally invariant along the two remain- 
ing axes. 

The small-angle magnetization dynamics of a multi- 
layer in an all-parallel configuration are described in 
terms of local deviations from the equilibrium: Ui{t) — 
m.i{t) — m(°^. For long- wavelength excitations, it may 
be approximated as a continuous function u(a::, t) of the 
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coordinate x normal to the interfaces. For a uniaxial 
effective field i?off — -ffcffm''*'\ the spin- wave dynamics 
obey the differential equation 

dtu =m'-°-' X [wqU — LUxb^dxxM 

+adtu-a'b^dxx^tn] , (87) 

where we used the quantities defined for magnetic 
bilayers in Sec. IV. Al iv^ — iJ/iMgd), a' — 
/i7ReAj^|j^|p/(47rMsdS'), and ujq = jHch. The second 
term on the right-hand side of Eq. H87|) is due to the static 
exchange interaction mediated by quantum- well states of 
the spacer layers, and in the last term we recognize the 
dynamic coupling induced by the spin pumping. The sec- 
ond spatial derivatives reflect simply the difference of the 
spin currents through two consecutive normal spacers in 
the continuum limit. The static Heisenberg coupling can 
be interpreted as the superlattice equivalent of the bulk 
exchange-stiffness parameter A, which for the superlat- 
tice is given by ^ = Jb^ /d. Both oj^ and a' depend on the 
normal- interlayer thickness L. It follows from Eq. H87|) 
that the small- momentum, k <C b~^, spin- wave excita- 
tions of the superlattice, propagating perpendicular to 
the interfaces, oc exp{i[kx — u!{k)t\}, obey a dispersion 
relation 



Loik) 



Up + {bk^uJx 
l + i[a+ {bkfa'] 



When fc — > 0, uj(k) reduces to the Larmor frequency wq 
of the individual magnetic layers because the static and 
dynamic exchange couplings vanish when the consecutive 
magnetic layers move coherently in phase, as explained 
in Sec. IV. Al Eq. (|88|l holds for momenta comparable to 
b~^ when bk is replaced by 2sin(&fc/2). 

The situation is very different for an 
antiferromagnetically-aligned superlattice, which is 
the lowest-energy state when, for example, J < and 
Hf.({ = 0. In this case. 



uj{k) = 



2a + {bkfa' ± ^Aa^ - (6fc)2(l 



1 



Aaa' -f a'{bky 



(89) 

where plus and minus signs refer, respectively, to the 
modes with antisymmetric and symmetric dynamics in 
adjacent layers for overdamped motion, and to the right- 
and left-propagating modes when the real part of uj{k) 
is significant. Note that now uj^ < 0, so that Imw > 0, 
as required for a stable configuration. In the absence of 
bulk magnetization damping, a = 0, Eq. (|89|l reduces to 



uj(k) 



±{bk)ujx 
1 ± i{bk)a' 



(90) 



with linear dispersion and damping at small k. Eqs. H89|) 
and H9U|I can also be generalized to large momenta by 
replacing bk with 2sin(&fc/2). Notice that in Eqs. (|87|l . 



(|88|l , and (|90|l , the dynamic coupling modifies the damp- 
ing similarly to the way the static coupling affects the 
excitation frequency of the magnetic superlattice. Crys- 
tal and shape anisotropics on top of the simple effective 
fields assumed above, might become important in real 
structures, and their inclusion is straightforward. 

FMR experiments access the multilayer dynamics from 
the sample surface down to the microwave skin depth 
Askin ~ 100 nm (which is even smaller in BLS). Therefore 
only modes with momenta k > tt/ Askin can be measured. 
Since the skin depth decreases when the temperature is 
lowered, both the FMR frequency and the damping in a 
ferromagnetically-aligned multilayer should grow roughly 
as l/Agjjjjj cx r, the momentum scattering time in the 
normal skin-effect regime. Such studies can thus give 
information about the temperature-dependent scattering 
in superlattices; a weak temperature dependence could 
indicate that scattering is dominated by structural dis- 
order. Inelastic neutron-scattering spectroscopy may be 
useful in elucidating the collective dynamics in thick mul- 
tilayers, especially if supported by elastic neutron scat- 
tering (Fitzsimmons et al., 2004,) to probe the magnetic 
profile in the superlattice. 



C. Large-angle motion in biased spin valves 

Perpendicular spin valves, i.e., Fs|N|F;i trilayer pillar 
structures with layer thicknesses down to a few mono- 
layers and submicron lateral dimensions, are ideal sys- 
tems to study precession and switching phenomena in 
magnetic heterostructures. By attaching contacts on the 
outer sides an electric bias can be applied perpendicu- 
lar to the interface planes. In many experiments, F^ is 
a "soft" ferromagnetic film with a magnetization that 
can change easily, whereas ¥h is a "hard" magnetic layer 
whose magnetization is assumed to be stationary. The 
relevant variable is then the time-dependent magnetiza- 
tion of the soft layer. The soft layer can be excited by a 
current bias or an applied rf magnetic field (or both) . In 
the second case (realized as, e.g., an isolated magnetic bi- 
layer of Secs. lIV.'O and rV'^ . ¥h can be pin ned by an ex- 
chang e bias or surface magnetic anisotropy l|Urban et al\ . 
l200l|f . In the former case, the magnetization of layer Yh 
may be also rendered less sensitive to a given spin torque 
simply by growin g it much thicker tha n F^ or by a resis- 
tance anisotropy l)Kovalev et a/.l 120021) . For a sufficiently 
thick spacer N, the static interaction between the ferro- 
magnetic layers can be disregarded, while the dynamic 
coupling induced by the spin pumping may be still siz- 
a ble, see Sec.IV.Al 

ISlonczewskil Il996l and iBereeil Il996l predicted that 
spin valves should display novel time-dependent effects. 
ISlonczewskl 119961 realized that a current flowing through 
a spin valve causes a spin transfer through the non- 
magnetic spacer, i nduci ng spin torques on the ferro- 
magnets. iBergeil Il996l predicted that the two ferro- 
magnets should interact even without an applied elec- 
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trie current, resulting in a significant contribution to 
the Gilbert damping of the magnetization dynamics. He 
further showed that an electric current can excite zero- 
momentum spin- waves in the ferromagnet, a n idea that 
was l ater supported experimentally, see, e.g.. lTsoi et al\ . 
I2nnnl The condit ion for the resonant spin-wave emis- 
sion ijBergeil Il996tl is similar to the criterion for magne- 
tization reversal by Slonczcwski, 1996, 1999., who treated 
the Gilbert damping par ameter as a phenomenological 
constant. iBerged 120011 however found a dependence 
of the damping parameter in spin valves on the rela- 
tive magnetization angle. Some of Berger's and Slon- 
czewski's results as well as the underlying theoretical 
models are thus not consistent with each other. The spin- 
pumpi ng concept unified the seminal work of these pio- 
neers l|Tserkovnvak et oil . l2003b) as explained in the re- 
mainder of this subsection. We will calculate the critical 
current bias for the low-temperature magnetization in- 
stability and the configuration-dependent Gilbert damp- 
ing parameter. The treatment here is limited to the 
macrospin model, i.e., to small systems in which the mag- 
netic layers are monodomain ferromagnets characterized 
by two magnetization vectors. If com bined with micro- 
magnetic simulations ijLee et a/.U2004^ , the full range of 
the precession and switching dynamics can be studied in 
principle, however. 

Consider the system sketched in Fig.[T7| The Fs|N|F,i 
trilayer is sandwiched between two normal-metal con- 
tacts sustaining a charge-current bias Ic- The soft layer 
Fg magnetization mi then starts moving from its equi- 
librium direction at a critical value that depends on the 
applied magnetic field. Thermal activati on facilitates 
current-induced magnetization switching ijMvers et all 
linn^), but we focus here on the low-temperature regime. 
The spin torque on the magnetization of F^ in the pres- 
ence of a spin current Isi flowing from F^ into the normal 
spacer is given by Eq. H14(l . The spin current 
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(91) 



consists of the dynamic-exchange current 1°™^ induced 
by the spin pumping (|38|l and of the current driven 
by an applied current bias. The former is responsible 
for a dynamic coupling between the ferromagnets, see 
Sec. IV.AI and, as we discuss in the following, can be in- 
terpreted as a viscous friction term that stabilizes the 
relative magnetization configuration of the spin valve 
against the torques exerted by I^^^ or an applied mag- 
netic field. In high-density metallic systems, the applied 
voltages and spin accumulations are safely smaller than 
the Fermi energies, which means that we are in the linear- 
response regime and both spin currents may be calculated 
independently of each other. Spin pumping in the out- 
ward direction, i.e., into the external connectors, would 
only increase the intrinsic damping coefficient by a con- 
stant value, as discussed in Secs. lIV.Xl and HV.BI and can 
thus be added trivially to the intrinsic damping a. Spin 
pumping into the normal spacer, which gives rise to I^f'^'^, 
requires more attention. 




FIG. 17 Schematic of a current-biased spin valve (a "dy- 
namic" version of Fig. 0. The symbols are explained in the 
text. 



We start by considering the spin current H38|l pumped 
into the spacer by a time-dependent mi (t) in the absence 
of an applied current bias, Ic = 0. The following assump- 
tions are convenient and realistic: (i) The magnetic films 
are sufficiently thicker than Age, so that the transmis- 
sion contribution to the spin-pumping parameter A^'^ can 
be disregarded, (ii) the mixing conductance g^^ is real- 
valued, (iii) two magnetic films including the interfaces 
have the same conductances, (iv) the normal spacer is 
ballistic, and (v) the spin-flip processes are restricted to 
the ferromagnetic layers, with the spin-diffusion length 
Asd ^ Asc. (Note that in this subsection Agd denotes the 
spin-diffusion length inside the ferromagnets for longitu- 
dinal spin transport along the magnetization direction). 
The spin currents out of the magnetic layers into the 



normal spacer are then given by I^. 
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(93) 



where the last equation is derived using Eqs. 
and a circuit-theory analysis discussed in Sec. III.AI as- 
suming vanishing charge current. Here, m^ is the ith- 
layer magnetization direction and A/x^j = /x^jy ~ ^J■sF^^^^ 
is the spin- accumulation difference across the NjF^ in- 
terface. The magnetization-precession period is typically 
much longer than the electron dwell time in metallic spac- 
ers. Assuming weak spin-flip scattering in N, conserva- 
tion of angular momentum then implies that IsiH-Is2 — 0. 
It is now straightforward to calculate the exchange spin 
current which is given by I^i'^'^ = Isi when Ic = 0. 

The longitudinal component of the spin accumulation 
/Xg^ can penetrate into the ferromagnets on the scale 
of the spin-diffusion length Agd, whereas the transverse 
component vanishes on the shorter scale of Age near the 
interface. The longitudinal spin accumulation in the fer- 
romagnet, ^sFi, and thereby I^^'', can be obtained for a 
given fi^ff from the diffusion equation for the (longitudi- 
nal) spin transport in the ferromagnets, similarly to the 
normal-layer spin diffusion discussed in Sec. IIV.BI To 



34 



be specific, we take both the charge and spin currents 
to vanish on th e outer boundaries of F., a nd Ff^. It can 
then be shown ijTserkovnvak et all . I2fl08b^ that the lon- 
gitudinal spin current flowing into a ferromagnetic slab 
of thickness d is governed by an effective conductance g* 
defined by 



5* 



gsd tanh((i/Asd) 



(94) 



where g^d = (Ve^)(S'/Asd)(2crV^)/(o-T + a^) and a" is 
the spin-s conductivity of the ferromagnetic bulk, so that 



r(o) 



-gl^im X fi^j^ X m,] 



(95) 

The transverse spin current is determined simply by gj^ , 
since we have taken Age to be the shortest relevant length 
scale in the problem. Note that g* ^ when d <C Agd, 
i.e., when the spin-flip relaxation vanishes, or when the 
ferromagnet is halfmetallic, so that it completely blocks 
the longitudinal spin flow for a vanishing charge flow. A 
new parameter 



(96) 



characterizes the asymmetry of the absorption of trans- 
verse vs longitudinal spin currents. Putting together 
Eqs. 1)92(1 and ((95|l and demanding conservation of an- 
gular momentum in the spacer (i.e., I^i -I- 1^2 = 0), one 
arrives (after some algebra) at 
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m2 — vm.1 cosf 



1 — j/^ cos"^ 9 

(97) 

where cosO = mi • m2. Since the normal spacer was 
taken to be ballistic, /x^,^ is uniform. The exchange 
spin current would otherwise be somewhat suppressed 
by the spacer diffuse scattering, which can be taken 
into account easily by solving the spin-diffusion equa- 
tion as in Sec. 



IV.BI if necessary. Let us estimate typi- 
cal values of v for sputtered Cu|Co and Cu|Py systems 
at low temperatures, taking d = b nm. The main dif- 
ference between the two material combinations is the 
spin-diffusion length in the ferromagnets: Co has a rel- 
atively long Asri ~ 60 nm, whi l e Ash ~ 5 nm is very 
short in Pv llBass and Pratt. Jri Il999t iFert and Pirauxl 
^99; Pi raux et a/.l . Il996f) . Using known values f or spin- 
dependent conductivi t ies llBass and Pratt. Jr.l Il999t 
iFert and Pirauxl Il999l: iPiraux et all Il996j) . we obtain 
<7sd/<5 ~ 2.7 nm~^ for Co and 16 nm~^ for Py. 
2^TT^ii I j^gll) I s « 20 nm-2 for the Cu|Co interface. 
Tabled we may expect the value for Cu|Py to be similar. 
With ql^ I S ^i 27 nm-^ f or CuIC o, Table|I| and 15 nm'^ 
for Cu|Py l|Bauer et aZ.Ll2003b^ . one flnds u « 0.98 for 
Cu|Co and v « 0.33 for Cu|Py. 

The magnetization dynamics in the absence of an ap- 



plied bias are determined by substituting I^'^'^ into equa- 



tion H14|l . which thus becomes inconsistent with a con- 
stant effective Gilbert parameter. We now analyze the 
configuration dependence of the damping in more detail, 
which can be measured, in principle, by the FMR line 
width broadening at high rf intensities (and therefore fi- 
nite "precession cones"). For mi precessing around m2. 



mi X I'f^ X mi = - 1 , , 
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The angular dependence of the additional Gilbert damp- 
ing p arameter due to the excha nge spin current then 
reads ()Tserkovnvak et all l2003h|) 



am 
a'(0) 



= 1 - 
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(99) 



where a'(0) = jhgl^ / (SnMsdS) is the damping enhance- 
ment in a coUinear configuration, see Eq. H83|l . For small 
angles, 6* w 0, Eq. (|?^ can be rewritten as 
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where 
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(101) 



Eq. 1(1 00(1 was also obtained bv lBergeil l200ll in the limit 
of small precession angles, however, for a definition of s 
with s (X Tst differing from Eq. 1(101(1 . Expressions 1(99(1 
and ifTUni) are plotted in Fig. E 




FIG. 18 Solid lines display the precession-cone angle depen- 
dence of the Gilbert damping parameter, Eq. 19911 . and the 
dotted lines show the extrapolated small-angle approxima- 
tion 1(100(1 . The lowest lines are representative for Co and the 
uppermost ones for Py, assuming magnetic-film thickness of 
5 nm. Fe and Ni are expected to be intermediate between 
these two cases. 

As mentioned above, i/ is close to 0.98 for cobalt, so 
that the lower solid line in Fig. ^| represents the damp- 
ing for Co according to Eq. 1(99(1 . s = 333 is found by 



35 



iBergeil I2OO1I for Cu|Co with Co 1.5 nm thick, which is 
remarkably similar to the estimate based on v for this 
thickness when substituted into Eq. 1)101(1 . The damping 
is thus predicted to be significantly reduced for precession 
angles deviating only slightly from the coUinear configu- 
rations (we expect this conclusion to be also true for Fe 
and Ni). Modeling the magnetization dynamics with a 
constant damping parameter is therefore not allowed for 
sufficiently thin magnetic layers. For permalloy, on the 
other hand, the magnetization damping is expected to 
remain significant at all angles, see the upper solid line 
in Fig. 1181 This implies that the magnetization reversal 
by an applied magnetic field should be faster in Py than 
in Co spin valves. 

If mi rotates around an axis perpendicular to xn.2, i.e., 
only the relative angle 6 changes, then l?"™? _L and 
Eq. (123 reduces to I^f = IP™p/2. The corresponding 
damping a'(0) thus has an angle-independent enhance- 
ment with respect to the intrinsic Gilbert damping. 

An applied current bias is an additional instrument to 
control the magnetization dynamics. When the conduc- 
tance parameters of the spin valve are mirror symmetric, 
the bias-induced spin transfer I^}'^'' is coplanar with the 
magnetization directions and can be written as 
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for a given applied bias. This follows from expanding the 
spin current as 

jbias ^ /^;^(cos6')mi -|-/22(cos6l)m2 -f /i2(cos6')mi x m2 

(103) 

and noting that I^i'*^(mi, m2) = Ig{'*'^(m2, mi), by sym- 
metry, which implies that /n = /22 and /12 — 0. The 
electric current necessary to induce a given spin current 
jbias (depends on 9 and can be calculated readily by the 
circuit theory summarized in Sec. IT01 Eqs. IHJ), (EI), 
(|97|l and the form of the bias current completely deter- 
mine the dynamics of mi{t). The exchange induced by 
the spin pumping causes relaxation toward an equilib- 
rium configuration, while the bias current can either re- 
lax or excite the magnetization, depending on the sign of 
the current, as discussed below. 

Near the parallel configuration, 6* w 0, Eq. H97|l simpli- 
fies to Isf^*^ = I51 "^^/2. Let mi circularly precess around 
a fixed m2 with the frequency w (determined by effective 
fields): mi x rhi = wmi x m2 x mi. The total projected 
spin current has the Gilbert form: 
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Instability is reached when the effective Gilbert-damping 
coefficient becomes negative. The critical bias is thus 
given by 



s.crit 



An 



(105) 



where a is the intrinsic Gilbert constant. Neglecting the 
first term in t he brackets, one obtain s a result analo- 
gous to that of lSlonczewskii Il99(l Il999t while neglecting 
the second term leads to a conditio n similar to t he res- 
onant spin- wave emission criterion (^Bergeil Il996|) . The 
spin-pumping contribution (the first term) is compara- 
ble with the intrinsic damping (the second term) for 
transition-metal films with thickness d of several nanome- 
ters, but dominates for very thin films. When the in- 
stability is reached, the trajectory of mi(i) can become 
very complicated and could possibly lead to a magne- 
tization reversal to a different (meta)stable configura- 
tion. A more complete discussion of spin torques and 
macrospin switch ing dynamics in asymmetric spin valves 
can b e found in iBrataas et all l2005at iManschot et all 
12004 Nonlinear large-angle dynamics and field-induced 
switching in the presence of spin pumping were dis- 
cussed bv lKim and Chappertl. 12005* within the macrospin 
model. Nonuniform dynamics have been studied by 
iBrataas et all l2005bl Experimental evidence for a spin- 
pumping-induced da mping in current -driven spin valves 
has been reported bv iKrivorotov et' al. . 2005. The im- 
portance of the spin-pumping contribution to the critical 
current (I105II h as been investigated experimentally by 
ISun et oill2005l for Co|Cu|Co spin valves. 



VI. LINEAR-RESPONSE APPROACH 

A. Heterostmctures 

This subsection is devoted to an alternative descrip- 
tion of the spin emission by a dynamic ferromag- 
netic magnetization embedded in a conducting nonmag- 
netic matrix, as in F ig. |21 which was put forward by 
SimMigk_andJIeinricl, , 12003 and further elaborated by 
MillsL l2003t ISimanekL l2003l l2004 In Sec. |lll| we formu- 
lated such spin emission as a scattering pumping process, 
which requires the concept of waveguide leads for electron 
states that are reflected by or transmitted through a fer- 
romagnetic layer between normal reservoirs. The time 
dependence of the scattering matrix caused by a moving 
magnetization leads to the pumping of spin currents, and 
the corresponding loss of angular momentum by the fer- 
romagnet to an increased viscous damping of the magne- 
tization dynamics. Such language is stan dard in t he field 
of mesoscopic transport phenomena (^Beenakkeil [l997) , 
but the magnetism community is more accustomed to 
linear-response susceptibilities rather than scattering ma- 
trices. Of course, the final result should not depend on 
the formalism used, as verified in the following, but one 
may argue that the scattering theory has distinct advan- 
tages over the linear-response approach for the present 
problem. 

For a linear-response formulation, we proceed from a 
Hamiltonian for conducting electrons experiencing a fer- 
romagnetic exchange field as given by Eqs. H49|l . H50|l . 
When the exchange field Vx = —Kl/2 is uniform in- 
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side the ferromagnetic volume V and vanishes outside, 
we may rewrite Eq. (|50|) as 



[ [s(r) • M(r, i)] 
Jv 



(106) 



Here, s(r) is the spin-density operator for conduction 
electrons that are polarized by the exchange field with 
strength Q along the magnetization direction m = 
M/Afs. M is a collective property that is treated as 
a classical time-dependent potential. The Hamiltonian 
(|106|l follows, e.g., from a mean-field approximation of 
the s — d model. The exchange is then viewed as an ex- 
ternal potential (provided by the d electrons) that does 
not depend on the s-electron distribution. The spin- 
density-functional formulation of the mag netizatio n dy - 
namics in itinerant ferromagne ts, see, e.g.. iCaDelle etall . 
1200 It Iq ian and Vignalei I200I leads to the Hamiltonian 
Eq. H106|l in the local-density approximation. The effec- 
tive field lO due to the interaction (|106|l then reads 



H;f,(r,t) = ^(s(r)), 



(107) 



In the s — d model, this field corresponds to the reac- 
tion torque by the nonequilibrium conduction-electron- 
spin distribution inside the ferromagnet, which is excited 
b y the moving mag netization direction. 

ISimanek and Heinrich. 2003 suggested to calculate the 
reaction torque on the ferromagnetic magnetization di- 
rectly by using Eq. (|107|l in the LLG equation of motion. 
This appears to be very different from the spin-pumping 
picture that relies on the spin currents that are emit- 
ted at the interface to the normal metal. The two ap- 
proaches are however related by the continuity equation 
for electron spin dynamics, see, e.g.. lCapelle et aU . l200ll 
according to which the spin current equals the torque 
exerted on the magnetization M by the effective field 
(|107|l . up to a term given by the time derivative of the 
average spin density, (s)t, times volume. When (s)t fol- 
lows the magnetization adiabatically, the difference be- 
tween the spin current and torque is thus proportional to 
m. In the case of the s — d model, this has a physical 
meaning, corresponding to the torque required to change 
the angular momentum of the s-electron spin density. 
For Stoner ferromagnets, in which the conduction elec- 
trons are identical to the ones that make up the mag- 
netization, we have to calculate the spin flow emitted 
by the ferromagnet as an additional term in the LLG 
equation, while the reaction torque on the magnetiza- 
tion by the conduction electrons has no obvious physical 
meaning. When treating itinerant (Stoner) ferromagnets 
within the density-functional theory, the spin currents in- 
duced by the Kohn-Sham Hamiltonian (|106|l in principle 
differ from the physical ones inside or very close to the 
ferromagnet. However, this difference vanishes asymp- 
totically as a functio n of distance from the ferromagnet 
llCapelle et oi.Ll200lD . The emitted spin currents are thus 
similar to those in the s — d model, and their evaluation 



can thus be mapped on calculating corresponding reac- 
tion torques. 

A uniform small-angle dynamics (of a possibly large 
51) only weakly perturbs the system. The induced spin 
imbalance {Ss{r))^ can therefore be expressed in terms 
of the linear-response susceptibility of the unperturbed 
system 
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where &{t) is the Heaviside step function. The effective 
field due to the induced nonequilibrium spin density is 
then given by 
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where a, a' € {x, y, z} are the indices of the Cartesian 
axes and a stands for the corresponding unit vectors. 
Suppose now for simplicity that the system is invariant 
under spin rotations about the z axis, and consider small- 
angle magnetization dynamics near this axis, dm m — 
z. Substituting Eq. (|109|l into the LLG equation © then 
yields the following lowest-order dynamic term: 

^ , , , 751^ / . dm , dm\ 
-7mxmi,(t).l^^A,-+A,mx-j , (110) 

which is the most general adiabatic term f or axially- 
symmetric systems (|Mills and R.ezendS . l2003|) . Here, 
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are real numbers. Ai renormalizes the effective gyromag- 
netic ratio 7cff in Eq. © and A2 is a Gilbert-like damping 
parameter 



Q!eff 



dio 
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Eq. (|112|l can also be obtained by equating the en- 
ergy dissipated into the itinerant degrees of freedom 
by moving magnetization and the work done by an 
rf magnetic field applied against the effective viscous 
Gilbert term, at a steady magnetic precession, since 
liniij^o Rexs^s (^)/'-^ = in the case of spin-rotational 
symmetry around the z axis. For thermally-equilibrated 
s-electron subsystem, Eq. (I112|) implies aoff7off > 0, as 
required by the LLG phenomenology, see Sec. II. CI 

The damping constant a can generally be expressed in 
terms of the response function of the total magnetization 
(e.g., of the s and d electrons in the s — d model) by in- 
verting Eq. (|ll|l . In the present discussion, this problem 
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is self-consistently reduced to the simpler task of calculat- 
ing the spin response of s electrons to a time-dependent 
exchange field aligned with the d-electron magnetization. 

The similarity between the torque (|110|l and the spin- 
pumping current H38() should not be s urprising i n the 
light of the above discussion. Indeed, ISimaned l2003l 
explicitly demonstrated that the low-frequency linear- 
response and spin-pumping pictures lead to identical a^g 
for a (5-function magnetic layer embedded into a clean 
normal metal (which corresponds to an ideal spin sink 
for emitted spin currents in the spin-pumping language) . 
However, evaluating the linear-response correlation func- 
tions becomes very tedious for more realistic models, see, 
e.g., [Mills, 2003'. It is also not obvious how to treat the 
current-induced spin transfer and the spin pumping on 
the same footing by the linear-response formalism. On 
the other hand, the formulation in terms of the suscepti- 
bilities is complementary to the spin-pumping approach. 
It appeals to the intuition of many researchers in the 
magnetism community. Furthermore, it could be help- 
ful to obtain insights into problems that are hard to 
solve within the scattering theory. Examples of these are 
magnetic bilayers coupled by a strong static exchange in- 
teraction (see Sec. IVII.A.2II . strongly-correlated systems 
fSec. lVll.U|l . and the bulk damping of the magnetization 
dynamics fSec. IVI.B|) . 



B. Bulk clamping 

In the previous subsection, we considered coupling of 
the magnetization to itinerant electrons via an exchange 
interaction (|l()6|l . In particular, it was reasserted that 
when a ferromagnetic film is inserted into a nonmagnetic 
metal, the magnetization dynamics causes an emission of 
spins. In the presence of a spin sink outside the ferromag- 
net, this pumping leads to a net angular- momentum loss 
that is equivalent to an excess damping of the magneti- 
zation dynamics, as discussed in Sec. IIVI An analogous 
mechanism could be effective in bulk ferromagnets pro- 
vided that spin-relaxation decay channels exist that dis- 
sipate a spin accumulation created by the pumping. This 
happens in the presence of momentum-scattering mech- 
anisms, such as lattice impurities and phonons, and the 
accompanying (or band-structure) relativistic spin-orbit 
interaction. The inductive coupling of the magnetization 
to the conduction electrons that causes dissipation by 
eddy currents is less important in thin layered structures 
and is disregarded in the following. For simplicity, the 
discussion is restricted in the following to a mean-field 
s — d model. 

A mechanism of bulk magnetization damping bearing 
similarity with the spin-pumping picture in heterostruc- 
tures was proposed a long time ago by iMitchellL Il957t 
It involves a transfer of the angular momentum (and 
energy) of a nonequilibrium ferromagnetic configuration 
to the itinerant electrons via the exchange interaction, 
with a subsequent spin-orbit relaxation to the lattice. 



The consequences of such process for macrospin (long- 
wavelengt h) dynamics have been worked out for the s — d 
model by iHeinrich et all Il967l The s — d pictu re has 
been resurrected recentlv bv ISinova et all l2004l in or- 
der to address the magnetization relaxation in the fer- 
romagnetic semiconductor (Ga,Mn)As, in which ferro- 
magnetism originates from the free-hole-mediated ex- 
change interactio n between di lute, substitutional spin- 
5/2 Mn atoms l)Ohnd. Il999(l . It is possible to re- 
produce and gene r alize the results of IHeinrich et all 
119671: ISinova et all l2004l by calculating the magnetiza- 
tion damping in terms of the conduction-electron-spin 
dynamics i n a time-dependent excha nge field similarly to 
Sec. NTH fiTserkovnvak et all 120041) . 

Consider an sp—d model of a conducting ferromagnet, 
where the spin S of the itinerant s or p states (either 
electrons or holes) experiences an exchange field of mag- 
nitude Q along the uniform magnetization direction m of 
the localized d orbitals, as in Eq. (|106|l : 



H{t) ^Ho- nm{t) ■ S . 



(113) 



Here, Hq is a one-particle Hamiltonian refiecting the host 
band structure. In (Ga,Mn)As, substitutional Mn are 
paramagnetic acceptors that strongly interact with the 
free holes. Although the exchange field can be highly 
nonuniform on atomic scales, it is customary to start 
with a simplifying assumption that it may be smeared 
out. As before, we treat the magnetization m as a clas- 
sical and, on the relevant length scales of the carrier dy- 
namics, spatially-uniform variable. A suitable model for 
the valence bands of dilute p-doped semiconductor (e.g., 
GaAs, Si, or Ge) is the spherical Luttinger Hamiltonian 
for spin-3/2 holes: 



Hn = 



1 



2mp 



71 + ^1-2 ]P 



272 



pS 



, (114) 



where nie is the free-electr on mass and t,; are the so- 
called Luttinger parameters l|Luttingeii Il956|) . The spin- 
orbit term couples the hole momentum p with its spin 
S. The four-band model (|114|l . is valid when the hole 
Fermi energy is sufficiently smaller than the spin-orbit 
splitting of the semiconductor host between the spin- 
3/2 and spin-1/2 valence bands. Eq. (|114|l neglects 
corrections for lattices with broken inversion symmetry. 
Spin-1/2 electron systems with vanishing spin-orbit cou- 
pling are recovered by setting 72 = 0. Suppose the 
magnetization of the localized orbitals varies slowly in 
time (being uniform at all times), so that the time- 
dependent m modulates the Hamiltonian (|113|l adiabat- 
ically. This means that the system equilibrates on time 
scales faster than the motion of m and all the quan- 
tities parametrizing the carrier Hamiltonian stay con- 
stant. The observation of ferromagnetic resonance indi- 
cates that time-dependent long-range ferromagnetic or- 
der ind eed exists in thin film s of m agnetic transition 
metals l|Heinrich and CqchraE , Il993l) and semic onduc- 
tors (jGoennenwein et all 12005 ; iRappoport et a^J . l2004i|1 . 
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Eq. (|112|l may then be taken as a microscopic definition 
of the Gilbert-damping parameter. In the following, a is 
formally evaluated for electron and hole systems with an 
emphasis on its dependence on spin dephasing. (We will 
drop the subscripts "efF on a and 7 from the previous 
subsection.) 

We initially focus on a system without spin-orbit cou- 
pling in the band structure, 72 = 0. Suppose the trans- 
verse spin-density dynamics in the exchange field is 
described by the Bloch equation 



ds 
dt 



fls X m{t) 



som(i) 



(115) 



where the last term is a phenomenological relaxation due 
to impurities, parametrized by the transverse dephasing 
time T2. Let us assume that m(i) moves slowly on the 
scales of fi~^ and T2. It is then convenient to transform 
Eq. (|115() into a frame of reference (for spin variables) 
that moves together with m.{t). When, for example, m, 
at a given instant, rotates with frequency uj around the 
y axis, see inset of Fig. it is stationary in the rotat- 
ing frame at the cost of a new (Larmor) term ws x y on 
the right-hand side of Eq. H115|) , which corresponds to a 
magnetic field along the y axis. Since the motion is slow, 
it is sufficient to solve for s as the instantaneous station- 
ary state in the moving frame of reference. If m points 
along the z axis in the rotating frame, spin polarization 
along the x axis then exerts a damping torque on m that 
corresponds to 



7Xs 



nr. 



7*0 



M, 



1 + (f]T2)2 Ms ' 



(116) 



where Xs^sy is the stationary (real- valued) response func- 
tion in the rotating frame. (We recall that 7 is the gy- 
romagnetic ratio of d orbitals.) The calculation of the 
time-dependent linear response in the laboratory frame, 
Eq. I|112|) , is thus simplified to that of the static response 
in the rotating frame. Such a transformation is how- 
ever not possible in general, in particular for Hamilto- 
nians that are not spin-rotationally invariant, such as 
Eq. l(TB|) . Eq. l(TTC|) is plotted in Fig. [H The equihb- 
rium spin density Sq can be calculated from the specific 
form of the static Hamiltonian. a vanishes at both small 
and large spin-relaxation limits. 

The damping parameter 1)116(1 depends only on the 
ratio of spin-relaxation rate and exchange energy. The 
low spin-relaxation rate regime, a (x T^j"^; analogous 
to the spin-pumping damping of thin films in contact 
with a spin-sink conductor: The moving magnetization 
"pumps" spins into the itinerant carriers at a constant 
rate, which arc then relaxed with a probability cx 
before exchanging spins with the ferromagnet again. The 
difference is that now the spins are pumped into the 
internal conduction electrons of the ferromagnet rather 
than those of an externally attached metal. The other 
limit, a cx T2, can be understood by noting that for 
a very viscous dynamics of s{t), s(t) « So(t) — T2So(t) 




FIG. 19 Gilbert damping, Eq. (11161 . in units of Omax ~ 
7So/(2Ms) as a function of the normalized spin-relaxation 
rate. Inset: Geometry of the model. 



in the laboratory frame, a cx T2 then follows from 
the torque cx (s — Sq) x m. This is analogous to the 
"b reathing Fermi-surface" dampi ng mechanism discussed 
bv lKunes and Kambersk^ l2002l in the regime of fast re- 
laxation: The itinerant carriers try to lower their energy 
by rearranging themselves in the field of slowly vary- 
ing magnetization direction but lag behind with a short 
delay time determined by the relevant relaxation pro- 
cesses. In the presence of an anisotropic spin-orbit in- 
teraction in the fcrromagnct's crystal field, the breathing 
Fermi surface gives an additional contribution to damp- 
ing ijKunes and Kambc rskvl |206^ . see also a short dis- 



cussion in Sec. I VII. A. 21 In the present model described 
by equation H115|) . the "breathing" takes place in spin 
space. 

It is interesting to note that Eq. II116I) reduces to 



It IS mterestmg to note tnat n^q . Iliibll reauces t^ 
the result obtained bv lHeinrich et a/lll967l for the long 

wavelength magnon lifetime due to the s — d interaction 
with spin-1/2 conduction electrons in the random-phase 
approximation: 



7 [nm*kF/{4Tr'^h)] 



1 + (r!T2)2 



Ms 



(117) 



where m* is the band-structure mass and kp the Fermi 
wave vector. Here, it was assumed that hft ^ Ep (the 
Fermi energy). The quantity in the square brackets is 
jus t the total carrier spin density. Eq. (|117|l was used 
by llngvarsson et al\ . l2002l in order to explain the mea- 
sured damping in thin permalloy films, which scaled lin- 
early with the film resistivity p. A linear relation be- 
tween the two quantities is expected when fl 3> T^^, 
so that a cx T^^, and using Drude formula p cx r"^ (t 
is the transport mean free time) and assuming T/T2 = 
const which depends on material and scattering-impurity 
type but not on scattering rate ijAbrikosov and Gor'kovl 
Unlike for ferromagnetic semiconductors, how- 
ever, the use of the s — d model for itinerant ferromagnets 
like transition metals is questionable, since the separate 
treatment of magnetic and conduction electrons is un- 
physical. 

Let us turn to the magnetization relaxation in hole- 
doped magnetic semiconductor (Ga,Mn)As. Eq. H116|) 
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may be used to obtain a rough estimate of the damp- 
ing coefficient: The largest achievable value of amax — 
jso/{2Ms) occurs when the holes are fully polarized, giv- 
ing ttmax ~ 0.1 — 0.3, roughly one third the ratio of the 
(spin-3/2) hole to the substitutional (spin-5/2) Mn con- 
centrations. For realistic samples with a spin polariza- 
tion of the order of unity, therefore, a,nax ^ 0.1. The 
damping a is further suppressed by the factor a/amax = 
2riT2[l-|-(17T2)2]-i < 1. For clean bulk samples of GaAs 
the spin-relaxation time is ~ 100 fs ijHilton and Tan j . 
12003). For approx imately 5% Mn doping, Hfl ^ 0.1 eV 
l^ietl aZl . l200lj) . so that ^ 10. This corresponds 
to the a oc T,^^ regime with a ^ 0.01. Reduced spin- 
relaxation times should thus result in a larger damping. 
Experimentally, the impurity scattering is likely to be the 
easiest parameter to vary in order to tune a to a desired 
value. The strong spin-orbit coupling 72, however, makes 
the validity of the phcnomcnological equation l|115l) , and 
thus the result (|116|l . questionable for the hole system. 
Besides, a strong crystal anisotropy would require a fur- 
ther refinement of the analysis. Returning to our basic 
equation H112|l and inserting the response function for a 
noninteracting system yields 

a lim - J2 mS.\j)f [/FD(e.) - /FDfe)] 

Ms V i^^O LO ^ — ' 
13 

X S{huj + Si - Ej) , (118) 

where i, j label the one-particle eigenstates of the sam- 
ple with volume V. When the wave vector k is con- 
served, Y.^j /V = E<T<T' / d^k/{2TT)^, where a, a' label 
spin states. For a perfect crystal, therefore, a vanishes, 
unless there is a finite Fermi-surface area with spin de- 
generate states. Introducing lattice defects leads to a 
finite a, see, e.g.. ISinova et gj ], l2004l (these authors did 
not include important vertex corrections in the response 
function, however). 

In the above analysis, we have assumed a coherent mo- 
tion of the ferromagnetic magnetization without specify- 
ing the source of excitation. An FMR magnetic field 
(with a large dc and small rf components), for example, 
can be included explicitly into the Hamiltonian H113|) of 
the itinerant carriers. The results for the Gilbert damp- 
ing will stay unaffected, however, as long as the exchange 
energy Ml is much larger than the carrier Zeeman split- 
ting in the applied field and the ferromagnetic magneti- 
zation is mainly carried by the localized orbitals. (The 
energy pumped into the carrier-magnetization dynamics 
by the rf field must be taken into account otherwise.) 
Inhomogeneities in the bulk magnetization are not im- 
portant on the length scales set by the transverse spin- 
relaxation rate and the precession frequency in the ex- 
change field. The bulk spin dynamics discussed in this 
subsection have no effect on the spin pumping into ad- 
jacent conductors discussed in Sec. IIIII as long as the 
transverse spin-relaxation rate is small compared to the 
exchange precession frequency. In the opposite, rather 
unrealistic, limit, the s-electron spin dynamics are locked 



to the d-electron magnetization motion, suppressing spin 
leakage (pumping) into adjacent normal conductors. 



VII. MISCELLANEOUS 
A. Quantum-size effects 

1. Ultrathin magnetic layer 

As explained in Sec. IIIII the spin pumping by a mag- 
netic layer in contact with normal metals is governed 
by the parameter A^^ = g^^ - t'"^^, Eq. When 
the magnetic-film thickness d exceeds the spin-coherence 
length (|13|l . d ^ Age, t'^^ vanishes and A^^ is given sim- 
ply by the interfacial spin-mixing conductance g^^. In 
this subsection, on the other hand, we focus on the regime 
in which d is smaller or of the order of Age, i.e., for thick- 
nesses of a few monolayers in the case of transition-metal 
ferromagnets. In this limit, the coherence between up- 
and down-spin states in the ferromagnet leads to a thick- 
ness dependence of A^^, and thus of the spin pumping 
and magnetization torque (|27|l . The decoherence of the 
orbital wave function due to inelastic scattering processes 
is disregarded, assuming its characteristic length scale is 
much longer than Age. 

As in most of this review, we discuss in the follow- 
ing layered structures. A large lateral area S renders 
mesoscopic phenomena such as the Coulomb blockade ir- 
relevant, and we focus on quantum-coherence effects due 
to small layer thicknesses. For a study of magnetiza- 
tion dynamics in magnetic nanoparticles in contact with 
large reservoirs, see, e.g.. IWaintal and Brouweil l2003t 
IWaintal and Parcolletll2fm 

The linea r -respo nse framework 

ijSimanek and Heinrichl l2003() . see also Sec. IVII has 
been used to calculate the enhanced Gilbert damping 
of finite-thickness ferromagnetic films embedded into a 
conducting medium by Mills, 2003. For an idealistic 
model of isotropic band structure with spin-dependent 
exchange potential, he found that ultrathin films display 
oscillatory damping (as a function of thickness) due 
to quantum-size effects. However, by first-principles 
scattering matrices based on realistic electronic struc- 
tures computed in the l ocal spin-density approximation, 
IZwierzvcki et aO . l2005l have shown that quantum-size 
oscillations are much smaller than those obtained by 
iMild. I2OO3I and small amounts of dis order suppress the 
remaining oscillations even further. IZwierzvcki et aO . 
.200 5 also found that the spin-pumping torque is of the 
Gilbert-damping form, with only a very small correction 
to the gyromagnetic ratio. 

The above observations have been made for cop- 
per/cobalt and gold/iron heterostructures, of which 
we will restrict our discussion to the former. Real- 
istic band-structure and simple-model calculations on 
Cu|Co|Cu(lll) trilayers are compared in Fig. 1201 where 
the real part of Gi, — {e^ /h)t^^ /S is plotted as a func- 
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tion of thickness d of the magnetic layer (measured in 
atomic monolayers). Solid circl es in fipjure | 20 | show the 
result of the ab initio calculation llZwierzvcki^ al\ . \2m!t\ 
without impurities and for specular (k|| -conserving) in- 
terfaces. The smooth solid lines represent the calcula- 
tions for the isotropic free-electron model. For exchange 
splittings A = 2, 4, 6 eV, the amplitude of oscillation is 
much larger and the decay is much slower in the model 
than in the more realistic first-principles calculations. As 
might be expected, increasing the exchange splitting from 
2 to 6 eV leads to a shorter period and more rapid decay 
of the oscillations. In order to mimic the parameter-free 
ab initio result, an exchange splitting in the range of 
10 eV would be needed (not shown in the figure), how- 
ever. Such a large value can be justified neither on the- 
oretical nor experimental grounds. This discrepancy il- 
lustrates the difhculty of quantitatively representing the 
complex electronic structure of transition metals by sim- 
ple model Hamiltonians. 
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--- A=4eV 
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FIG. 20 The real part of G]^ = {e^/h)t^^/S as a function of 
the magnetic-layer thickness d calculated from first principles 
(•) for realistic multiband electronic structures and that cal- 
culated for an isotropic free-electron gas with Fermi energy 
ep = 7 eV (chosen to obtain the correct value of the Sharvin 
conductance i n Cu) and various choic es of the exchange split- 
ting A. From lZwierzvcki et alll2005l 

Figure 12] shows = (e^ /h)g^^ / S and G*^^ calcu- 
lated from first principles in the presence of disorder mod- 
eled by a monolayer of 50% alloy added on each side of 
the magnetic layer. The thickness d in this case is defined 
as that of the remaining clean ferromagnetic layer. The 
disorder strongly quenches the amplitudes of the oscilla- 
tions as a function of d, so that G|| is practically constant 
at the level of its asymptotic (i.e., single-interface) value. 
For G||, the oscillations do not vanish completely, but 
their amplitude is substantially reduced to values that 
are negligibly small compared to ReG^j^ for all but the 
thinnest magnetic layers. Diffusive scattering in the bulk 
of the magnetic layer, which has not been considered here 
explicitly, should have similar effects. 

Both G|| and G*^ are governed by (i) the matching 
of the normal-metal and ferromagnetic-metal states de- 
scribed by the scattering matrix of the isolated inter- 
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FIG. 21 Spin-mixing conductances of a Cu|Co|Cu(lll) struc- 
ture with disord ered interfaces as a fun ction of the Co-layer 
thickness. From lZwierzvcki et oi.L 1200,'tI . 



face and (ii) the phases accumulated by the electrons on 
their quantum-coherent propagation through the mag- 
netic layer. The interface determines the amplitudes of 
the oscillations and, the asymptotic value of G^^, whereas 
the bulk term is responsible for the oscillation period. It 
is instructive to interpret the transmission and reflection 
coefficients of the finite-size magnetic layer in terms of 
multiple scattering between the two interfaces. The task 
is simplified by the simple Fermi surface of Cu, which 
corresponds to only one left- and right-going state at 
the Fermi energy for each value of k|| and spin. The 
sums over states in Eqs. H23|l and 1)24(1 therefore reduce 
to integrations over the two-dimensional Brillouin zone 
(2DBZ) involving the complex-valued functions r''^(k||) 
and t°'(k||). To lowest order in the number of scattering 
processes and dropping explicit reference to kn. 



t° 



where t 
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= [ti, . . . jf^y is a vector of transmission 
coefficients between a single propagating state in the nor- 
mal metal and a set of states in the ferromagnet, A'^ is a 
diagonal matrix of phase factors e^^^^"^ {j is an index of 
the states in the ferromagnet), r'^^ff is a scalar reflec- 
tion coefhcient for states incoming from the normal metal 
and Tp^p is a square matrix describing reflection at the 
ferromagnetic side. The set of states in the ferromagnet 
consists of both propagating and evanescent states. The 
contribution of the latter decreases exponentially with 
the thickness of the layer. 

Let us first analyze the thickness dependence of t^K 
Substituting Eq. H119() , the summation in Eq. (|24() is car- 
ried out over terms containing phase factors e*^'^*^"'^^^''''. 
Because of the large differences between majority and mi- 
nority Fermi surfaces of the ferromagnet, this typically 
leads to rapidly oscillating terms that to a large extent 
cancel in the 2DBZ integration over ky. In th e spirit 
of the theory of interlayer exchange coupling l)Brunol 
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ll99Rt IStiles and Zangwilli 120021) . long-ranffe contribu- 
tions must originate from the vicinity of points for which 
^k|| — k'jj^) ~ 0, corresponding to a constant phase 
in the integrand of Eq. (|24|l . The damped oscillations 
around zero are therefore caused by the stationary-phase 
Fermi-surface calipers. 

By substituting Eq. ^(T^ into g^^, Eq. two 
thickness-independent contributions can be identified. 
The first, coming from integrating the Snn' term in 
Eq. H23|l . is just the number of transport channels 
in the normal metal, i.e., the dimensionless Sharvin 
conductance. The second contribution comes from 
the T^N^^T^N^N term and provides an interface-specific 
correction to the first contribution. The thickness- 
dependent contributions contain, to lowest order, phase 
factors e'-'-^'i^'^^'^^'^'^ and e~*'^''^J-+''^'j-)''. Just as in the case 
of t^-'-, their 2DBZ integrals have an oscillatory character, 
with periods governed by different Fermi-surface calipers. 
These oscillations now occur around the constant value 
set by the first two contributions. The asymptotic value 
of g^^ is clearly the reflection mixing conductance of a 
single F|N interface. 

In metallic films, the electronic structure of all but the 
outermost atomic layers is practically identical to that of 
the bulk material. The period of oscillations of g^^ and 
as a function of the magnetic-layer thickness d is thus 
a bulk property of the magnetic layer. The amplitudes, 
on the other hand, involve the interfacial scattering co- 
efficients introduced in Eqs. (|119|l and H12U|) . When the 
quantum-size oscillations are small, ~ g^-'^, where 
g'^'^ is a property of the N|F interface instead of the en- 
tire structure. Furthermore, g^^ can be estimated by the 
Sharvin conductance of the normal metal, Eq. (|64l) . The 
results of single-interface calculations are listed in Ta- 
ble ^ for clean and disordered interfaces of Cu|Co and 
Au|Fe material combinations. The disorder in Table ^ 
was modeled by 2 monolayers of 50% alloy instead of a 
single monolayer in the present subsection. In spite of 
this difference, the values are practically identical to the 
asymptotic ones seen in Fig. [2] for the Cu|Co combina- 
tion. 



2. Ultrathin normal spacer 

Let us now turn to a discussion of the magnetization 
dynamics of two monodomain magnetic layers separated 
by an ultrathin normal-metal spacers with quantum-well 
states that penetrate into and couple the ferromagnets. 
As explained in Sec. lI.BI the free energy F of the system 
depends on the relative angle between the two magne- 
tizations even in the absence of magnetic anisotropics. 
The dependence of F(Mi), i — 1,2, on the magnetic 
configuration corresponds to nonlocal effective fields ® 
exerting torques on the magnetizations. In the following, 
we comment on the dynamic component of the exchange 
interaction in time-dependent problems, cf. the semiclas- 
sical dynamic exchange interaction discussed in section 



m 

Consider for simplicity an s — d model with nonin- 
teracting s electrons, where the magnetic d orbitals are 
coupled to itinerant electrons via a mean-field exchange 
interaction. The transverse component of effective field 
O entering the LLG equation ||SJ| for the ith magnetic 
moment is then given by 

H.W=-^^a™, (i?(m.))t 

= - ^^y dm,^en{m.i)n^{m.r,t), (121) 

where H{nii) is the Hamiltonian for itinerant electrons, 
parametrized by the magnetization directions m,i, and ()t 
denotes the (quantum-mechanical) expectation value at 
time t. The sum on the second line of Eq. H121|l runs over 
all eigenstates of H{m.i), where k labels both the orbital 
and spin degrees of freedom, is the occupation number 
corresponding to the many-body state at time t and is 
the energy of the Kth eigenstate. Setting the many-body 
ensemble at time t to its equilibrium value determined by 
mi{t) reproduces the Landau-Lifshitz (LL) definition ||HJ), 
leading to dissipationless trajectories [assuming there arc 
no other sources of damping, i.e., a = in Eq. ||SJ|]. Such 
an approximation thus captures only the static exchange 
component. 

In reality, ()t lags behind its instantaneous equilibrium 
value. The corresponding nonequilibrium component 
of (mi , t) in momentum spa ce reflects the "breath- 
ing F ermi surface" discussed bv lKunes and Kamberskvi 
120021 in the context of transition-metal bulk magneti- 
zation damping in th e presence of crystal anisotropy. 
iHeinrich et a^J . l2003bl conjectured that such mechanism 
may also play a role in the magnetic dynamics of F|N|F 
structures: The modulation of the exchange energy 
stored in the normal spacer may cause an additional 
damping through the time lag in the itinerant-electron 
response. 

For a spin-rotationally-invariant system, the effective 
field (|121(l reduces to 

Hr^(i)--^^E-^-^-."«("^^'^)- (122) 

The bulk damping coefhcient pi6|) determined by the 
transverse spin-relaxation time T2 in a model discussed 
in Sec. IVI.BI where the average electron spin density 
lagged behind that corresponding to the instantaneous 
equilibrium, can be understood to arise from such an 
effective field. Spin-orbit interaction in the presence 
of a crystal field in bulk ferromagnets and/or the ex- 
change coupling through the normal spacer modulate 
^^(mi) via the time-dependent rrii, leading to an ad- 
ditional dynamic contribution to the effective field, as 
follows from Eq. (|121|l . This contribution has a partic- 
ul arly simple form in the limit of short relaxation times 
T ()Heinrich et all l2003bt iKorenman and Prangel Il972t 
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iKunes and Kamberskvll2002|) : 

MsiVi ^ drai \dva..i dt J 

(123) 

where £f is the Fermi energy and assuming low tem- 
perature. It is clear that this effective field results in a 
(tensor) Gilbert damping. 

In disordered structures and/or thick spacer layers 
with vanishing static exchange coupling between the 
magnetic films, Eq. H121(l should in principle reproduce 
the dynamic exchange coupling discussed in Sec. In 
other words, the spin pumping captures the semiclassi- 
cal component of the time-dependent exchange coupling 
between ferromagnetic films when the static contribution 
vanishes. 



is 

L,a^J2^aa'iM)^, (124) 
aa' 

with constraints on the form of the 3x3 tensor Aaa' (M), 
a, a' = 1,2,3, for a given crystal symmetry. In particu- 
lar, Is^a may now have a component along the magne- 
tization direction. It is always the case, for example, 
when the ferromagnetic exchange field varies in magni- 
tude as well as in direction. By conservation of the to- 
tal angular momentum, the spin-pumping current H124|l 
is accompanied by a torque on the magnetization, and 
also transferred into the orbital angular momentum as 
well as lattice torque. Note that in the presence of mag- 
netic anisotropies, the LLG equation of motion is a ten- 
sor equation in which the scalar Gilbert parameter a is 
replaced by a 3 x 3 tensor l)Mills and Rezend3 . l2003|) . 



B. Spin-orbit coupling 

The derivation of the spin-pumping current in Sec. IIIII 
relies on Eqs. H21|l and H22|l which relate the scattering 
matrix in spin space, Snn',u', for given channel and lead 
indices, n, n' and /, I', to the magnetization direction 
m. For systems isotropic in spin space, the scattering 
matrix depends on m only through the simple projection 
(|22|l . In transition- metal ferromagnets, this is a good 
approximation, since their exchange splitting is by far the 
largest relevant energy scale. A large spin-orbit coupling 
in the electronic structure, such as in p-doped magnetic 
semiconductors like (Ga,Mn)As, on the other hand, can 
considerably modify the spin-pumping currents. 

In ferromagnets with spin-orbit interaction, the 
rotating-frame analysis of Sec. 1111.51 becomes tedious by 
the need to apply the transformation to the orbital as well 
as spin degrees of freedom. Crystal anisotropies or even 
the presence of planar interfaces would make such an ap- 
proach impractical. The adiabatic-pumping formalism of 
Sec, nil. A1 still applies, however. One may in general cal- 
culate the tensor current (|45|l in terms of the emissivity 
(|46|l . when the dependence of the full scattering matrix 
on the magnetization direction is known. Of course, the 
2x2 matrices have to be generalized to (25+1) x (25+1) 
for spin-5 carriers. In the case of spin-orbit coupling in 
the nonmagnetic leads, however, such tensor currents are 
in general different from spin currents, since spin is not a 
good quantum number for transverse quantum channels. 
This complicates the matter beyond the scope of this ar- 
ticle. In the simplest case of weak spin-orbit coupling 
within the ferromagnets, Eq. (|21|l holds approximately 
near a given m with weakly m-dependent scattering co- 
efficients s^'--'-^ for spin-1/2 carriers. The adiabatic spin- 
pumping current is then given by Eq. (|38|) in terms of a 
possibly m-dependent A^^. 

For ferromagnets with a strong spin-orbit coupling, the 
most general form for the adiabatic spin-pumping current 



C. Inhomogeneous magnetization dynamics 

We have so far restricted our attention to spatially- 
uniform magnetization dynamics. Recall, for example, 
Eq. H75|) for spin-pumping-enhanced Gilbert damping in 
F|N heterostructures. Several authors recently discussed 
possibilities to access also nonuniform spin-wave modes, 
^olianski and Brouwcr, 2004 showed that for sufficiently 
large perpendicular-current densities, a single thin-film 
ferromagnet sandwiched between diffuse normal-metal 
contacts becomes unstable with respect to transverse (to 
the layering direction) excitations of finite-wavelength 
spin waves when the source and drain contacts are asym- 
metric, and only for one direction (determined by the 
asymmetry) of the current bias. Their calculation is 
based on the magneto electronic ci r cuit t heory and adi- 
abatic spin pumping. IStiles et all \2004 considered the 
magnetic instability in diffuse N|F|N structures in the 
limit of thicker F layers that can undergo longitudi- 
nal (along the layering direction) as well as transverse 
magnetization dynam ics, allowing to relax t he contact- 
asymmetry condition. lOzvilmaz et "all. 120041 reported an 
experimental study of the current-induced excitations in 
Cu|Co|Cu nanopillars, qualitatively confirming the theo- 
retical ideas. iBrataas et a^J . l2005d : lOzvilmaz et all . l2005l 
investigated current-driven spin-wave instabilities also 
in spin valves, in which they compete with the coher- 
ent macrospin motion. The detailed discussion of the 
current-induced instabilities is beyond the scope of the 
present review. We would nonetheless like to outline how 
the Gilbert damping l(75|) is affected by a weak transverse 
spin-wave excitation in a single magnetic film in contact 
with a diffuse metal. 

Consider a nonuniform transverse excitation of the 
thin-film magnetization in the limit of small amplitudes, 
so that spin waves at different wavelengths do not couple. 
It is then sufficient to study excitations at a single wave- 
length with the magnetization direction deviating from 
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its equilibrium value mo as 

m(r) — mo = 6m cos (q^ ■ p + ujt + (p) 



(125) 



where ip is an arbitrary phase. The spatial-position vec- 
tor r is decomposed here into a coordinate along the 
layering direction, x, and a transverse vector in the in- 
terface plane, p. The wave vector of transverse 
spin waves is parallel to the F|N interface and the am- 
plitude Sm does not depend on x. The derivation of 
the effective Gilbert damping for transverse spin waves 
(jPolianski and Br ouwci. 2QQ£) closely resembles that in 
section IIV.BI for the macrospin dynamics. Let us con- 
sider a normal layer capping one side of the ferromagnet. 
The presence of two normal layers sandwiching a mag- 
netic film thicker than its coherence length Age simply 
doubles the effect. For an unbiased structure, there is 
no charge current or voltage imbalance as long as the 
dynamics are slow on the characteristic spin-relaxation 
time scales, uj <C , which is assumed in the following. 
In order to find the enhanced Gilbert damping, the dif- 
fusion equation for spins in the normal metal l|65() must 
be solved with boundary condition IjtiTI) . In contrast to 
section HV.BI the spin-pumping current now depends on 
the transverse coordinate p. The solution of the diffusion 
equation for the spin accumulation in the normal-metal 
film with thickness L as a function of distance x from the 
F|N interface is [cf. Eq. ^] 

Psix,p)^cosh {x-L)/x'-^f^ /sinh L/X^f^ 



2X 



(cff) 
sd 



f^gjj cos {q^- p + ujt + if) I, , (126) 

where the effective spin-diffusion length for the transverse 
mode is defined by 



A, 



(cff) 
sd 



1 + (Asdq_L) 



(127) 



which reduces to the usual Agd, Eq. (|70|l . for uni- 
form dynamics. By a calculation similar to that of 
section IIV.BI the wave- vector dependence of the en- 
hanced Gilbert damping that generalizes Eq. l(75|l reads 



l|Polianski and Brouweifl2004|) : 



Gcff(q±)-G = 



tanh 



Thick layers with L 3> x'f^^ are thus the best spin sinks 
for a given material composition and spin-wave wave vec- 
tor. 

There is a crossover in the behavior of the enhanced 
Gilbert damping when the wavelength is comparable to 
the spin-diffusion length: In the long- wavelength limit, 
Asw ^ 27rAsd (Asw = 27r/g_L), the result reduces to 
that of a uniformly-precessing ferromagnet; for short- 
wavelength excitations, Asw <C 27rAsd, the damping de- 
pends on the wavelength corresponding to the reduced 



4:TTV 

(128) 



effective spin-diffusion length A^^*^^ « Asw/(27r). In the 
latter regime, numerical estimates for transition-metal 
ferromagnets in contact with simple normal metals, in 
the spirit of Sec. ITyTBl give 



Goff(Asw) — G ' 



1 + [4A/Asw tanh(27rL/Asw)] 



(129) 



where A is the transport mean free path. We thus find 
that a normal metal is always a good spin sink in the limit 
Asw ^ A < i, independently of the spin- relaxation rates, 
in stark contrast to the long- wavelength result l|7t)|) . This 
can be understood referring to the discussion of the dy- 
namic exchange coupling in Sec. IV. Al and noticing that 
the pumping and backfiow reabsorption of spins are sep- 
arated in space by distances of the order or larger than 
A, corresponding to regions of the magnetic layer with 
dynamics that are incoherent upon averaging over vari- 
ous diffuse paths when Asw <C A. The damping of each 
magnetic region is therefore affected by the spin pumping 
with vanishing spin backfiow, rendering the normal metal 
a good spin sink. Consequently, the efficiency of the nor- 
mal metal as a spin sink increases for shorter-wavelength 
spin- wave excitations and longer mean free paths. In par- 
ticular, the Gilbert damping can be enhanced even when 
there is no enhancement of the Gilbert damping for long- 
wavelength excitations. This general conclusion can also 
be extended to magnetic films that are inhomogeneous 
or disordered on length scales shorter than the trans- 
port mean free path in nonmagnetic buffers. In addi- 
tion, the normal metal may become an efficient spin sink 
for both short- and long-wavelength spin waves, when 
the frequency of these excitations becomes larger than 
the normal-metal scattering rate, a regime not explicitly 
treated in this review. (The interfacial spin-pumping cur- 
rent can still be evaluated by the adiabatic formalism, as 
long as the frequency remains much smaller than the fer- 
romagnetic exchange energy.) Finally, we remark that 
quite generally, in the limit when the normal metal be- 
comes effectively an ideal spin sink, the mixing conduc- 
tance that determines the strength of the spin-pumping 
current should not be renormalized (in the sense of the 
Sec. III.BI discussion) . In particular, it is g^^ and not g^^ 
that enters Eq. (|129|l when the denominator on the right- 
hand side is close to unity. In the intermediate spin-sink 
regime, one has to extend the discussion of Sec. III.Bl to 
laterally-inhomogeneous systems. 



D. Electron-electron Interactions 

The appropriate framework for describing metallic 
magnetism, including the id transition metals, is band 
theory (Kiiblcr, 200Q(), treating the electron-electron in- 
teraction in a mean-field approximation. For qualita- 
tive purposes, this comes down to a simple Stoner or 
s ~ d model Hamiltonian having parabolic free-electron 
dispersion for the conduction electrons with parametrized 
masses and exchange splittings. The local-spin-density 
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approximation (LSDA) of density-functional theory is a 
very successful framework for treating itinerant electron 
systems from first principles. Electron-correlation effects 
are taken into account by the exchange-correlation poten- 
tials, be it in the lo cal approximation . This approach has 
been taken by, e.g.,|ziSiSSl],M 

in calculating 

the scattering matrix and the conductance parameters. 

ISimanek . .2003; .Simanek and HeinrichL l2003l raised the 
question of possible enhancement of spin pumping in 
metallic F|N heterostructures by electron-electron inter- 
actions in the normal metals. A potential candidate 
would be Pd as a normal metal with an interaction- 
enhanced magnetic susceptibility. Pd is "nearly" ferro- 
magnetic, causing, e.g., giant moments around magnetic 
impurities. A ferromagnetic film in contact with such a 
material can thus induce magnetic moments renormaliz- 
ing the e xchange p o tentia l felt by electrons at the F|N 
contact. ISimanekL 12003^ 'Sim anek and HeinrichL l2003l 
considered a problem of an ultrathin "^-function" mag- 
netic layer embedded in a nonmagnetic material with 
a large Stoner enhancement of the spin susceptibility. 
Treating the ferromagnetic exchange field as a perturba- 
tion felt by normal-metal electrons, one can obtain a sig- 
nificant enhancement of its effective mean-field profile by 
electron-electron interactions. This in turn can consider- 
ably increase the s pin- mixing cond uctances that govern 
the spin pumping ijSimanekL l2003|) . Such perturbative 
analysis, however, overestimates the effect of electron 
correlations on spin-mixing conductances and thereby 
spin pumping in transition-metal heterostructures. It has 
been explained in Sec. IVII.A.ll that mixing conductance 
g^^ (computed nonperturbatively by density-functional 
theory) of even the thinnest (and more so the thicker) 
magnetic films is close to the normal-metal Sharvin con- 
ductance, regardless of the possible Stoner enhancement 
|Zwicrzvcki et at, 2005). The mixing conductance that 
determines the strength of the spin pumping is thus not 
expected to be correlated with the normal-metal spin sus- 
ceptibility, as already emphasized in Sec. IIV.AI 

In order to understand that a Stoner enhancement 
does not directly affect the spin-pumping enhancement, 
it is convenient to perform the rotating-frame analy- 
sis of Sec. IIII.BI that is also valid in the presence of 
electron-electron interactions, which are invariant to the 
rotating-frame transformation. The precessing ferromag- 
net changes the polarization in the normal metal along 
the axis of rotation by an amount which corresponds to 
the spin-imbalance potential fis — huj that is determined 
by the precession frequency lo, irrespective of the mate- 
rials under consideration. In particular, the steady-state 
spin-accumulation due to the spin pumping does not 
depend on the magnetic susceptibility. Since it is /x^ that 
drives the nonequilibrium spin flow and not the spin den- 
sity (that is affected by the Stoner enhancement), the 
arguments given before in this article remain valid for in- 
teracting systems treated within the mean-field picture. 

The electron-electron effects such as an enhanced 
normal-metal susceptibility are implicitly and nonper- 



turbatively taken into account in self-consistent ab initio 
calculations of heterostructures . First-princip les calcu- 
lations of the scattering matrix ijZwierzvcki et al. 2005), 
fully include the electron-correlation effects which arc dif- 
ficult to capture by a perturbative formalism. Although 
many researchers in the field of solid-state magnetism are 
not familiar with scattering theory, it appears to be the 
most natural language to pursue the study of coupled 
transport and magnetization dynamics. 

Whereas electron-electron effects beyond the mean- 
field LSDA are probably small for transition metals, this 
does not mean that they can be neglected in other sys- 
tems. For example, a suppression of spin pumping by the 
correlation effects in one-dimen sional metals ( "Luttinger 
liquids" ) has been predicted by iBena and" Balentsll200l 
These authors found that the spin pumping by a moving 
fcrromagnct into a Luttinger liquid through a tunnel bar- 
rier is given by the same expression as for noninteracting 
electrons, viz., Eq. (|38|l . but with parameters Al^ and 
A}^ vanishing as a power of temperature at low temper- 
atures with an exponent characteristic for the Luttinger- 
liquid zero-bias anomaly in tunneling density of states. 
Correlation effects might also become important if one 
wishes to describe spin pumping close to the ferromag- 
netic critical temperature, to understand significant devi- 
ations from the macrospin model for the magnetization, 
and to quantitatively model the spin-pumping parameter 
for interfaces with strongly-correlated materials. 



VIII. SUMMARY AND OUTLOOK 

In this review, we presented a coherent picture of the 
nonlocal magnetization dynamics in heterostructures of 
ferromagnets and nonmagnetic conductors. It is based 
on the assumption of semiclassical transport in the bulk 
materials that is valid for diffuse and chaotic systems, 
as well as a separation of time scales of the electronic 
and magnetic degrees of freedom. Interfaces are treated 
as sharp quantum-mechanical boundary conditions for 
electron distribution functions and nonequilibrium trans- 
port. Except for the phenomcnological treatment of spin- 
flip scattering processes, the theory is derived from first 
principles. The main subject in this context was the 
concept of spin pumping due to moving magnetization 
vectors. The magnetization dynamics affected by the 
spin-transfer torque in the presence of an electrical bias 
should be treated on equal footing with the spin pump- 
ing. The crucial material parameter is the spin-mixing 
conductance that can be computed from ab initio elec- 
tronic band structures. 

Several phenomena can be predicted or explained 
based on the basic formalism. One of them is the in- 
creased Gilbert damping of thin magnetic films in good 
electrical contact with normal metals that efficiently dis- 
sipate spin angular momentum. In more complex mag- 
netic structures, we predict an interplay between the spin 
pumping and magnetization torques that is truly nonlo- 
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cal, i.e., depends on the entire spin-coherent volume of 
the sample (determined by the spin-flip diffusion length) . 
Novel collective effects appear when different magnetic 
elements in a spin-coherent circuit or device resonate at 
nearby frequencies. A moving ferromagnet can be used 
as a source that pumps spin currents into normal met- 
als or semiconductors, that leads to a spin accumulation 
determined by the ferromagnetic-resonance frequency. 

Although some basic principles are rather subtle, the 
final formalism is easy to use. ft can often be mapped 
on an equivalent circuit model that is governed by a 
few parameters that can be either fitted to experi- 
ments or computed from first principles. Such calcula- 
tions can be used, e.g., to estimate a nd op timize crit- 
ical magnetization-switching currents ijManschot et all 
l2004j) . We focused attention on quasi-one-dimensional 
configurations and the macrospin model for the magneti- 
zation. Generalization of the formalism to include spin- 
wave excitations in the ferromagnet or inhomogeneous 
spin currents have also been illustrated fSec. lVII.C|l . In- 
tegration of micromagnctic simulations with the trans- 
port equations based on boundary conditions at the in- 
terfaces as described here might be necessary to improve 
the accuracy of first-principles modeling. 

In the review, we had mainly metallic heterostruc- 
tures with transition-metal ferromagnets in mind. But 
since the approach is quite general, we should by no 
means exclude other materials. We already speculated in 
Sec. IVI.51 that the formalism can be used to understand 
the Gilbert damping in magnetic semiconductors. There 
is little doubt that a modeling of t he current-induced 
switching of (Ga,Mn)As observed by IChiba et aR l2004l 
requires the concept of magnetization torque. The 
spin pumping into carbo n nanotubes as investigated by 
iBena and Balentsl. l200i is partly suppressed by correla- 
tion effects. The spin pumping by ferromagn etic super- 
conductors l)Brataas and Tserkovnvail l2004|) is entan- 
gled with Cooper-pair pumping and depends on the spin- 
pairing symmetry of the superconducting state. Spin 
dynamics in heterostructures of high- density magnets in 
contact with doped semiconductors ijBauer et all |-jiMi|[ 
I2OO5) are another promising playground for the formal- 
ism described here. Besides, it is possible to general- 
ize the spin-pumping picture to other symmetry-broken 
heter ostructures with adiabatica l ly- var ying order param- 
eters (jTserkovnvak and Brataa4l2005(l . 

A critical parameter in the nonlocal magnetization dy- 
namics is the spin-flip diffusion length in the normal- 
metal components, that can be of the order of microns 
even at room temperature, see, e.g.. lJedema et aL . 20021 
The nonlocal effects in structures based on transition 
metals could therefore be robust and observable at am- 
bient temperatures. We therefore believe that the for- 
malism discussed here should be useful to understand, 
compute, and design magnetic-device operations. An ex- 
ample is the possibility of "a engineering" based on the 
increase of the Gilbert damping in ultrathin magnetic 
films, e.g., by an order of magnitude by just evaporat- 



ing a few monolayers of Pt on top of it. We also expect 
that with decreasing size of magnetic circuits and de- 
vices, the dynamic coupling discussed in this review will 
become even more relevant. It might lead to effects like 
cross talk between magnetic elements and excess noise. 
On the other hand, proper engineering of phenomena like 
the nonlocal dynamic locking might lead to an increased 
stability against external perturbations as well. 

Theoretical challenges for the future include a proper 
treatment of spin-orbit interactions, the coupling of mag- 
netic degrees of freedom to the lattice, and effects beyond 
the semiclassical regime. 
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